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Abstract 

The Smart Grid is not just about the digitahzation of the Power Grid. In its more visionary acceptation, it 
is a model of energy management in which the users are engaged in producing energy as well as consuming it, 
while having information systems fully aware of the energy demand-response of the network and of dynamically 
varying prices. A natural question is then: to make the Smart Grid a reality will the Distribution Grid have 
to be updated? We assume a positive answer to the question and we consider the lower layers of Medium 
and Low Voltage to be the most affected by the change. In our previous work, we have analyzed samples 
of the Dutch Distribution Grid [58' and we have considered possible evolutions of these using synthetic 
topologies modeled after studies of complex systems in other technological domains |59j . In this paper, we 
take an extra important further step by defining a methodology for evolving any existing physical Power 
Grid to a good Smart Grid model thus laying the foundations for a decision support system for utilities and 
governmental organizations. In doing so, we consider several possible evolution strategies and apply then to 
the Dutch Distribution Grid. We show how more connectivity is beneficial in realizing more efficient and 
reliable networks. Our proposal is topological in nature, and enhanced with economic considerations of the 
costs of such evolutions in terms of cabling expenses and economic benefits of evolving the Grid. 

Keywords: Power Grid, Distribution Grid, Smart Grid, Complex Network Analysis, Decision Support Systems 

1 Introduction 

The Power Grid traces its roots in the late IS*'' century. Then it involved large production facilities and a 
hierarchical Grid composed of High, Medium and Low Voltage to transport and distribute the power to the 
end user. Such infrastructure and its principles are still valid today. As a sector it has traditionally being a 
(state-owned) monopoly. Nowadays there is a strong trend for innovation driven by the need of making the 
infrastructure more reliable, open and to accommodate for new technology, both in term of renewable energy 
generation and digital infrastructure. The technological innovation and the political pressure for a new Power 
Grid often go under the name of Smart Grid, which though has no unique definition |59j . One of the prominent 
aspects that the Smart Grid will enable is the participation of small producers in the energy market. It is already 
possible today, for an end user to produce energy with small-scale production units such as solar panels, small 
wind turbines, and micro-CHP units and sell it back to the unique energy distributor. But more can be achieved 
if the producers would have access to an open energy market where they can auction their over-production and 
buy in a truly competitive market the energy for their needs. Local energy production and distribution will 
change the traditional way power systems have been considered so far: the Low Voltage layer of the Grid will 
go from being a passive end receiving energy from the upper layers, to an active segment of multi-directional 
energy flows thanks for the pervasive distribution of sustainable energy generators. But this trend will inevitably 
affect the physical structure of the Low Voltage Grid, also topologically. The reason for this is that the cost of 
distribution will be of primary importance in enabling or repressing the local market of energy. We also consider 
that the design of the infrastructure will need to take into account the new paradigm of local energy energy 
generation and distribution. 

Assuming the need for an enhancement of the physical Low Voltage Grid, our aim is to realize a decision 
support system to guide the distribution operators, policy makers, and utilities to evaluate scenarios of network 
improvement and to realize distribution Grids that are more efficient and that facilitate (from an economical point 
of view) the delocalized distribution of energy. We propose a process to analyze, design, and adapt distribution 
networks based on statistical models of the Power Grid as a weighted network. A visual representation of the 



process we propose is presented in Figure [T] In the figure, several phases and input are considered to plan the 
evolution of the infrastructure where a local energy exchange is the guiding goal. It starts with a pre-processing 
phase where the input data of the Grid is converted into a graph; the output of this initial phase is a Power Grid 
graph. The following phase consists of the analysis of the topological properties characterizing the graph. The 
output of this phase consists of a set of values representing the metrics related to the Power Grid that influence 
the price of electricity {a and /? metrics in the figure). The process continues with the generation of a network 
model. The number of nodes and edges of this reference model are provided according to the targets for the 
cost-related parameters {a and f3) and the will to invest of the stakeholders. Based on the theoretical model 
identified, the physical network under assessment is then fitted to a topological structure similar to the one of 
the model. Several solutions are provided that differ in the topology and the a and /3 metrics. All these solutions 
are then input to a computer assisted decision support system that presents the benefits/costs of the evolution of 
the network; an expert is involved in the selection of the evolution to be implemented among the most promising 
candidates built by the computer. Once the decision is made, the adaptation of the physical Grid can begin. 

Power Grid Human 




Figure 1: Engineering process for Medium and Low Voltage Grid optimized for prosumer-based energy exchange. 

We covered the initial steps of the process in our previous works. In |58j . we performed topological analysis of 
the Dutch Medium and Low Voltage Power Grid. We developed a set of metrics (a and /3 metrics in the figure) 
based on weighted topological properties to assess the influence on the cost of electricity distribution for a set of 
real Dutch Medium and Low Voltage Grids. In the following [SS], we have evaluated known (reference) models of 
technological networks (such as the Web, the Internet, etc.) to evaluate how they would perform for local energy 
exchange. We found that an increase in connectivity from the typical current value of average degree of two to 
higher values such as four and higher is beneflcial in improving the efficiency and reliability of the network. We 
also found that the small-world model coming from the social sciences with average degree (k) w 4 provides the 
right compromise between performance improvement and the thrift in the realization of this connectivity (i.e., 
cost for cabling) . Although the generation from scratch of a network topology is interesting from a modeling and 
theoretical perspective, it is not common to design a Distribution Grid from scratch. In this paper, we take the a 
practical step in evaluating how to evolve the current Distribution Grid to a more interconnected network taking 
into account the existing network and the physical constraints. We keep the basis of our statistical approach 
to weighted network evolution and apply it to Dutch samples of the Medium and Low Voltage networks. We 
not only assess the benefits in terms of pure topology, but also the costs for cabling the evolved networks, and 
the benefit in terms of costs of electricity distribution. The are several novelties with this approach. First, we 
use statistical graph models as a tool for designing networks and not just for analyzing the existing; second, we 
focus on the Distribution Grid rather than the High Voltage Grid that is usually the target of the analysis in 
the literature; third, we provide a decision support system that works in the large for the network and not for 
individual subcomponents of the Grid; fourth, we provide a tool that can help electricity Distribution companies 
in the design and expansion of Power Grid networks suited for local-scale energy exchange. 

The rest of the paper is organized as follows. We open by analyzing the motivations for a new energy 
landscape and the required changes to the current Grid, Section [2] The basic Graph Theory background is 
presented in Section [Sj Section |4] describes the evolution strategies followed in upgrading the samples of the 
Distribution Grid. The analysis of the results is presented in Section [5] while an overall discussion comparing the 
evolution strategies is provided in Section [6j Section [7] takes into account benefits and costs of evolution of the 
Dutch Grid samples. Section |8] reviews the main approaches to Electrical Grid and System design and evolution, 
while Section [9] provides a conclusion of the paper. [Appendix A| is included to provide extra details about the 
topological metrics that are assessed in the evolution process. 



2 The Need for a New Grid 

From the executive summary of Grid 2030 — A National Vision for Electricity s Second 100 Years:' "Americas 
electric system, the supreme engineering achievement of the 20th century," is aging, inefficient, and congested. 
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and incapable of meeting the future energy needs" ; "Unprecedented levels of risk and uncertainty about future 
conditions in the electric industry have raised concerns about the ability of the system to meet future needs" ; 
"There are several promising technologies on the horizon that could help modernize and expand the Nation's 
electric delivery system, relieve transmission congestion, and address other problems in system planning and 
operations" ; "The revolution in information technologies that has transformed other "network" industries in 
America (e.g., telecommunications) has yet to transform the electric power business"; and finally "It is becoming 
increasingly difficult to site new conventional overhead transmission lines, particularly in urban and suburban 
areas experiencing the greatest load growth" [7^. Here are all the ingredients that already ten years ago set the 
stage for a Grid that needed to evolve and become smarter. The quoted comments refers mainly to the High 
Voltage Power Grid , but the same conclusion can be drawn for the Low Voltage Grid. Actually, the High Voltage 
Grid is already a quite "Smart" with SCADA and EMS systems dealing with the control and communication of 
the system. On the other hand, the Distribution infrastructure, that is usually neglected in the big picture, is 
the layer of the Grid that has less ICT technology in place so far [84] . A key aspect is that the Medium and Low 
Voltage Grid will be the ones that will have to accommodate the renewable sources coming from the distributed 
generation paradigm [35]. We therefore consider that the same problems and challenges that were envisioned 10 
years ago for the Transmission Grid now need to be considered for the Distribution infrastructure. 

Such requirement for the evolution in the Distribution Grid is driven by two inter-related aspects: on the one 
hand, the pressure for unbundling the energy sector, and, on the other hand, the availability and affordability 
of small-scale renewable-based energy production units. Unbundling proposes to get rid of the monopoly system 
that has dominated the energy world so far by providing the possibility of competition in the energy market 
for production, transmission, distribution, and retail. The underlying aim of unbundling is the one of providing 
better services and tariffs for the end-user and promoting innovation and new investment in a traditionally 
sluggish (by definition) sector of the economy dominated by the demand following mindset [3T] . At the extreme 
of the unbundling process is the idea where potentially everybody can produce energy and participate as a seller 
on a free energy market [77j . This last aspect is the linking point with renewable energy production: nowadays 
photovoltaic panels, small-wind turbines are affordable for everybody and often incentivized by governments' 
policies [33] . The step is small to envision a future Grid where everybody can sell the surplus of energy not used 
at home in a market where it is traded as a commodity by software agents embedded in the future generation 
of Smart Meters. In such a context, with many small-scale producers and still without an efficient and cheap 
energy storage technology, a local energy exchange at the neighborhood or municipal level between end-users is 
foreseeable and desirable. Micro-grids increased performance in terms of reduced losses and power quality have 
been successfully tested [HI [M] , but little attention has been devoted to the network topology of these type of 
Grids. 

A future with plenty of prosumers that produce energy and sell or share it at the level of neighborhoods will 
affect the Distribution Grid. The change from a passive-only Grid to a Smart Grid [15] will require to rethink 
the role of the Medium and Low Voltage Grid and the design principles and techniques that have guided its 
development so far. In our study of considering how to evolve or adapt the current Distribution Grid to a Smart 
Grid, we resort to Complex Network Analysis not only to analyze the existing, but also to drive the design of the 
next generation Grid. Complex Network Analysis ( CNA J is a branch of Graph Theory taking its root in the early 
studies of Erdos and Renyi 30 on random graphs and considering statistical structural properties of evolving 
very large graphs. Taking its root in the past. Complex Network Analysis is a relatively young field of research. 
The first systematic studies appeared in the late 1990s [53] [75J [71 S] having the goal of looking at the properties 
of large networks with a complex systems behavior. Afterwards, Complex Network Analysis has been used in 
many diverse fields of knowledge, from biology [H] to chemistry [3^], from linguistics to social sciences [75], 
from telephone call patterns [1] to computer networks [31] and Web [3l [27] to virus spreading [43j |20l [35] to 
logistics [m [38l m] and also inter-banking systems [12] . Men-made infrastructures are especially interesting to 
study under the Complex Network Analysis lenses, especially when they are large scale and grow in a decentralized 
and independent fashion, thus not being the result of a global, but rather of many local autonomous designs. 
The Power Grid is a prominent example. In this work we consider a novel approach both in considering Complex 
Network Analysis tools as a design instrument (i.e., CNA-related metrics are used in finding the most suited 
Medium and Low Voltage Grid for local energy exchange) and in focusing on the Medium and Low Voltage 
layers of the Power Grid. In fact, traditionally, Complex Network Analysis studies applied to the Power Grid 
only evaluate reliability issues and disruption behavior of the Grid when nodes or edges of the High Voltage layer 
are compromised. 

3 Graph Theory Background 

The approach used in this work to model the Power Grid and its evolution is based on Graph Theory and 
Complex Networks. Here we recall the basic definitions that we use throughout the paper and refer to standard 
textbooks such as [IHKII] for a broader introduction. First, we define a graph for the Power Grid [5Dj . 
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Definition 1 (Power Grid graph). A Power Grid graph is a graph G{V,E) such that each element Vi € V is 
either a substation, transformer, or consuming unit of a physical Power Grid. There is an edge eij = {vi, Vj) € E 
between two nodes if there is physical cable connecting directly the elements represented by Vi and Vj . 

Next, we associate weights to the edges representing physical cable properties (e.g., resistance, voltage, sup- 
ported current flow). 

Definition 2 (Weighted Power Grid graph). A Weighted Power Grid graph is a Power Grid graph Gyj{V,E) 

with an additional function / : _E — > K associating a real num,ber to an edge representing the physical property of 
the corresponding cable (e.g., the resistance, expressed in Ohm, of the physical cable). 

A first classification of graphs is expressed in terms of their size. 

Definition 3 (Order and size of a graph). Given the graph G the order is given by N = \V\, while the size is 
given by M = \E\. 

From order and size it is possible to have a global value for the connectivity of the vertexes of the graph, 
known as average node degree . That is <k>= To characterize the relationship between a node and the 
others it is connected to, the following properties provide an indication of the bond between them. 

Definition 4 (Adjacency, neighborhood and degree). // e^.y E is an edge in graph G, then x and y are 
adjacent, or neighboring, vertexes, and the vertexes x and y are incident with the edge e^^y The set of vertexes 
adjacent to a vertex x gV, called the neighborhood ofx, is denoted by T^. The number d{x) = ir^,! is the degree 
ofx. 

A measure of the average 'density' of the graph is given by the clustering coefficient, characterizing the extent 
to which vertexes adjacent to any vertex v are adjacent to each other. 

Definition 5 (Glustering coefficient (CC)). The clustering coefficient 7„ ofVy is 

^ _ l^(r.)l 

(2) 

where \E{Ty)\ is the number of edges in the neighborhood of v and i^^) total number of possible edges in 

Ty. 

This local property of a node can be extended to an entire graph by averaging over all nodes. 

Another important property is how much any two nodes are far apart from each other, in particular the 
minimal distance between them or shortest path. The concepts of path and path length are crucial to understand 
the way two vertexes are connected. 

Definition 6 (Path and path length). A path of G is a subgraph P of the form: 

V{P) = {xo,Xi,...,Xl}, E{P) = {{xo,Xi),{xi,X2),...,{xi-i,Xl)}. 

such that V{P) C V and E{P) C E. The vertexes xq and xi are end- vertexes of P and I = \E{P)\ is the length 
of P. A graph is connected if for any two distinct vertexes Vi,Vj G V there is a finite path from Vi to Vj. 



Definition 7 (Distance). Given a graph G and vertexes Vi and Vj, their distance d{vi,Vj) is the minimal length 
of any — Vj path in the graph. If there is no Vi — Vj path then it is conventionally set to d{vi, Vj) =00. 



Definition 8 (Shortest path). Given a graph G and vertexes Vi and Vj the shortest path is the the path corre- 
sponding to the minimum of to the 5ef {|Pi|, IP2I, • • • , |-Pfe|} containing the lengths of all paths for which Vi and 



Vj are the end-vertexes. 



A global measure for a graph is given by its average distance among any two nodes. 
Definition 9 (Average path length (APL)). Let Vi G V be a vertex in graph G. The average path length for G 

Lai, 

-^--= iV-(iV-l) ^'^^^"^^^ 
where d{vi,Vj) is the finite distance between Vi and Vj and N is the order of G. 
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Definition 10 (Characteristic path length (CPL)). Let Vi £ V be a vertex in graph G, the characteristic path 
length for G, Lcp is defined as the median of dy. where: 



is the mean of the distances connecting Vi to any other vertex Vj in G and N is the order of G. 



To describe the importance of a node with respect to minimal paths in the graph, the concept of betweenness 
helps. Betweenness (sometimes also referred as load) for a given vertex is the number of shortest paths between 
any other nodes that traverse it. 

Definition 11 (Betweenness). The betweenness b{v) of vertex v G V is 



where ast{v) is 1 if the shortest path between vertex s and vertex t goes through vertex v, otherwise and Ugt is 
the number of shortest paths between vertex s and vertex t. 

Looking at large graphs, one is usually interested in global statistical measures rather than the properties of a 
specific node. A typical example is the node degree, where one measures the node degree probability distribution. 

Definition 12 (Node degree distribution). Consider the degree k of a node in a graph as a random variable. 
The function 

Nk = {veG: d{v) = k} 
is called probability node degree distribution. 

The shape of the distribution is a salient characteristic of the network. For the Power Grid, the shape 
is typically either exponential or a Power-law [71 O 1581 166j . More precisely, an exponential node degree (fc) 
distribution has a fast decay in the probability of having nodes with relative high node degree. The relation: 

P(fc) = ae'^'' 

follows, where a and /3 are parameters of the specific network considered. On the contrary, a Power-law distri- 
bution has a slower decay with higher probability of having nodes with high node degree. It is expressed by the 
relation: 

P{k) = ak-'^ 

where a and 7 are parameters of the specific network considered. We remark that the graphs considered in the 
Power Grid domain are usually large, although finite, in terms of order and size thus providing limited and finite 
probability distributions. 

A Graph can also be represented as a matrix, typically an adjacency matrix. 

Definition 13 (Adjacency matrix). The adjacency matrix A — A{G) — {ai,j) of a graph G of order N is the 
N X N matrix given by 

^ _ 1 1 if {vt,Vj) e E, 
1 otherwise. 

We have now provided the basic definitions needed to present the modeling tools for the Power Grid evolutions. 



4 Evolving the Current Power Grid 

We start by taking samples of Low Voltage and Medium Voltage from the Dutch Power Grid. Tables [T] and [3] 
summarize the main facts of these samples (described in greater detail in [5S]). The first column of each table 
represents the identifier of the sample, column two and three provide the order and size of the sample; these 
two values are used to compute the average degree that is shown in the fourth column. Column five and six 
give an impression on the the effort to reach one node from any other in the network through the average 
path length and the characteristic path length. A measure of local clustering is given in column seven with the 
clustering coefficient metric. One notes a low average node degree (k) « 2 both for the Low Voltage networks and 
the Medium Voltage networks. Besides the difference in the order and size between the two types of networks 
(generally the Medium Voltage network samples are bigger), one sees that the Low Voltage samples have a mostly 
a null clustering coefficient, while the Medium Voltage networks present a small, but at least significant, value. 
This difference is explained in the different topology and purpose of the networks: a radial structure with no 
clustering to distribute electricity to the end-users (Low Voltage) and a more meshed structure for the Medium 
Voltage that shows small clustering values. 
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Present study 


Random Graph 


iU 


Order 


oize 


Avg. 

d 






7 




r^r>T 


7 


1 


17 


18 


2.118 


3.398 


3.313 


0.00000 


1.427 


1.688 


0.13726 


2 


15 


16 


2.133 


3.086 


3.000 


0.00000 


2.319 


2.358 


0.00000 


3 


24 


23 


2.087 


4.499 


4.228 


0.00000 


3.127 


3.091 


0.05508 


4 


30 


29 


1.933 


4.545 


4.449 


0.00000 


1.860 


2.242 


0.05778 


5 


188 


191 


2.032 


17.726 


17.878 


0.00000 


3.846 


4.345 


0.00532 


6 


10 


9 


1.800 


2.423 


2.223 


0.00000 


0.978 


1.167 


0.26667 


7 


63 


62 


1.968 


5.204 


5.404 


0.00000 


2.514 


2.904 


0.03175 


8 


28 


27 


1.929 


4.784 


5.000 


0.00000 


2.553 


2.945 


0.04762 


9 


133 


140 


2.105 


11.543 


11.366 


0.01112 


3.702 


4.172 


0.01482 


10 


124 


138 


2.226 


8.053 


7.070 


0.00869 


3.010 


3.540 


0.02914 


11 


31 


30 


1.935 


4.353 


4.357 


0.00000 


1.590 


1.969 


0.07475 



Table 1: Low Voltage samples from the northern Netherlands Power Grid compared with Random graphs of the 
same size. 



Network 
sample 


Order 


Size 


Avg. be- 

tweenness 


Avg. 

betw/order 


Coeff. 
varia- 
tion 


1 


17 


18 


32.933 


0.074 


0.773 


2 


15 


16 


25.231 


0.053 


0.887 


3 


24 


23 


70.286 


0.295 


0.643 


4 


30 


29 


81.167 


0.309 


1.153 


5 


188 


191 


2928.227 


13.494 


1.207 


6 


10 


9 


9.000 


0.047 


1.291 


7 


63 


62 


255.016 


0.288 


2.091 


8 


28 


27 


102.143 


0.279 


1.301 


9 


133 


140 


1355.953 


6.220 


1.534 


10 


124 


138 


771.911 


3.840 


1.351 


11 


31 


30 


139.677 


0.691 


1.265 



Table 2: Betweenness for Dutch Low Voltage samples. 



Next, we consider evolutions starting from the Dutch samples, that is, adding cables according to several 
strategies of network growth. We break the evolutions into four groups of edge growth: increments of 25%, 50%, 
75%, and 100%. The choice of stopping at 100% is performed based on the results of ^59j, where it is shown that 
an average node degree of 4 has the right balance of improved network qualities and costs of network evolution. 
We consider several strategies for evolving the graph by adding more links, namely: 

• Assortativity. A network is assortative if nodes having similar characteristics or properties are connected 
one another [56j. The property we consider is that of node degree and then take two strategies: 

— High degree nodes are connected one another. The process starts with considering a set of nodes 
with the highest equal node degree and connect them together. The process goes on considering the 
next set of nodes with equal high degree in the order of rank and so on. 

— Low degree nodes are connected one another. The process goes on as for the the high degree strategy, 
but node are linked starting from the couples with lowest degree. 

• Dissortativity is the opposite of assortativity, that is, a network is dissortative if nodes having different 
characteristics or properties are connected together. Following this strategy, nodes with highest node degree 
are linked to nodes with lowest node degree. 

• Triangle closure is based on the principle of increasing the clustering coefficient of the network. At each 
step, a node is selected at random and for each pair of its neighbors an edge is added between them, if not 
already present. 

• Least distance gives priority to the connection of nodes that are geographically closer to each other. This 
strategy can minimize the costs of cabling since such costs are directly proportional to the length of cables. 
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Random Graph 


ILJ 


Order 


oize 


Avg. 
d 






7 




r^r>T 


7 


1 


444 


486 


2.189 


11.033 


10.858 


0.00537 


5.547 


6.163 


0.00333 


2 


472 


506 


2.144 


17.095 


17.174 


0.01360 


5.039 


5.700 


0.00106 


3 


238 


245 


2.059 


11.715 


11.580 


0.00000 


3.558 


4.234 


0.00595 


4 


263 


288 


2.190 


12.775 


12.311 


0.01118 


5.046 


5.368 


0.01080 


5 


217 


229 


2.111 


10.321 


10.241 


0.00140 


4.894 


5.391 


0.00121 


6 


191 


207 


2.168 


9.288 


8.990 


0.00296 


4.616 


5.079 


0.00225 


7 


884 


1059 


2.396 


9.817 


9.527 


0.00494 


5.440 


6.010 


r\ r\r\'i 

0.00170 


8 


366 


382 


2.087 


15.113 


14.546 


0.00000 


4.691 


5.249 


0.00405 


9 


218 


232 


2.128 


10.850 


10.915 


0.00000 


5.454 


5.856 


0.00539 


10 


201 


204 


2.030 


15.742 


15.257 


0.00166 


4.898 


5.503 


0.00491 


11 


202 


213 


2.109 


13.504 


12.891 


0.00140 


4.801 


5.217 


0.08750 


12 


464 


499 


2.151 


13.144 


12.703 


0.00036 


4.718 


5.390 


0.00209 



Table 3: Medium Voltage samples from the northern Netherlands Power Grid compared with Random graphs of 
the same size. 



Network 


Order 


Size 


Avg. be- 


Avg. 


Coeff. 


sample 






tweenness 


betw/order 


varia- 












tion 


1 


444 


486 


4329.054 


9.750 


2.050 


2 


472 


506 


5087.728 


10.779 


1.704 


3 


238 


245 


1910.757 


8.028 


1.566 


4 


263 


288 


1237.711 


4.706 


1.517 


5 


217 


229 


3169.571 


14.606 


1.743 


6 


191 


207 


3870.640 


20.265 


1.432 


7 


884 


1059 


7755.542 


8.773 


2.875 


8 


366 


382 


5136.520 


14.034 


1.691 


9 


218 


232 


1244.663 


5.709 


1.544 


10 


201 


204 


3613.691 


17.979 


1.173 


11 


202 


213 


2690.183 


13.318 


1.331 


12 


464 


499 


3424.602 


7.381 


1.687 



Table 4: Betweenness for Dutch Medium Voltage samples. 



• Random is based on the random selection of nodes to attach edges. At each step of the growth process, 
a pair of distinct nodes are randomly selected and an edge between them is added. 

For every strategy, if two nodes already have an edge that connects them the edge is not added and the evolution 
strategy continues. In fact, in the graph models we only allow a single edge between a pair of nodes, if not 
already present. 

5 Prom the Current Distribution Infrastructure to the Smart Grid 

We adapt current physical networks according the strategies described in Section [4] and we analyze the obtained 
graphs according to a set of metrics that provide a view of efficiency of the whole network and its adequacy 
for local energy distribution. Such metrics consider the path length properties of the graph, the presence of 
cliques at local scale (i.e., clustering coefficient), the presence of critical nodes that manage the majority of paths 
(i.e., betweenness). Robustness of the network to failures concerning its connectivity (RobN) is evaluated by 
computing the the order of the maximal connected component (MCC) of the graph when nodes (20% of the 
initial order of the graph) are removed randomly or in a targeted way by focusing on the nodes with higher 
degree first, then taking the average between the two values to have an overall estimation of the disruption. The 
last metric we analyze deals with the redundancy of paths (APLiQth)] which takes into account the increase 
in the average path length when the 10*'' shortest path is computed, therefore this metric provides a measure 
of the additional effort required to benefit of alternative paths than the optimal one. For a more detailed 
dissertation over the metric analyzed we refer to Appendix A To study the effects of the strategy, we implement 



the strategies described in a software based on the JAVA graph library JGraphT (http://www.jgrapht.org/). 
The same library suite has been used to compute the metrics just described. The only metric computed with a 
different software library is the 'betweennees' one. For this computation the Stanford Network Analysis Project 
(SNAP - http://snap.staiiford.edu/) software library has been used since it leverages on the the algorithm 
developed by Brandos [T3] with optimal performance. To perform the generation and computation of the metrics 
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we used a PC with Intel Core2 Quad CPU Q9400 2.66GHz with 4GB RAM. The Operating system is based 
on the Linux kernel 2.6.32 with a 4.4.3 GCC compiler and JAVA framework 1.6. The versions of JGraphT and 
SNAP software libraries used are respectively vlO. 10.01 and vO.8.1. Next we present the results of the generation 
and metrics evaluation. 

5.1 Evolution of Medium Voltage Distribution Grids 

We start by considering the Medium Voltage Grid and apply the evolution strategies presented in the previous 
section one at the time. 



Assortative high node degree evolution 

The results for evolution according of the assortative strategy is summarized in Tabic [5] For each sample of the 
Medium Voltage Grid (column one) the values for the main topological quantities: order and size in columns 
three and four, respectively; average node degree in the fifth column; the characteristic path length is reported 
in column six; the clustering coefficient follows in column seven; robustness is shown in eight column; the cost in 
term of redundant path length closes the data series (column nine). For the path length, we note that the increase 
in the connectivity is extremely beneficial. This is notable already with the addition of 25% of links, where on 
average about 50% of the path length is reduced. The improvement is then minor when more link are added, in 
fact, when the number of links is doubled the reduction of the characteristic path length is around 60%, on average. 
The increment in the transition +75%/+100% of the number of links brings just a gain in the reduction of the 
characteristic path length of 3%. This tendency of saturation in the reduction of the path length is explained by 
the small-world phenomenon that arises when a sufficient connectivity threshold is reached. The improvement for 
the clustering coefficient is significant, reaching values that are even two orders of magnitude higher. In general, 
the clustering coefficient has a linear improvement in the evolution steps considered, some of the samples show 
however a behavior that is similar to a logistic shape[^ as shown in Figure [3j The improvement in robustness 
tends to double in values. An exception is sample #7 that had already a good initial value for robustness. The 
addition of edges according to this strategy is mainly beneficial in the random node attack case, while the fraction 
of the metric that considers targeted attacks against the most connected nodes is marginally affected. In fact, 
the ranking of the nodes is almost untouched reinforcing the node degree of those nodes that had already an 
high degree at the beginning of the evolution process. As explained by Newman [55], assortative networks that 
link high degree nodes have a sort of redundancy in the main cluster. The redundancy of new connections is 
not particularly helpful in improving the resilience of the network since the nodes with high connectivity already 
form a cluster structure. The addition of edges between the high degree nodes does not increase the number 
of new potential target nodes. This is empirically evident by the results of the robustness metric presented in 
Figures [2a| and |2b[ showing the evolution of the robustness metric in the random and targeted node removal 
situations, respectively. One sees that the increase in connectivity is beneficial to contrast random attacks and 
the first step of link addition is the most beneficial. Then the disruption of the network does not benefit from 
the additional connectivity anymore. The lack of benefit from the additional connectivity is emphasized in the 



resilience against targeted attacks (cf. Figure 2a) where almost all samples do not have improvements in their 
performance with just two exceptions. 

Considering the values for the redundant path robustness, the same considerations on the characteristic path 
length apply. There is a reduction of more than 50% already when the networks are evolved with a 25% increase 
in the size. An exception is sample #2 that shows a consistent improvement also in the later stages of evolution, 
especially between step two and step three where the average redundant path length from about ten, declines to 
a value slightly higher than six. 

Table |6] shows for each sample of the Medium Voltage Grid (column one) the values of metrics related to 
betweenness. Columns three and four show order and size, then average betweenness is provided in column 
five, while a value of average betweenness normalized by the order of the graph is presented in column six as it 
provides a normalized value. The statistical coefficient of variation which is shown in the seventh column. 

The assortative strategy involving the nodes with highest node degree is beneficial in reducing the average 
betweenness of all the Medium Voltage samples. Already from the first step of the network evolution the 
betweenness on average reduces to almost 50% of the original value. Some samples (i.e., samples #5 and #6) 
reach even higher reduction of up to 80%. The same trend is followed by the normalized value of average 
betweenness divided by the order of the network. Already in the first step no sample exceeds seven for the 
betweenness to order ratio. Both metrics improve of up to 60% (on average) in the last step of the evolution. 
Considering the variability of betweenness in each of the samples we note a general increase. In particular, only 
four of the twelve samples show a decrease when the connectivity in the network is double of the original size of 
the network. Such behavior is due to a slower decrease of the standard deviation of betweenness compared to its 
average. This strategy of adding the connection between the nodes that already have the highest connectivity, that 



' http : //en . ulklpedla . org/wiki/Logistic_curve 
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Figure 4: Evolution of robustness sub-metrics. 



are usually the nodes that also have highest betweenness, reduces the number of "bottleneck" nodes. However, 
the strategy of adding links is not helpful for substantially reducing the variability of betweenness. Evidence to 
this claim comes from the fact that the median of the betweenness in each of the samples at every state of the 
evolution process is zero, that is, the majority of nodes are terminal nodes that are not involved in any path 
between other nodes. 



Assortative low node degree evolution 

Table [t] contains for each sample of the Medium Voltage Grid (column one) the values for main topological 
quantities: order and size in columns three and four, respectively; average node degree in the fifth column; the 
characteristic path length is reported in column six; the clustering coefficient follows in column seven; robustness 
is shown in eight column; and the cost in term of redundant path length closes the data series (column nine). 

The increase in the connectivity with the assortative low node degree strategy is beneficial for the path 
length already with the addition of 25% of links reducing about 45% the initial value of the original graph. The 
improvement is lower when more and more links are added. In fact, when the number of link is doubled compared 
to the original size, the reduction of the characteristic path length is around 47% on average. In particular, the 
improvement comparing all the evolution steps is around 0.9, which shows a limit in the benefits achievable by 
the increased connectivity. This effect is explained by the small- world phenomenon that arises when a sufficient 
connectivity threshold is obtained. Once this threshold is obtained, the subsequent addition of edges following 
the same strategy have a reduced effect on the improvement of the property. Following this evolution strategy, 
we see that a saturation effect arises after the first two steps of the growth process since the characteristic 
path length then has no more significant improvement. The improvement for the clustering coefficient increases 
substantially compared to the initial values. It reaches values that are even three order of magnitude higher than 
the initial situation. In general, the improvement in the clustering coefficient tends to have a logistic trend in the 
evolution step considered, as shown in Figure [5] The improvement in robustness is in general three times higher 
compared to the initial value for the same samples. Three out of the twelve samples have values for robustness 
higher than 0.6. Even in this evolution scenario, sample #7 poses an exception whose improvement are limited 
compared to the other samples. The addition of edges according to this strategy is beneficial in contrasting the 
effects of both random attacks and targeted attacks against the most connected nodes. This evolution strategy 
tends to change the hierarchy (in terms of node degree) of the nodes, adding more connections between the nodes 
that are less connected. Therefore, these nodes assume more importance in terms of node degree compared to 
the initial situation, giving more homogeneity in the degree. We see in this strategy how the new connectivity 
is particularly beneficial in contrasting targeted attacks. The new established connections tend to bond nodes 
with low degree (that anyway become small hubs of the network) thus improving the resilience of the network. 



This is empirically evident by the results of the robustness metric (Figures 4a and 4b ) that show the evolution 
of the robustness metric in the random and targeted node removal situations, respectively. One sees that the 
increase in connectivity is beneficial to contrast random attacks that reach a plateau around 0.6-0.7 after the 
second evolution step and more beneficial against the attacks that target high degree nodes: the majority of the 
samples experiences an improvement in the metric about one order of magnitude. 

Considering the values for the redundant path robustness the same considerations done for the characteristic 
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Sample ID 


Network 
type 


Order 


Size 


Avg. 
deg. 


CPL 


CC 


Removal 
robust- 
ness 

(.RobN) 


Redundancy 
cost 


original order +25% 


1 


MV 


444 


607 


2.734 


4.784 


0.03797 


0.330 


7.871 


2 


MV 


472 


632 


2.678 


9.203 


0.03242 


0.247 


11.476 


3 


MV 


238 


306 


2.571 


6.243 


0.02486 


0.295 


9.549 


4 


MV 


263 


360 


2.738 


5.973 


0.03196 


0.346 


8.388 


5 


MV 


217 


286 


2.636 


4.704 


0.03072 


0.329 


7.855 


6 


MV 


191 


258 


2.702 


4.889 


0.02907 


0.341 


7.455 


7 


MV 


884 


1323 


2.993 


7.312 


0.02723 


0.318 


9.094 


8 


MV 


366 


477 


2.607 


6.532 


0.02451 


0.303 


9.120 


9 


MV 


218 


290 


2.661 


4.806 


0.03351 


0.346 


7.452 


10 


MV 


201 


255 


2.537 


6.360 


0.02425 


0.294 


8.943 


11 


MV 


202 


266 


2.634 


5.669 


0.03563 


0.314 


8.449 


12 


MV 


464 


623 


2.685 


4.871 


0.01950 


0.319 


7.884 


original order +50% 


1 


MV 


444 


729 


3.284 


4.734 


0.07307 


0.346 


6.305 


2 


MV 


472 


759 


3.216 


8.304 


0.05193 


0.264 


9.975 


3 


MV 


238 


367 


3.084 


5.405 


0.05780 


0.320 


6.995 


4 


MV 


263 


432 


3.285 


4.676 


0.07645 


0.367 


6.000 


5 


MV 


217 


343 


3.161 


4.486 


0.06157 


0.363 


6.228 


6 


MV 


191 


310 


3.246 


4.805 


0.06290 


0.349 


6.425 


7 


MV 


884 


1588 


3.593 


4.496 


0.06287 


0.,363 


6.848 


8 


MV 


366 


573 


3.131 


5.638 


0.05272 


0.324 


7.972 


9 


MV 


218 


348 


3.193 


4.749 


0.06138 


0.354 


6.732 


10 


MV 


201 


306 


3.045 


6.145 


0.05243 


0.296 


8.195 


11 


MV 


202 


319 


3.158 


5.502 


0.06155 


0.318 


8.171 


12 


MV 


464 


748 


3.224 


4.847 


0.04752 


0.343 


6.728 


original order +75% 


1 


MV 


444 


850 


3.829 


4.705 


0.09481 


0.346 


6.267 


2 


MV 


472 


885 


3.750 


4.327 


0.07763 


0.352 


6.311 


3 


MV 


238 


428 


3.597 


5.,344 


0.07648 


0.,320 


6.980 


4 


MV 


263 


504 


3.833 


1.622 


0.10929 


0.382 


6.018 


5 


AiV 


217 


100 


:-!.()iSV 


1.119 


0.08276 


0.387 


r,.s:',:.i 


6 


MV 


191 


362 


3.791 


4.532 


0.08257 


0.349 


6.049 


7 


MV 


884 


1853 


4.192 


4.384 


0.10213 


0.378 


5.738 


8 


MV 


366 


668 


3.650 


5.573 


0.07256 


0.,326 


7.627 


9 


MV 


218 


406 


3.725 


4.712 


0.07867 


0.356 


6.740 


10 


MV 


201 


357 


3.552 


5.9,35 


0.07022 


0.293 


8.019 


11 


MV 


202 


372 


3.683 


5.435 


0.08035 


0.314 


7.380 


12 


MV 


464 


873 


3.763 


4.809 


0.06792 


0.338 


6.641 


original order +100% 


1 


MV 


444 


972 


4.378 


4.646 


0.10837 


0.346 


6.195 


2 


MV 


472 


1012 


4.288 


4.297 


0.10749 


0.378 


5.684 


3 


MV 


238 


490 


4.118 


5.135 


0.09390 


0.321 


6.773 


4 


MV 


263 


576 


4.380 


4.595 


0.12805 


0.436 


5.476 


5 


MV 


217 


458 


4.221 


4.347 


0.10451 


0.397 


5.580 


6 


MV 


191 


414 


4.335 


4.268 


0.10310 


0.348 


5.947 


7 


MV 


884 


2118 


4.792 


4.375 


0.12869 


0.,366 


5.410 


8 


MV 


366 


764 


4.175 


5.538 


0.08809 


0.320 


7.673 


9 


MV 


218 


464 


4.257 


4.664 


0.09572 


0.355 


6.227 


10 


MV 


201 


408 


4.060 


5.610 


0.08821 


0.294 


6.741 


11 


MV 


202 


426 


4.218 


4.391 


0.13484 


0.366 


6.249 


12 


MV 


464 


998 


4.302 


4.798 


0.08111 


0.339 


6.488 



Table 5: Metrics for assortative high node degree strategy Medium Vohage samples evolution. 
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Sample ID 


Network 
type 


Order 


Size 


Avg. be- 
tweenness 


Avg. 

betw/order 


Coeff. 
vEiria- 
tion 


original order +25% 


1 


MV 


444 


607 


1696.993 


3.822 


4.047 


2 


MV 


472 


632 


3290.807 


6.972 


1.651 


3 


MV 


238 


306 


1048.901 


4.407 


1.803 


4 


MV 


263 


360 


1236.901 


4.703 


1.909 


5 


MV 


217 


286 


857.71 


3.953 


2.937 


6 


MV 


191 


258 


750.095 


3.927 


2.181 


7 


MV 


884 


1323 


5826.819 


6.591 


2.067 


8 


MV 


366 


477 


1989.602 


5.436 


1.901 


9 


MV 


218 


290 


905.139 


4.152 


2.4 


10 


MV 


201 


255 


1144.505 


5.694 


1.524 


11 


MV 


202 


266 


1026.579 


5.082 


1.756 


12 


.MV 


464 


623 


2039.336 


4.395 


4.074 


original order- +50% 


1 


MV 


444 


729 


1673.842 


3.77 


2.707 


2 


MV 


472 


759 


2827.063 


5.99 


1.707 


3 


MV 


238 


367 


886.658 


3.725 


1.67 


4 


MV 


263 


432 


886.347 


3.37 


2.294 


5 


MV 


217 


343 


789.065 


3.636 


2.236 


6 


MV 


191 


310 


732.148 


3.833 


1.643 


7 


MV 


884 


1588 


3243.961 


3.67 


6.472 


8 


MV 


366 


573 


1748.53 


4.777 


1.85 


9 


MV 


218 


348 


884 


4.055 


1.817 


10 


MV 


201 


306 


1097.948 


5.462 


1.221 


11 


MV 


202 


319 


965.401 


4.779 


1.724 


12 


MV 


464 


748 


2020.101 


4.354 


2.556 


original order +75% 


1 


MV 


444 


850 


1659.788 


3.738 


2.195 


2 


MV 


472 


885 


1577.847 


3.343 


4.268 


3 


MV 


238 


428 


867.261 


3.644 


1.499 


4 


MV 


263 


501 


870.669 


3.311 


1.825 


5 


AIV 


217 


100 


775.35."') 


3.573 


1.782 


6 


MV 


191 


362 


688.18 


3.603 


1.624 


7 


MV 


884 


1853 


3111.479 


3.52 


4.838 


8 


MV 


366 


668 


1734.622 


4.739 


1.674 


9 


MV 


218 


406 


869.273 


3.987 


1.571 


10 


MV 


201 


357 


1061.907 


5.283 


1.135 


11 


MV 


202 


372 


940.802 


4.657 


1.683 


12 


MV 


464 


873 


1998.422 


4.307 


2.111 


original order +100% 


1 


MV 


444 


972 


1639.753 


3.693 


1.964 


2 


MV 


472 


1012 


1544.805 


3.273 


3.358 


3 


MV 


238 


490 


827.586 


3.477 


1.618 


4 


MV 


263 


576 


858.372 


3.264 


1.625 


5 


MV 


217 


458 


754.514 


3.477 


1.615 


6 


MV 


191 


414 


657.291 


3.441 


1.592 


7 


MV 


884 


2118 


3099.314 


3.506 


3.815 


8 


MV 


366 


764 


1719.62 


4.698 


1.606 


9 


AlV 


218 


16 1 


857.522 


3.931 


1.458 


10 


Ai\" 


201 


108 


1015.278 


5.051 


1.11 


11 


MV 


202 


426 


935.381 


4.631 


1.679 


12 


MV 


464 


998 


1989.684 


4.288 


1.875 



Table 6: Betweenness for assortative high node degree strategy Medium Voltage samples evolution. 
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path length apply. There is already a reduction of 45% compared to the initial value of the samples already when 
the networks are evolved with a 25% increase in the size of the original graph. 
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Figure 5: Evolution of the clustering coefficient metric. 



Table [S] shows for each sample of the Medium Voltage Grid (column one) the values of metrics related to 
betweenness. In addition to order and size (columns three and four), average betweenness is provided in column 
five, while a value of average betweenness normalized by the order of the graph is computed in column six. A 
measure of the statistical variation of betweenness is the coefficient of variation which is shown in the seventh 
column. 

Considering betweenness, the assortative strategy involving the nodes with lowest node degree is beneficial in 
reducing the average betweenness of all the Medium Voltage samples. Already from the first step of the network 
evolution the betweenness on average reduces more than 38% of the original value. However samples #4 and 
#9 have in the first step a slight increment in the average betweenness. The same reduction trend is followed 
by the normalized value of average betweenness divided by the order of the network; already in the first step no 
sample exceeds 7.5 for the betweenness to order ratio. The benefits for both these metrics improve up to 41% 
(on average) in the last step of the evolution. Considering the variability of betweenness in each of the samples, 
we note a general increase for this metric in the first two steps of the evolution, while the tendency is inverted in 
the last two steps. In particular, only sample #8 shows a slight increase when the connectivity in the network is 
double the original size of the network. Such behavior is due to an initial increase in the standard deviation of 
betweenness when only few edges are added, then after the second step of the evolution the tendency is inverted 
and the standard deviation decreases. This strategy of adding the connection, in fact, provides more connections 
between the nodes that have small connectivity. These are usually the nodes at the periphery of the network 
and do not have paths traversing them. A proof is that the median of the betweenness for each sample in the 
first two stages of evolution is zero, then in the later stages some samples (four out of fourteen) have a non-zero 
median that is a sign that all nodes are more evenly involved in the paths between other nodes. 



Triangle closure evolution 

Table |9] contains for each sample of the Medium Voltage Grid (column one) the values for the main topological 
quantities: order and size in columns three and four, respectively; average node degree in the fifth column; the 
characteristic path length is reported in column six; the clustering coefficient follows in column seven; robustness 
is shown in eight column; the cost in term of redundant path length closes the data series (column nine). 

The primary focus of the triangle closure strategy is to improve the clustering coefficient of the network and, 
as a side effect, the other metrics benefit from the additional 'local' links (the links added connect the neighbors 
of a node). The increase in the connectivity with such a strategy is beneficial for the path length: the addition of 
25% of links reduces the path length around 12% on average, the trend is slightly sub-linear and the final step of 
the evolution provides a shrinking of the characteristic path length measure to about 38%. The best improvement 
is naturally obtained for the clustering coefficient. The average of the clustering coefficient reaches almost 0.5 
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Sample ID 


Network 
type 


Order 


Size 


Avg. 
deg. 


CPL 


CC 


Removal 
robust- 
ness 

(.RobN) 


Redundancy 
cost 


original order +25% 


1 


MV 


444 


607 


2.734 


8.037 


0.05039 


0.,324 


10.170 


2 


MV 


472 


632 


2.678 


6.146 


0.04055 


0.357 


9.194 


3 


MV 


238 


306 


2.571 


6.171 


0.02847 


0.331 


9.713 


4 


MV 


263 


360 


2.738 


6.115 


0.04448 


0.355 


8.514 


5 


MV 


217 


286 


2.636 


5.667 


0.02814 


0.373 


8.608 


6 


MV 


191 


258 


2.702 


5.379 


0.03156 


0.365 


8.331 


7 


MV 


884 


1323 


2.993 


7.439 


0.05290 


0.404 


8.827 


8 


MV 


366 


477 


2.607 


7.130 


0.03298 


0.337 


10.249 


9 


MV 


218 


290 


2.661 


6.661 


0.04338 


0.351 


10.134 


10 


MV 


201 


255 


2.537 


6.650 


0.04379 


0.310 


9.466 


11 


MV 


202 


266 


2.634 


6.774 


0.04825 


0.310 


9.734 


12 


MV 


464 


623 


2.685 


7.254 


0.04071 


0.352 


9.779 


original order +50% 


1 


MV 


444 


729 


3.284 


7.860 


0.08023 


0.374 


9.941 


2 


MV 


472 


759 


3.216 


6.142 


0.10215 


0.420 


7.791 


3 


MV 


238 


367 


3.084 


6.105 


0.08438 


0.442 


8.317 


4 


MV 


263 


432 


3.285 


6.103 


0.09884 


0.431 


7.649 


5 


MV 


217 


343 


3.161 


5.620 


0.09096 


0.422 


7.236 


6 


MV 


191 


310 


3.246 


5.332 


0.09208 


0.430 


7.121 


7 


MV 


884 


1588 


3.593 


7.376 


0.07033 


0.408 


8.576 


8 


MV 


366 


573 


3.131 


7.068 


0.08246 


0.358 


9.747 


9 


MV 


218 


348 


3.193 


6.569 


0.08068 


0.424 


10.086 


10 


MV 


201 


306 


3.045 


6.565 


0.08028 


0.359 


8.689 


11 


MV 


202 


319 


3.158 


6.714 


0.08388 


0.353 


8.946 


12 


MV 


464 


748 


3.224 


7.211 


0.08205 


0.366 


9.542 


original order +75% 


1 


MV 


444 


850 


3.829 


7.795 


0.09130 


0.373 


10.521 


2 


MV 


472 


885 


3.750 


6.125 


0.12608 


0.634 


7.455 


3 


MV 


238 


428 


3.597 


6.046 


0.10408 


0.507 


8.178 


4 


MV 


263 


504 


3.833 


6.069 


0.11968 


0.636 


7.619 


5 


AiV 


217 


100 


3.()8V 


r>.5!)3 


0.11971 


0.588 


(i.77() 


6 


MV 


191 


362 


3.791 


5.279 


0.11311 


0.501 


7.030 


7 


MV 


884 


1853 


4.192 


7.328 


0.07592 


0.405 


8.261 


8 


MV 


366 


668 


3.650 


7.049 


0.10329 


0.405 


9.180 


9 


MV 


218 


406 


3.725 


6.498 


0.09285 


0.425 


8.571 


10 


MV 


201 


357 


3.552 


6.4,30 


0.09449 


0.363 


8.551 


11 


MV 


202 


372 


3.683 


6.654 


0.09751 


0.363 


8.521 


12 


MV 


464 


873 


3.763 


7.195 


0.09838 


0.364 


9.375 


original order +100% 


1 


MV 


444 


972 


4.378 


7.761 


0.09597 


0.371 


9.426 


2 


MV 


472 


1012 


4.288 


6.113 


0.13987 


0.633 


6.909 


3 


MV 


238 


490 


4.118 


5.977 


0.11510 


0.503 


7.432 


4 


MV 


263 


576 


4.380 


5.966 


0.12978 


0.656 


7.346 


5 


MV 


217 


458 


4.221 


5.560 


0.13315 


0.635 


7.478 


6 


MV 


191 


414 


4.335 


5.200 


0.12489 


0.609 


7.037 


7 


MV 


884 


2118 


4.792 


7.258 


0.07806 


0.407 


8.759 


8 


MV 


366 


764 


4.175 


6.758 


0.11224 


0.407 


8.902 


9 


MV 


218 


464 


4.257 


6.394 


0.10102 


0.426 


8.791 


10 


MV 


201 


408 


4.060 


6.340 


0.10233 


0.361 


8.074 


11 


MV 


202 


426 


4.218 


6.604 


0.10718 


0.354 


9.461 


12 


MV 


464 


998 


4.302 


7.188 


0.10677 


0.364 


8.800 



Table 7: Metrics for assortative low node degree strategy Medium Voltage samples evolution. 
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Sample ID 


Network 


Order 


Size 


Avg. be- 


Avg. 


Coeff. 




type 






tweenness 


betw/order 


varia- 














tion 


original order +25% 


1 


iVi V 


444 


607 


oozy .ozi 


7 AOS 


i .y^o 


2 


MV 


472 


632 


zouu.ooo 


o.o± 1 


497 


Q 
O 


MV 


238 


306 




A Q77 


2 264 


4 


MV 


263 


360 


1 ^91 RAO 




9 nQi 


o 


MV 


217 


286 


lUOO.OU / 


A 


2 147 




MV 


191 


258 


QQ7 f\Q^ 


4 644 


9 nni 


7 


MV 


884 


1323 




fi Rn7 




8 


MV 


366 


477 


2393.787 


6.54 


2.241 


9 


MV 


218 


290 


1251.694 


5.742 


1.708 


10 


MV 


201 


255 


1212.602 


6.033 


1.527 


11 


MV 


202 


266 


1232.495 


6.101 


1.63 


12 


.MV 


464 


623 


3134.5 


6.755 


2.236 






original order 


+50% 






1 


MV 


444 


729 


3260.014 


7.342 


1.992 


2 


MV 


472 


759 


2501.373 


5.3 


2.294 


3 


MV 


238 


367 


1165.202 


4.896 


1.785 


4 


MV 


263 


432 


1314.276 


4.997 


1.556 


5 


MV 


217 


343 


1044.607 


4.814 


1.555 


6 


MV 


191 


310 


874.358 


4.578 


1.604 


7 


MV 


884 


1588 


5781.723 


6.54 


2.091 


8 


MV 


366 


573 


2362.72 


6.456 


1.936 


9 


MV 


218 


348 


1228.775 


5.637 


1.43 


10 


MV 


201 


306 


1189.092 


5.916 


1.239 


11 


MV 


202 


319 


1218.99 


6.035 


1.437 


12 


MV 


464 


748 


3102.316 


6.686 


1.853 


original order +75% 


1 


MV 


444 


850 


3228.822 


7.272 


1.905 


2 


MV 


472 


885 


2489.82 


5.275 


1.978 


3 


MV 


238 


428 


1151.462 


4.838 


1.623 


4 


MV 


263 


501 


1301.374 


4.948 


1.316 


5 


AlV 


217 


100 


1036.017 


1.771 


1.313 


6 


MV 


191 


362 


861.737 


4.512 


1.462 


7 


MV 


884 


1853 


5734.982 


6.488 


1.994 


8 


MV 


366 


668 


2350.599 


6.422 


1.821 


9 


MV 


218 


406 


1209.014 


5.546 


1.311 


10 


MV 


201 


357 


1164.48 


5.793 


1.13 


11 


MV 


202 


372 


1201.455 


5.948 


1.349 


12 


MV 


464 


873 


3088.206 


6.656 


1.701 


original order +100% 


1 


MV 


444 


972 


3213.273 


7.237 


1.864 


2 


MV 


472 


1012 


2481.614 


5.258 


1.746 


3 


MV 


238 


490 


1137.444 


4.779 


1.516 


4 


MV 


263 


576 


1283.073 


4.879 


1.207 


5 


MV 


217 


458 


1025.907 


4.728 


1.215 


6 


MV 


191 


414 


844.295 


4.42 


1.406 


7 


MV 


884 


2118 


5673.277 


6.418 


1.941 


8 


MV 


366 


764 


2240.144 


6.121 


1.985 


9 


AlV 


218 


16 1 


1186.383 


5.142 


1.278 


10 


Ai\" 


201 


108 


1110.022 


5.675 


1.072 


11 


MV 


202 


426 


1188.545 


5.884 


1.305 


12 


MV 


464 


998 


3080.978 


6.64 


1.619 



Table 8: Betweenness for assortative low node degree strategy Medium Voltage samples evolution. 
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when the edges of the networks are doubled compared to the original size; for some samples the improvement is 
more than three orders of magnitude. The improvement in the clustering coefficient tend to have a logistic trend 
in the evolution steps considered (Figure |6|. In general, robustness doubles when the size of the graph doubles, 
in particular, one notices a sharp increase in this metric between the first and the second evolution step. In 
addition, it is interesting to note that for some samples the highest value of robustness is reached not when the 
connectivity is the highest, but in intermediate steps of evolution. For example, sample ^1 has better robustness 
in evolution step two than step three with values of 0.368 and 0.363; for sample #2, the highest value is reached 
in step three of the evolution process. An aspect that needs to be highlighted is the effect of this evolution on 
robustness against random and targeted attacks. The addition of edges to contrast the first type of attacks is 



always beneficial, increasing sub- linearly while more edges are added (Figure 7a); on the other hand, to contrast 



the second type of attack the maximum robustness is obtained in the second step of edges addition (Figure 7b ) 



For the redundant path robustness, the same considerations done for the characteristic path length apply. There 
is a reduction of 28%, compared to the initial value of the samples already when the networks are evolved with 
a 25% increase in the size. The maximal reduction in the redundant average path length is obtained when the 
edges are doubled with the average 10*'' path that is about 55% less than the initial value. 



0.60000 



0.50000 



0.40000 



0.30000 



0.20000 



0.10000 



0.00000 

Original 




• #1 

*(« 

■M-S7 
■•■(IIO 

■#■((11 



Evol +25% 



Evol +50% 



Evoi +75% 



Evol +100% 



Figure 6: Evolution of the clustering coefficient metric. 



Table 10 contains for each sample of the Medium Voltage Grid (column one) the values of metrics related to 
betweenness. In addition to order and size (columns three and four), average betweenness is provided in columns 
five, while a value of average betweenness normalized by the order of the graph is computed in column six in 




*»2 0.26000 

V #3 



Evol +76% EvoH-100% 




Original 



(a) Random node removal robustness. (b) Targeted node removal robustness. 

Figure 7: Evolution of robustness sub-metrics. 
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order to compare the different samples. A measure of the statistical variation of betweenness is the coefficient of 
variation which is shown in the seventh column. 

Considering betweenness, the triangle closure strategy involving the nodes with highest node degree is ben- 
eficial in reducing the average betweenness of all the Medium Voltage samples. The reduction in the average 
betweenness is around 45% compared to the original value. In the first step of the network evolution, the between- 
ness on average reduces about 20%. Some samples (i.e., samples #5 and#6) reach even higher reduction already 
in the first phase of the evolution (45% and 66% respectively). The same trend is followed by the normalized 
value of average betweenness divided by the order of the network; already in the first step no sample exceeds 
seven for the betweenness to order ratio. Considering the variability of betweenness in each of the samples we 
note a general increase for this metric. All the twelve samples present an increase in the coefficient of variation 
that increments at each stage of the evolution. Such behavior is due to a slower decrease of the standard deviation 
of betweenness compared to its average. This strategy of adding the connection, by providing more connections 
between neighbors of a node, has less effect on the nodes at the edge of the network which have a marginal role 
in providing shortest path management. A proof is that the median of the betweenness in each of the samples at 
every state of the evolution process is zero, that is the majority of nodes are terminal nodes that are not involved 
in any path between other nodes. 



Dissortative node degree evolution 



Table 11 contains for each sample of the Medium Voltage Grid (column one) the values for main topological 
quantities: order and size in columns three and four, respectively; average node degree in the fifth column; the 
characteristic path length is reported in column six; the clustering coefficient follows in column seven; robustness 
is shown in eight column; the cost in term of redundant path length closes the data series (column nine) . 

The primary focus of the dissortative node degree strategy is to connect nodes with small node degree with 
nodes with high nodes degree. A general consideration that can be made analyzing the data in Table [TT] is 
that the most of the improvement for all the metrics is obtained after the third step of evolution, the last 
increment in connectivity is scarcely beneficial. Considering the characteristic path length metric, the increase 
in the connectivity with such an evolution strategy is beneficial. The samples score for this metric all below 6.3 
showing an improvement of 55% in general when the number of edges is doubled. As already remarked, the final 
evolution step is slightly beneficial, providing a reduction in characteristic path length of just 1% compared to 
step three. The clustering coefficient evolution for such strategy provides benefits that are not more than two 
orders of magnitude. It is also interesting to remark that according to this evolution strategy for many samples 
the peak value of the clustering coefficient is not obtained when the maximum amount of edges are added, but 
in the first or second step of the evolution. In fact, more connectivity between heterogeneous nodes (in terms of 
node degree) leads to a situation where nodes with a big neighborhood of nodes in which the new node is not 
likely have other connections, therefore reducing the clustering coefficient of the whole network. The graphical 
representation of the evolution of the clustering coefficient is shown in Figure |10[ In general, robustness triples 
when the size of the graph doubles, in particular a notable increase takes place until the third evolution step. 
As mentioned for other strategies, the addition of edges provides benefits in dealing with random attacks, but 
for some samples (i.e., samples #2, #5, #6) the improvement is particularly high also for the targeted attacks 



having the two values that compose the robustness metric (cf. Appendix A) that score almost equally around 
0.5-0.6. Newman [SS' explains that assortative networks that tight high node degree have a sort of redundancy 
in the main cluster that connects the nodes with high degree, while this is absent in a dissortative network. In 
the evolution following this strategy we see that the new connectivity is particularly beneficial in contrasting the 
targeted attacks. The new established connections tend to bond nodes with low degree to the already highly 
connected nodes. On the one hand, this strategy reinforces the connectivity of already established hubs; on the 
other hand, it creates small hubs at the periphery of the network since the nodes with smallest connectivity tend 
to be more and more connected with the central nodes while more new edges are attached. Therefore, these 
new external and redundant hubs are new targets for the high node removal policy once they have sufficient 
connectivity. This improvement in the reliability is empirically evident by the results of the robustness metric 



(Figures 9a and 9b ) that show the evolution in the random and targeted node removal situations, respectively. 
One sees that the increase in connectivity is beneficial to contrast random attacks. A plateau for random failures 
around 0.6-0.7 is reached after the second evolution step. More benefits are obtained against the attacks that 
target high degree nodes: about half of the samples experience an improvement in the metric about one order of 
magnitude. 

Considering the values for the redundant path robustness, the same considerations for the characteristic path 
length apply. There is already a reduction of 50%, compared to the initial value of the samples, already when 
the networks are evolved with a 25% increase in the size of the original graph. The additional connectivity that 
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Sample ID 


Network 
type 


Order 


Size 


Avg. 
deg. 


CPL 


CC 


Removal 
robust- 
ness 

(.RobN) 


Redundancy 
cost 


original order +25% 


1 


MV 


444 


607 


2.734 


11.178 


0.08388 


0.231 


14.166 


2 


MV 


472 


632 


2.678 


15.929 


0.15681 


0.191 


20.221 


3 


MV 


238 


306 


2.571 


11.768 


0.18787 


0.199 


15.655 


4 


MV 


263 


360 


2.738 


11.271 


0.16888 


0.274 


14.735 


5 


MV 


217 


286 


2.636 


9.093 


0.15887 


0.269 


14.314 


6 


MV 


191 


258 


2.702 


7.763 


0.17561 


0.268 


10.813 


7 


MV 


884 


1323 


2.993 


8.552 


0.19566 


0.355 


11.051 


8 


MV 


366 


477 


2.607 


13.453 


0.18984 


0.177 


16.282 


9 


MV 


218 


290 


2.661 


9.9,33 


0.11586 


0.324 


19.088 


10 


MV 


201 


255 


2.537 


13.400 


0.19964 


0.185 


15.747 


11 


MV 


202 


266 


2.634 


11., 376 


0.21345 


0.229 


15.254 


12 


MV 


464 


623 


2.685 


11.848 


0.14181 


0.245 


15.371 


original order +50% 


1 


MV 


444 


729 


3.284 


10.025 


0.19503 


0.368 


11.951 


2 


MV 


472 


759 


3.216 


13.900 


0.30131 


0.267 


17.160 


3 


MV 


238 


367 


3.084 


9.848 


0.26126 


0.,321 


12.735 


4 


MV 


263 


432 


3.285 


9.538 


0.26645 


0.,330 


11.308 


5 


MV 


217 


343 


3.161 


8.412 


0.28220 


0.406 


11.819 


6 


MV 


191 


310 


3.246 


7.116 


0.22559 


0.334 


8.789 


7 


MV 


884 


1588 


3.593 


7.871 


0.33041 


0.394 


9.865 


8 


MV 


366 


573 


3.131 


11.775 


0.30838 


0.248 


13.960 


9 


MV 


218 


348 


3.193 


8.323 


0.27785 


0.296 


9.878 


10 


MV 


201 


306 


3.045 


11.570 


0.32079 


0.229 


12.207 


11 


MV 


202 


319 


3.158 


10.634 


0.31854 


0.252 


13.553 


12 


MV 


464 


748 


3.224 


10.748 


0.27124 


0.328 


12.741 


original order +75% 


1 


MV 


444 


850 


3.829 


8.796 


0.32400 


0.363 


10.642 


2 


MV 


472 


885 


3.750 


12.008 


0.41309 


0.,350 


13.903 


3 


MV 


238 


428 


3.597 


8.880 


0.40270 


0.288 


10.186 


4 


MV 


263 


504 


3.833 


8.950 


0.37425 


0.370 


10.482 


5 


AiV 


217 


100 


3.()iS7 


7.191 


0.12(i01 


0.301 


8.353 


6 


MV 


191 


362 


3.791 


6.126 


0.34172 


0.313 


7.584 


7 


MV 


884 


1853 


4.192 


7.323 


0.39933 


0.385 


8.965 


8 


MV 


366 


668 


3.650 


10. ,363 


0.38451 


0.257 


11.336 


9 


MV 


218 


406 


3.725 


7.297 


0.36753 


0.,397 


8.770 


10 


MV 


201 


357 


3.552 


10.185 


0.42392 


0.279 


11.888 


11 


MV 


202 


372 


3.683 


9.963 


0.35298 


0.284 


11.178 


12 


MV 


464 


873 


3.763 


9.477 


0.34785 


0.362 


11.815 


original order +100% 


1 


MV 


444 


972 


4.378 


7.869 


0.39146 


0.385 


9.194 


2 


MV 


472 


1012 


4.288 


10.839 


0.53045 


0.342 


11.835 


3 


MV 


238 


490 


4.118 


8.086 


0.56850 


0.358 


9.286 


4 


MV 


263 


576 


4.380 


8.057 


0.44110 


0.414 


9.389 


5 


MV 


217 


458 


4.221 


6.597 


0.50752 


0.389 


7.457 


6 


MV 


191 


414 


4.335 


5.700 


0.46025 


0.361 


7.443 


7 


MV 


884 


2118 


4.792 


6.745 


0.43981 


0.417 


8.210 


8 


MV 


366 


764 


4.175 


9.393 


0.51501 


0.313 


10.760 


9 


MV 


218 


464 


4.257 


6.569 


0.45466 


0.442 


7.729 


10 


MV 


201 


408 


4.060 


9.190 


0.55617 


0.325 


9.556 


11 


MV 


202 


426 


4.218 


8.338 


0.47119 


0.332 


11.194 


12 


MV 


464 


998 


4.302 


8.486 


0.41405 


0.405 


9.871 



Table 9: Metrics for triangle closure strategy Medium Voltage samples evolution. 
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Sample ID 


Network 


Order 


Size 


Avg. be- 


Avg. 


Coeff. 




type 






tweenness 


betw/order 


vEiria- 














tion 


original order +25% 


1 


iVi V 


444 


607 


000 ^ . i 10 


0.0^0 


9 0^9 


2 


MV 


472 


632 


Aa7-\ qco 

^±\) ( -L .OOZ 


y . y ( 


1 849 


Q 
O 


MV 


238 


306 




R Q1 


J- .00 1 


4 


MV 


263 


360 




A 1 nf> 


1 614 


o 


MV 


217 


286 




8 232 


1 724 




MV 


191 


258 


1 9Q7 74R 


fi 7Q4 
u. 1 y^ 


j-.yoo 


7 


MV 


884 


1323 






0. ±0^ 


8 


MV 


366 


477 


4258.659 


11.636 


1.916 


9 


MV 


218 


290 


1100.51 


5.048 


1.67 


10 


MV 


201 


255 


2982.928 


14.84 


1.303 


11 


MV 


202 


266 


2334.081 


11.555 


1.56 


12 


.MV 


464 


623 


3076.954 


6.631 


1.826 






original order 


+50% 






1 


MV 


444 


729 


3183.386 


7.17 


2.213 


2 


MV 


472 


759 


4035.469 


8.55 


1.969 


3 


MV 


238 


367 


1387.793 


5.831 


1.791 


4 


MV 


263 


432 


950.802 


3.615 


1.731 


5 


MV 


217 


343 


1625.935 


7.493 


1.817 


6 


MV 


191 


310 


1180.635 


6.181 


1.73 


7 


MV 


884 


1588 


6298.27 


7.125 


3.405 


8 


MV 


366 


573 


3625.364 


9.905 


1.777 


9 


MV 


218 


348 


898.163 


4.12 


1.848 


10 


MV 


201 


306 


2562.557 


12.749 


1.319 


11 


MV 


202 


319 


2194.406 


10.863 


1.582 


12 


MV 


464 


748 


2901.182 


6.253 


1.906 


original order +75% 


1 


MV 


444 


850 


2782.5 


6.267 


2.158 


2 


MV 


472 


885 


3468.145 


7.348 


2.085 


3 


MV 


238 


428 


1179.631 


4.956 


1.828 


4 


MV 


263 


501 


873.116 


3.32 


1.711 


5 


AIV 


217 


100 


1389. .383 


0.103 


1.939 


6 


MV 


191 


362 


1030.677 


5.396 


1.827 


7 


MV 


884 


1853 


5736.492 


6.489 


3.495 


8 


MV 


366 


668 


3097.301 


8.463 


1.936 


9 


MV 


218 


406 


809.837 


3.715 


1.818 


10 


MV 


201 


357 


2130.351 


10.599 


1.54 


11 


MV 


202 


372 


2060 


10.198 


1.505 


12 


MV 


464 


873 


2733.855 


5.892 


1.974 


original order +100% 


1 


MV 


444 


972 


2540.609 


5.722 


2.093 


2 


MV 


472 


1012 


3034.596 


6.429 


2.323 


3 


MV 


238 


490 


1065.207 


4.476 


2.108 


4 


MV 


263 


576 


807.934 


3.072 


1.719 


5 


MV 


217 


458 


1223.467 


5.638 


2.064 


6 


MV 


191 


414 


941.228 


4.928 


1.969 


7 


MV 


884 


2118 


5292.311 


5.987 


3.439 


8 


MV 


366 


764 


2757.96 


7.535 


2.146 


9 


AlV 


218 


16 1 


681.923 


3.142 


2.051 


10 


Ai\" 


201 


108 


1837. 15 1 


9.112 


1.703 


11 


MV 


202 


426 


1710.355 


8.467 


1.666 


12 


MV 


464 


998 


2545.442 


5.486 


2.042 



Table 10: Betweenness for triangle closure strategy Medium Voltage samples evolution. 
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(e) 4*'' stage of evolution (i.e., +100% edges). 
Figure 8: Stages of evolution of Medium Voltage sample #8 following the triangle closure strategy. 
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Original Evoi +25% Evoi +50% Evoi+75% Evol+100% Originai Evol+25% Evol+50% Evol+75% Evol+100% 



(a) Random node removal robustness. (b) Targeted node removal robustness. 

Figure 9: Evolution of robustness sub- metrics. 



0.10000 



0.09000 




Figure 10: Evolution of the clustering coefficient metric. 
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is provided in the last evolution step is marginally beneficial, as seen for the characteristic path length, providing 
just an additional 2% reduction in the path length. 

Table 12 contains for each sample of the Medium Voltage Grid (column one) the values of metrics related to 
betweenness. In addition to order and size (columns three and four), average betweenness is provided in columns 
five, while a value of average betweenness normalized by the order of the graph is computed in column six in 
order to compare the different samples. A measure of the statistical variation of betweenness is the coefficient of 
variation which is shown in the seventh column. 

Considering betweenness, the dissortative strategy involving the nodes with lowest node degree is substantially 
beneficial in reducing the average betweenness of all the Medium Voltage samples just in the first step of the 
evolution. In fact, with the addition of just a quarter of the initial number of links, the average betweenness 
reduces around 45% of the original value. The additional three evolution stages contribute modestly in further 
reducing this metric (just 5%). The same trend is followed by the normalized value of average betweenness 
divided by the order of the network; already in the first step just one sample exceeds 6 for the betweenness to 
order ratio. Considering the variability (i.e., coefficient of variation) of betweenness in each of the samples we 
note a general increase for this metric, only after the last stage of the evolution the coefficient of variation returns 
to values closer to the initial value. Such behavior is due to the decrease in the standard deviation of betweenness 
that is slower compared to the average betweenness. However, this strategy of adding connections between the 
nodes that are at the opposite in their degree, tends to make them more evenly involved in the shortest paths. 
A proof is that the median of the betweenness for seven of the twelve samples is higher than zero. 



Least distance evolution 



Table 13 contains for each sample of the Medium Voltage Grid (column one) the values for main topological 
quantities: order and size in columns three and four, respectively; average node degree in the fifth column; the 
characteristic path length is reported in column six; the clustering coefficient follows in column seven; robustness 
is shown in eight column; the cost in term of redundant path length closes the data series (column nine). 

The least distance strategy aims at connecting the nodes (i.e., substations in the physical Distribution Grid) 
that are geographically closer to each other, therefore there is no specific topological strategy in such a type of 
evolution. Considering the characteristic path length metric, the increase in the connectivity is beneficial. The 
samples on average reduce the characteristic path length by 33% after the doubling of the size of the graph. 
However, we note some variability: some samples improve more than 50%, while for the biggest sample (i.e., 
sample #7) the improvement is just 20%. Even in this type of evolution the tendency is that of saturation 
with further steps in connectivity that are less beneficial. The clustering coefficient evolution for such strategy 
provides benefits that are not more than the three orders of magnitude. The benefits obtained are comparable 
by those achieved with the triangle closure strategy at least in the first two stages of evolution; when more 
links are added the nodes that are connected are less and less belonging to the same neighborhood (i.e., nodes 
with longer distances are added), and this is the main difference at the final stage of evolution compared to the 
triangle closure strategy. For the metric concerning robustness, the average over all samples gives a doubling of 
the metric compared to its original value. The greatest improvement in robustness is achieved after the first two 
steps of the evolution (about 75% improvement), a lower increase is achieved in the following two steps (about 
30%). This strategy of evolution benefits both the random and targeted attacks towards nodes. In particular. 



in the forth stage of the evolution the components of the robustness metric (cf. Appendix A I have almost the 
same value, actually for some samples it is more damaging the random attack than the attack targeting the 
most connected nodes. Such situation is represented in Figures |12a| and |12b| Considering the values for the 
redundant path robustness, the same considerations done for the characteristic path length apply. The overall 
reduction in the redundant average path length is about 50% compared to the initial value of the samples. Just 
the first evolution step provides a reduction of more than 30%. Even in this path-related metric we see a certain 
variability by samples that have consistent variations (e.g., sample #3), and samples whose improvement are 
limited (e.g., sample #7). 



Table 14 contains for each sample of the Medium Voltage Grid (column one) the values of metrics related to 
betweenness. In addition to order and size (columns three and four), average betweenness is provided in columns 
five, while a value of average betweenness normalized by the order of the graph is computed in column six in 
order to compare the different samples. A measure of the statistical variation of betweenness is the coefficient of 
variation which is shown in the seventh column. 

The strategy of adding edges between nodes that are geographically closest provides small benefits in terms 
of reduction of betweenness compared to the other strategies proposed. Actually, in the first step of the evolution 
there are four samples (samples #2, #3, #9 and #12) whose betweenness increases. In the last step of the 
evolution, over all the samples, the reduction is about 20%, despite samples #4, #9 and #12 that have an 
increment in the betweenness metric. The same trend is followed by the normalized value of average betweenness 
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Sample ID 


Network 
type 


Order 


Size 


Avg. 
deg. 


CPL 


CC 


Removal 
robust- 
ness 

(.RobN) 


Redundancy 
cost 


original order +25% 


1 


MV 


444 


607 


2.734 


6.800 


0.00687 


0.,321 


9.311 


2 


MV 


472 


632 


2.678 


5.809 


0.01625 


0.,344 


8.604 


3 


MV 


238 


306 


2.571 


5.842 


0.00491 


0.319 


9.112 


4 


MV 


263 


360 


2.738 


5.580 


0.01966 


0.364 


8.399 


5 


MV 


217 


286 


2.636 


4.903 


0.01245 


0.333 


7.992 


6 


MV 


191 


258 


2.702 


4.979 


0.04375 


0.351 


7.648 


7 


MV 


884 


1323 


2.993 


6.866 


0.00906 


0.411 


8.235 


8 


MV 


366 


477 


2.607 


6.314 


0.03204 


0.315 


9.743 


9 


MV 


218 


290 


2.661 


5.970 


0.04170 


0.340 


8.548 


10 


MV 


201 


255 


2.537 


6.290 


0.01264 


0.310 


8.746 


11 


MV 


202 


266 


2.634 


6.413 


0.00800 


0.,326 


9.951 


12 


MV 


464 


623 


2.685 


6.470 


0.03983 


0.316 


9.295 


original order +50% 


1 


MV 


444 


729 


3.284 


6.341 


0.01046 


0.370 


9.079 


2 


MV 


472 


759 


3.216 


5.645 


0.01730 


0.425 


7.778 


3 


MV 


238 


367 


3.084 


5.586 


0.00701 


0.439 


7.748 


4 


MV 


263 


432 


3.285 


5.477 


0.04544 


0.432 


7.482 


5 


MV 


217 


343 


3.161 


4.801 


0.02574 


0.418 


7.258 


6 


MV 


191 


310 


3.246 


4.647 


0.08098 


0.436 


7.002 


7 


MV 


884 


1588 


3.593 


6.577 


0.01368 


0.409 


8.199 


8 


MV 


366 


573 


3.131 


6.178 


0.05844 


0.360 


8.085 


9 


MV 


218 


348 


3.193 


5.751 


0.03649 


0.435 


8.133 


10 


MV 


201 


306 


3.045 


6.240 


0.04287 


0.362 


8.330 


11 


MV 


202 


319 


3.158 


5.995 


0.00745 


0.355 


8.982 


12 


MV 


464 


748 


3.224 


6.323 


0.02381 


0.366 


8.452 


original order +75% 


1 


MV 


444 


850 


3.829 


6.246 


0.02326 


0.376 


9.187 


2 


MV 


472 


885 


3.750 


5.590 


0.01667 


0.635 


7.656 


3 


MV 


238 


428 


3.597 


5.308 


0.01065 


0.462 


6.847 


4 


MV 


263 


504 


3.833 


5.359 


0.05212 


0.634 


7.009 


5 


AiV 


217 


100 


:-!.()iSV 


l.()71 


0.02781 


0.583 


(i.9 19 


6 


MV 


191 


362 


3.791 


4.563 


0.07813 


0.500 


6.212 


7 


MV 


884 


1853 


4.192 


6.532 


0.02022 


0.409 


8.163 


8 


MV 


366 


668 


3.650 


6.012 


0.05044 


0.409 


7.997 


9 


MV 


218 


406 


3.725 


5.297 


0.02018 


0.434 


7.492 


10 


MV 


201 


357 


3.552 


6.000 


0.03812 


0.376 


7.988 


11 


MV 


202 


372 


3.683 


5.764 


0.02148 


0.357 


7.971 


12 


MV 


464 


873 


3.763 


6.056 


0.03175 


0.363 


7.598 


original order +100% 


1 


MV 


444 


972 


4.378 


6.089 


0.02662 


0.380 


8.296 


2 


MV 


472 


1012 


4.288 


5.089 


0.01730 


0.635 


7.110 


3 


MV 


238 


490 


4.118 


5.097 


0.02329 


0.479 


6.761 


4 


MV 


263 


576 


4.380 


5.011 


0.04146 


0.654 


6.320 


5 


MV 


217 


458 


4.221 


4.551 


0.04456 


0.585 


6.222 


6 


MV 


191 


414 


4.335 


4.547 


0.08621 


0.509 


6.516 


7 


MV 


884 


2118 


4.792 


6.278 


0.01899 


0.407 


7.475 


8 


MV 


366 


764 


4.175 


5.929 


0.02961 


0.404 


8.028 


9 


MV 


218 


464 


4.257 


5.189 


0.01336 


0.433 


6.919 


10 


MV 


201 


408 


4.060 


5.840 


0.02601 


0.383 


7.391 


11 


MV 


202 


426 


4.218 


5.637 


0.03385 


0.357 


7.656 


12 


MV 


464 


998 


4.302 


5.975 


0.02619 


0.367 


8.007 



Table 11: Metrics for dissortative node degree strategy Medium Voltage samples evolution. 
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Sample ID 


Network 


Order 


Size 


Avg. be- 


Avg. 


Coeff. 




type 






tweenness 


betw/order 


vEiria- 














tion 


original order +25% 


1 


iVi V 


444 


607 


Z 1 Do. ( DZ 


a 99^ 
D.ZZO 


9 9^4 


2 


MV 


472 


632 


99QO 1 QO 




O . tj 1 o 


Q 
O 


MV 


238 


306 


1 1 n'i «fio 

1 i-UO.OUO 




2 652 


4 


MV 


263 


360 


1161 52 


4 416 


z.oyo 


o 


MV 


217 


286 




A 9R8 


9 SQ7 
z.oy 1 




MV 


191 


258 


7AQ Q1 


^ Q9fi 


9 78*^ 
Z. 1 oo 


7 


MV 


884 


1323 






9 SI fi 

Z.O J.O 


8 


MV 


366 


477 


2084.738 


5.696 


2.577 


9 


MV 


218 


290 


1067.579 


4.897 


2.061 


10 


MV 


201 


255 


1125.408 


5.599 


1.729 


11 


MV 


202 


266 


1149 


5.688 


1.739 


12 


.MV 


464 


623 


2713.298 


5.818 


2.77 






original ofder 


+50% 






1 


MV 


444 


729 


2630.233 


5.924 


1.979 


2 


MV 


472 


759 


2178.939 


4.616 


2.644 


3 


MV 


238 


367 


1011.731 


4.251 


2.063 


4 


MV 


263 


432 


1124.138 


4.274 


1.948 


5 


MV 


217 


343 


894.523 


4.122 


2.235 


6 


MV 


191 


310 


725.905 


3.801 


2.11 


7 


MV 


884 


1588 


4969.261 


5.621 


2.495 


8 


MV 


366 


573 


2002.795 


5.472 


1.96 


9 


MV 


218 


348 


1019.053 


4.675 


1.648 


10 


MV 


201 


306 


1104.204 


5.494 


1.396 


11 


MV 


202 


319 


1063.747 


5.266 


1.517 


12 


MV 


464 


748 


2619.583 


5.646 


2.276 


original order +50% 


1 


MV 


444 


850 


2588.741 


5.83 


1.842 


2 


MV 


472 


885 


2141.579 


4.537 


2.301 


3 


MV 


238 


428 


970.009 


4.076 


1.763 


4 


MV 


263 


501 


1108.721 


4.216 


1.667 


5 


AIV 


217 


100 


850. 327 


3.91(> 


1.858 


6 


MV 


191 


362 


704.274 


3.687 


1.775 


7 


MV 


884 


1853 


4906.982 


5.551 


2.34 


8 


MV 


366 


668 


1936.04 


5.29 


1.778 


9 


MV 


218 


406 


947.675 


4.347 


1.478 


10 


MV 


201 


357 


1052.286 


5.235 


1.278 


11 


MV 


202 


372 


1016.879 


5.034 


1.394 


12 


MV 


464 


873 


2449.978 


5.28 


2.082 


original order +50% 


1 


MV 


444 


972 


2510.314 


5.654 


1.668 


2 


MV 


472 


1012 


2109.969 


4.47 


2.068 


3 


MV 


238 


490 


944.942 


3.97 


1.574 


4 


MV 


263 


576 


1058.61 


4.025 


1.58 


5 


MV 


217 


458 


842.523 


3.883 


1.708 


6 


MV 


191 


414 


690.221 


3.614 


1.569 


7 


MV 


884 


2118 


4818.952 


5.451 


2.128 


8 


MV 


366 


764 


1905.977 


5.208 


1.682 


9 


AlV 


218 


16 1 


925.225 


1.241 


1.391 


10 


Ai\" 


201 


108 


1018. 80() 


5.0(i9 


1.177 


11 


MV 


202 


426 


998.545 


4.943 


1.347 


12 


MV 


464 


998 


2408.728 


5.191 


1.986 



Table 12: Betweenness for dissortative node degree strategy Medium Voltage samples evolution. 
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(c) 2"'' stage of evolution (i.e., +50% edges) (d) 3'''* stage of evolution (i.e., +75% edges). 




(e) 4*'' stage of evolution (i.e., +100% edges). 
Figure 11: Stages of evolution of Medium Voltage sample #8 following the dissortative strategy. 
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•3 #6 0.40000 




EvoH-76% Evol-HOO% 



(a) Random node removal robustness. (b) Targeted node removal robustness. 

Figure 12: Evolution of robustness sub-metrics. 



divided by the order of the network; in the last step of evolution there are three sample which have a betweenness 
to order ratio which is higher than 9. Considering the variability of betweenness, we see a general increase in 
the coefficient of variation. Only sample #5 and sample #8 show a decrease for such a metric. This strategy 
of adding the connections does not clearly have a strategy from the topology point of view. Therefore, it makes 
sense that there is not a clear trend in the results of betweenness since every sample has different characteristics 
in the physical distance between the nodes (substations in the network). 



Random evolution 



Table 15 contains for each sample of the Medium Voltage Grid (column one) the values for main topological 
quantities: order and size in columns three and four, respectively; average node degree in the fifth column; the 
characteristic path length is reported in column six; the clustering coefRcient follows in column seven; robustness 
is shown in eight column; the cost in term of redundant path length closes the data series (column nine). 

The random evolution strategy is mainly intended to consider the difference between evolutions driven by a 
specific target or property of nodes and an evolution driven by mere random picking of nodes. Considering the 
characteristic path length metric the increase in the connectivity with such a strategy is extremely beneficial. 
The samples score for this metric below 5 after the third evolution step, with the exception of sample #12 which 
scores 5.014. The last step of connectivity addition brings the total benefit for reduction in characteristic path 
length to a value higher than 65%. The clustering coefficient evolution for such strategy provides limited benefits. 
In fact, to improve the clustering a special addition scheme needs to be followed targeting the neighbors of a 
node. The situation varies between the samples, but the majority of them improve in the clustering coefficient 
by one order of magnitude when the edges are double the initial number. Figure [14] shows that, although a 
general improvement of the metric, there is not a common trend or tendency for all the samples: some reach 
a peak in intermediate stages of evolution (e.g., sample #4), others a rapid increase (e.g., sample #12), others 
have a slightly decreasing (e.g., sample #2) or monotonously increasing (e.g., sample #7) clustering coefficient. 
Robustness is the metric that improves the most. Values for this metric of 0.7 are obtained already after the 
third evolution step. Edge addition following this strategy provides an increase in both the random and targeted 
component of this metric. In particular, the random addition of nodes tends to make the whole network more 
robust such that the effects of the random and targeted attack towards nodes have the same effect on the network. 
All the samples reach a value close to 0.75 that seems an upper bound in such conditions of order and size; we 
already noted such an upper bound behavior for various types of synthetic topologies in our previous work [55] . 
Considering the values for the redundant path robustness the same considerations done for the characteristic 
path length apply. Especially, the rates of reduction (in percentage) of the 10*'' average path length at each 
evolution step are exactly the same as those shown by the characteristic path length . When the number of edges 
is doubled the redundant path length for all the samples stays below 6.5 which is just 2 hops higher than the 
characteristic path length of the network. 



Table 16 contains for each sample of the Medium Voltage Grid (column one) the values of metrics related to 
betweenness. In addition to order and size (columns three and four), average betweenness is provided in columns 
five, while a value of average betweenness normalized by the order of the graph is computed in column six in 
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Sample ID 


Network 
type 


Order 


Size 


Avg. 
deg. 


CPL 


CC 


Removal 
robust- 
ness 

(.RobN) 


Redundancy 
cost 


original order +25% 


1 


MV 


444 


607 


2.734 


9.536 


0.13338 


0.293 


12.753 


2 


MV 


472 


632 


2.678 


14.726 


0.14207 


0.191 


17.222 


3 


MV 


238 


306 


2.571 


10.395 


0.14795 


0.236 


13.316 


4 


MV 


263 


360 


2.738 


10.431 


0.15275 


0.252 


13.923 


5 


MV 


217 


286 


2.636 


8.815 


0.08661 


0.308 


11.384 


6 


MV 


191 


258 


2.702 


8.132 


0.11575 


0.289 


10.958 


7 


MV 


884 


1323 


2.993 


8.527 


0.12721 


0.398 


10.388 


8 


MV 


366 


477 


2.607 


12.838 


0.14166 


0.184 


16.860 


9 


MV 


218 


290 


2.661 


9.606 


0.13026 


0.242 


12.764 


10 


MV 


201 


255 


2.537 


12.570 


0.13361 


0.204 


15.755 


11 


MV 


202 


266 


2.634 


9.284 


0.12775 


0.286 


13.146 


12 


MV 


464 


623 


2.685 


10.884 


0.11793 


0.278 


13.828 


original order +50% 


1 


MV 


444 


729 


3.284 


8.647 


0.21344 


0.404 


10.725 


2 


MV 


472 


759 


3.216 


13.190 


0.19715 


0.289 


15.654 


3 


MV 


238 


367 


3.084 


9.876 


0.20415 


0.276 


11.840 


4 


MV 


263 


432 


3.285 


9.565 


0.22611 


0.,388 


11.919 


5 


MV 


217 


343 


3.161 


8.352 


0.16083 


0.342 


10.662 


6 


MV 


191 


310 


3.246 


7.768 


0.15339 


0.329 


10.455 


7 


MV 


884 


1588 


3.593 


8.097 


0.18306 


0.463 


9.674 


8 


MV 


366 


573 


3.131 


12.177 


0.19587 


0.224 


15.782 


9 


MV 


218 


348 


3.193 


9.184 


0.20451 


0.301 


12.243 


10 


MV 


201 


306 


3.045 


8.985 


0.17791 


0.414 


11.746 


11 


MV 


202 


319 


3.158 


8.296 


0.18381 


0.397 


10.626 


12 


MV 


464 


748 


3.224 


10.014 


0.20459 


0.386 


12.681 


original order +75% 


1 


MV 


444 


850 


3.829 


8.336 


0.27250 


0.468 


10.525 


2 


MV 


472 


885 


3.750 


12.644 


0.25483 


0.333 


13.812 


3 


MV 


238 


428 


3.597 


9.120 


0.23636 


0.,346 


11.277 


4 


MV 


263 


504 


3.833 


8,912 


0.28127 


0. 112 


10.163 


5 


AiV 


217 


100 


:-!.()iSV 


7. DOS 


0.23789 


0. ill) 


10.177 


6 


MV 


191 


362 


3.791 


6.979 


0.22299 


0.451 


9.251 


7 


MV 


884 


1853 


4.192 


7.849 


0.24197 


0.581 


9.411 


8 


MV 


366 


668 


3.650 


11.858 


0.23799 


0.,357 


14.784 


9 


MV 


218 


406 


3.725 


8.760 


0.25211 


0.497 


11.866 


10 


MV 


201 


357 


3.552 


8.145 


0.20495 


0.488 


10.736 


11 


MV 


202 


372 


3.683 


7.326 


0.21618 


0.517 


9.444 


12 


MV 


464 


873 


3.763 


9.635 


0.25206 


0.400 


11.813 


original order +100% 


1 


MV 


444 


972 


4.378 


8.023 


0.31095 


0.497 


9.485 


2 


MV 


472 


1012 


4.288 


11.792 


0.31725 


0.461 


13.000 


3 


MV 


238 


490 


4.118 


8.842 


0.26450 


0.498 


11.043 


4 


MV 


263 


576 


4.380 


8.588 


0.33434 


0.447 


9.674 


5 


MV 


217 


458 


4.221 


7.505 


0.29041 


0.458 


9.504 


6 


MV 


191 


414 


4.335 


6.574 


0.28652 


0.470 


8.522 


7 


MV 


884 


2118 


4.792 


7.565 


0.29363 


0.713 


9.048 


8 


MV 


366 


764 


4.175 


11.340 


0.28634 


0.441 


14.256 


9 


MV 


218 


464 


4.257 


8.233 


0.29032 


0.564 


10.460 


10 


MV 


201 


408 


4.060 


7.605 


0.23156 


0.575 


9.656 


11 


MV 


202 


426 


4.218 


6.963 


0.26151 


0.676 


9.056 


12 


MV 


464 


998 


4.302 


9.472 


0.28112 


0.493 


11.417 



Table 13: Metrics for least distance strategy Medium Voltage samples evolution. 
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Sample ID 


Network 


Order 


Size 


Avg. be- 


Avg. 


Coeff. 




type 






tweenness 


betw/order 


vEiria- 














tion 


original order +25% 


1 


iVi V 


444 


607 


OOO / .DZ 


O.UIO 


9 1 AS 


2 


MV 


472 


632 




1 '\ 9Qzl 


1 70zl 


Q 
O 


MV 


238 


306 




R QS1 


1 

J- .oo ( 


4 


MV 


263 


360 


99Q7 1 ^7 


S 7^4 


1 419 


o 


MV 


217 


286 


1791 047 


7 Q'^1 


1 




MV 


191 


258 


1 41 R 741 


7 41 7 


1 09*^ 
j-.yzo 


7 


MV 


884 


1323 




7 704 


z.y J.O 


8 


MV 


366 


477 


4475.009 


12.227 


1.566 


9 


MV 


218 


290 


1840.431 


8.442 


1.547 


10 


MV 


201 


255 


2398.327 


11.932 


1.275 


11 


MV 


202 


266 


1741.949 


8.624 


1.368 


12 


.MV 


464 


623 


4814.58 


10.376 


2.003 






original ofder 


+ 50% 






1 


MV 


444 


729 


3221.546 


7.256 


2.191 


2 


MV 


472 


759 


5636.531 


11.942 


1.758 


3 


MV 


238 


367 


1999.372 


8.401 


1.577 


4 


MV 


263 


432 


2120.025 


8.061 


1.495 


5 


MV 


217 


343 


1627.028 


7.498 


1.638 


6 


MV 


191 


310 


1344.032 


7.037 


1.993 


7 


MV 


884 


1588 


6475.508 


7.325 


3.024 


8 


MV 


366 


573 


4287.187 


11.714 


1.552 


9 


MV 


218 


348 


1768.297 


8.111 


1.571 


10 


MV 


201 


306 


1655.439 


8.236 


1.305 


11 


MV 


202 


319 


1539.98 


7.624 


1.374 


12 


MV 


464 


748 


4523.666 


9.749 


2.081 


original order +50% 


1 


MV 


444 


850 


3068.074 


6.91 


2.263 


2 


MV 


472 


885 


5388.36 


11.416 


1.783 


3 


MV 


238 


428 


1905.408 


8.006 


1.646 


4 


MV 


263 


501 


1956.884 


7.441 


1.41 


5 


AIV 


217 


100 


151i..3()i 


7.117 


1.68 


6 


MV 


191 


362 


1239.704 


6.491 


2.067 


7 


MV 


884 


1853 


6210.444 


7.025 


3.074 


8 


MV 


366 


668 


4115.775 


11.245 


1.58 


9 


MV 


218 


406 


1687.799 


7.742 


1.625 


10 


MV 


201 


357 


1535.122 


7.637 


1.344 


11 


MV 


202 


372 


1358.508 


6.725 


1.423 


12 


MV 


464 


873 


4329.609 


9.331 


2.14 


original order +50% 


1 


MV 


444 


972 


2937.585 


6.616 


2.309 


2 


MV 


472 


1012 


5016.908 


10.629 


1.896 


3 


MV 


238 


490 


1832.179 


7.698 


1.651 


4 


MV 


263 


576 


1852.488 


7.044 


1.484 


5 


MV 


217 


458 


1467.346 


6.762 


1.744 


6 


MV 


191 


414 


1168.878 


6.12 


2.098 


7 


MV 


884 


2118 


5956.625 


6.738 


3.085 


8 


MV 


366 


764 


3973.55 


10.857 


1.583 


9 


AlV 


218 


16 1 


1602.22 


7.35 


1.689 


10 


Ai\" 


201 


108 


1119.810 


7.21.3 


1.29 


11 


MV 


202 


426 


1296.051 


6.416 


1.408 


12 


MV 


464 


998 


4229.495 


9.115 


2.099 



Table 14: Betweenness for least distance strategy Medium Voltage samples evolution. 
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(a) Physical sample. (b) 1"* stage of evolution (i.e., +25% edges). 




(c) 2"*^ stage of evolution (i.e., +50% edges) (d) 3'"'^ stage of evolution (i.e., +75% edges), (e) 4*'* stage of evolution (i.e., +100% edges). 

Figure 13: Stages of evolution of Medium Voltage sample #8 following least distance strategy. 
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Figure 14: Evolution of the clustering coefficient metric. 



order to compare the different samples. A measure of the statistical variation of betweenness is the coefficient of 
variation which is shown in the seventh column. 

Random addition of edges in the networks provides the highest improvement in the betweenness metric. 
Average betweenness reduces to about 67% in the final step of evolution. Already from the first step of the 
network evolution, the betweenness on average reduces to slightly more than 40% of the original value. The only 
exception is sample ^4 that shows an increment in the first evolution stage. The same trend is followed by the 
normalized value of average betweenness divided by the order of the network; already in the first step only two 
samples exceeds six, and in the final stage no sample exceeds 3.6. Considering the variability of betweenness, in 
each of the samples we note a decrease for this metric. In particular, in the last stage of the evolution, all samples 
for the coefficient of variation are below the unity with the exception of sample #7 that scores just below 1.1. 
The random addition of edges tend to even split the betweenness among all the nodes of the network and this 
is shown by the sharp fall of standard deviation of betweenness which is responsible for a coefficient of variation 
that is below the unit in the last stage of the evolution. However, the highest fraction of nodes for all the samples 
has zero betweenness, the median in each step for each sample is zero. 



5.2 Evolution of the Low Voltage Distribution Grids 

The Low Voltage Grid present some peculiarities caused by their practical use of being the terminals of the 
Grid. The majority of the networks (nine out of eleven samples) have null clustering coefficient, therefore it is 
difficult to analyze the improvement compared to the original as using percentage values as done for the Medium 
Voltage networks. Another peculiarity of this network is the limited number of redundant paths (other than 
the shortest path) to connect any two nodes. In fact, paths in such networks especially the small samples are 
limited and mainly fixed by the radial topology of the network. Therefore, we consider for the evolution of such 
a metric only those samples for which there is a significant increase for the 10*'* average path length compared 
to the characteristic path length in the physical (not yet evolved) samples. In addition, we remark that for these 
samples there are no information available to us concerning the geographical position of the nodes and usually 
the samples cover a geographically very limited area with a small size and order of the corresponding graph. 
Therefore, we could not perform the distance based evolution of the Low Voltage network samples. 



Assortative high node degree evolution 



Table 17 contains for each sample of the Low Voltage Grid (column one) the values for main topological quantities: 
order and size in columns three and four, respectively; average node degree in the fifth column; the characteristic 
path length is reported in column six; the clustering coefficient follows in column seven; robustness is shown in 
eight column; the cost in term of redundant path length closes the data series (column nine). 
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The increase in the connectivity with the assortative high node degree strategy is extremely beneficial for 
the path length: just by adding 25% of links it reduces, compared to the initial value of the original graphs, on 
average more than 30%. The improvement is then reduced in the other steps, reaching 50% when the number of 
link is doubled compared to the original size. The two biggest samples benefit substantially from just the first 
evolution step: sample #5 reduces the path from almost 19 to slightly more than 6, while sample #9 from more 
than 11 to 4.75. The improvement for the clustering coefficient is considerable: all the samples after the first step 
achieve a non null coefficient and the two samples that already have a significant clustering coefficient have a ten 
fold improvement. The improvement in robustness on average is about 70% at the final step of the evolution. 
The addition of edges according to this strategy is mainly beneficial in the random attacks towards nodes, while 
the fraction of the metric that consider targeted attacks against the most connected nodes is marginally affected. 
This evolution strategy, indeed, leave the hierarchy of the nodes almost untouched reinforcing the node degree of 
those nodesthat had already an high degree at the beginning of the evolution process. For some small samples, 
the improvement in the targeted attack robustness is not monotonic, but it increases, then decreases in one 
intermediate step and then it increases again. Considering the values for the redundant path robustness the 
same considerations done for the characteristic path length apply. There is a reduction more than 35% compared 
to the initial value of the samples already when the networks are evolved with a 25% increase in the size of the 
graph. The final improvement is about 63%. However, an exception can be noted: sample #2 shows an increase 
in the value for the redundant path in the first step of the evolution compared to the initial value. Such increase 
is due to the problems that affect this metric when small networks (in terms of order and size) are considered 
since they do not have redundant path between nodes; such redundant path emerge when more connectivity is 
added to the network. 

Table 18 contains for each sample of the Low Voltage Grid (column one) the values of the metrics related to 
betweenness. In addition to order and size (columns three and four), average betweenness is provided in column 
five, while a value of average betweenness normalized by the order of the graph is computed in column six in 
order to compare the different samples. A measure of the statistical variation of betweenness is the coefficient of 
variation which is shown in the seventh column. 

Considering betweenness, the assortative strategy involving the nodes with highest node degree is beneficial 
in reducing the average betweenness of all the Low Voltage samples. The benefit that is provided after the four 
evolution steps is more than 50%. Considering the two biggest samples (i.e., sample #5 and #9) the benefits 
reach the value of 70% at the end of the evolution process, and a value closer to 65% already after the first 
step. The same trend is followed by the normalized value of average betweenness divided by the order of the 
network. Considering the variability of betweenness in each of the samples, we note a general increase for this 
metric. In particular, only one of the eleven samples presents a decrease when the connectivity in the network 
is double the original size of the network. Such behavior is due to a slower decrease of the standard deviation 
of betweenness compared to its average. This strategy of adding the connections by providing more connections 
between the nodes that already have the highest connectivity and usually highest betweenness, reduces the 
number of "bottleneck" nodes. However, this strategy is not helpful for substantially reducing the variability of 
betweenness. In fact, the median of the betweenness for the biggest samples is always zero. This aspect reinforces 
the knowledge that the majority of nodes are terminal nodes that are not involved in any path between other 
nodes. Only small samples have non zero values in the later stages of the evolution process. 



Assortative low node degree evolution 



Table 19 contains for each sample of the Low Voltage Grid (column one) the values for the main topological 
quantities: order and size in columns three and four, respectively; average node degree in the fifth column; the 
characteristic path length is reported in column six; the clustering coefficient follows in column seven; robustness 
is shown in eight column; the cost in term of redundant path length closes the data series (column nine) . 

The increase in the connectivity with the assortative low node degree strategy is beneficial for the path length. 
The first two evolution steps in particular, already provide for a reduction in the characteristic path length of 
more than 42%, while the final step reaches a 50% improvement. Even with this evolution strategy, the two 
biggest samples benefit substantially: sample #5 reduces the path from almost 19 to less than 6, while sample 
^9 from more than 11 to 4. Once again with this evolution strategy the saturation tendency in the reduction 
of the path length arises. This tendency is a consequence of the small-world phenomenon that arises when a 
sufficient connectivity threshold is reached. Once this threshold is obtained, the subsequent addition of nodes has 
a reduced effect on the improvement of the property. In this case, we see that this phenomenon arises when the 
first two steps of the growth process are obtained since the characteristic path length then has no more significant 
improvement. The improvement for the clustering coefficient is considerable, but it takes three evolution steps 
before all the samples have non-zero clustering. In the final evolution step the clustering values are all higher 
than 0.1. Robustness almost triples, and eight out of the eleven samples are around or higher 0.7 for this metric. 
The addition of edges according to this strategy is beneficial almost only against the random attacks towards 
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nodes; while the component of the metric that consider targeted attacks against the most connected nodes is 
only partially affected. Considering the values for the redundant path robustness, the same considerations done 
for the characteristic path length apply. There is a reduction of more than 50% compared to the initial value 
when two steps of evolution are completed. The final improvement is of about 63%. 

Table 20 contains for each sample of the Low Voltage Grid (column one) the values of metrics related to 
betweenness. In addition to order and size (columns three and four), average betweenness is provided in columns 
five, while a value of average betweenness normalized by the order of the graph is computed in column six in 
order to compare the different samples. A measure of the statistical variation of betweenness is the coefficient of 
variation which is shown in the seventh column. 

The assortative strategy involving the nodes with lowest node degree is beneficial in reducing the average 
betweenness of all the Low Voltage samples. The benefit that is provided after the four evolution steps is more 
than 60%. Considering the two biggest samples (i.e., samples #5 and #9) the benefits reach up to a value more 
than 70% at the end of the evolution process. The same trend is, of course, followed by the normalized value of 
average betweenness divided by the order of the network. Considering the coefficient of variation for betweenness, 
the general trend for the samples is an increase compared to the original values. Actually in the first step of 
evolution six of the eleven samples experience a decrease in the coefficient of variation, but at the end of the 
process only four keep the reduction trend. This behavior is due to a slower decrease of the standard deviation of 
betweenness compared to its average. This strategy of adding the connectivity provides more connections between 
the nodes that have small connectivity that are usually the nodes at the periphery of the network and therefore 
that do not have paths traversing them in the initial sample. A proof is that the median of the betweenness 
for each sample in the first stages of evolution is zero, then in the last stage of evolution all the samples have a 
non-zero median that is a sign that all nodes are more evenly involved in path between other nodes. 



Triangle closure evolution 



Table 21 contains for each sample of the Low Voltage Grid (column one) the values for the main topological 
quantities: order and size in columns three and four, respectively; average node degree in the fifth column; the 
characteristic path length is reported in column six; the clustering coefficient follows in column seven; robustness 
is shown in eight column; the cost in term of redundant path length closes the data series (column nine) . 

As we already noted in the analysis of the Medium Voltage networks, the triangle closure strategy aims at 
incrementing the clustering coefficient. Therefore it is not surprising that the improvement in the path length 
are below 30% at the final stage of evolution. The increase in the improvement is on average linear between the 
various steps. For instance, the two biggest samples improve only marginally in the characteristic path length 
with values at the end of the 4*'* evolution step of almost 15 and 8 for samples #5 and #8, respectively. The 
evolution of clustering coefficient is the primary goal of the evolution, the values that are obtained in the first 
evolution step are already around 0.1, reaching for the majority of the samples a value near 0.5 in the last 
stage. The improvement for the clustering coefficient is considerable: all the samples after the first step achieve 
a non null coefficient and the two samples that already have a significant clustering coefficient have a ten fold 
improvement. The improvement in robustness is below 70% at the final step of the evolution. The addition 
of edges according to this strategy is beneficial both against random attacks towards nodes and also towards 
targeted attacks against the most connected nodes. Actually, the improvement are much higher in tolerating 
the targeted attacks than the improvement observed against the random failures. Considering the values for 
the redundant path availability, the same considerations done for the characteristic path length apply. There is 
a reduction more than 50% compared to the initial value of the samples when the last stage of edge addition 
is reached, even if the real improvement is already reached after the third stage. However, an exception can 
be noted: sample #2 shows an increase in the value for the redundant path in the first step of the evolution 
compared to the initial value. Such increase is due to the already mentioned problems that affect this metric 
when small networks (in terms of order and size) are considered since they do not have redundant path between 
nodes; such redundant path emerge when more connectivity is added to the network. 



Table 22 contains for each sample of the Low Voltage Grid (column one) the values of metrics related to 
betweenness. In addition to order and size (columns three and four), average betweenness is provided in columns 
five, while a value of average betweenness normalized by the order of the graph is computed in column six in 
order to compare the different samples. A measure of the statistical variation of betweenness is the coefficient of 
variation which is shown in the seventh column. 

The triangle closure strategy provides limited benefits in the reduction of average betweenness and the average 
betweenness to order ratio. The benefits are quantifiable in less than 40% after the whole evolution process. The 
improvements during the several evolution steps progress almost linearly with the addition of edges. Looking at 
the variability of betweenness, only few samples experience a reduction in the coefficient of variation at the end 
of the evolution process. In the first stage, the reduction compared to the original values take place only in six 
out of the eleven samples. The reason is once again the slower decrease of the standard deviation of betweenness 
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compared to its average. The standard deviation, in fact, slowly decreases, with a mode for the betweenness over 
the four evolution stages that is zero for all the samples; the only exception is sample ^1. 



Dissortative node degree evolution 



Table 23 contains for each sample of the Low Voltage Grid (column one) the values for the main topological 
quantities: order and size in columns three and four, respectively; average node degree in the fifth column; the 
characteristic path length is reported in column six; the clustering coefficient follows in column seven; robustness 
is shown in eight column; the cost in term of redundant path length closes the data series (column nine). 

Following the dissortative node degree strategy, the improvement considering the characteristic path length 
increases sub- linearly in the four steps of evolution. The first set of edge addition provides an average improvement 
around 35%, while in the final step the improvement reaches around 54% compared to the initial value. The 
two biggest samples (i.e., samples #5 and #8) reach values of path reduction that are higher than 65%. The 
evolution of clustering coefficient is less uniform than other strategies: the smallest samples in term of order 
reach an high level of improvement (i.e., samples #1, #2, and #3), while for the biggest samples the dissortative 
evolution scheme is less beneficial in the emergence of triangle structures in the network. The improvement in 
robustness is about 77% at the final step of the evolution. The trend in the four evolution steps is almost linear 
with an improvement at each step that is higher than 10%. The addition of edges according to this strategy 
is beneficial mainly against random attacks. The effect towards targeted attacks against the most connected 
nodes improves usually only in the final stages of the evolution (i.e. fourth step). Considering the values for 
the redundant path robustness one notes a double in the path reduction between the first and the last step of 
evolution. For the biggest samples, the difference between the value of characteristic path length and the value 
of 10*'' average path length is around or below three. However, an exception can be noticed: sample #2 shows 
an increase in the value for the redundant path in the first step of the evolution compared to the initial value. 
This increase is due to the problems that affect this metric when small networks (in terms of order and size) 
are considered since they do not have redundant path between nodes; such redundant path emerge when more 
connectivity is added to the network. 



Table 24 contains for each sample of the Low Voltage Grid (column one) the values of metrics related to 
betweenness. In addition to order and size (columns three and four), average betweenness is provided in columns 
five, while a value of average betweenness normalized by the order of the graph is computed in column six in 
order to compare the different samples. A measure of the statistical variation of betweenness is the coefficient of 
variation which is shown in the seventh column. 

The evolution of the betweenness metric realized by the dissortative strategy provides a quite good reduction. 
Overall the reduction is more than 65%, and the major achievement is gained in the first two steps (54% of 
improvement). This trend is followed both by the average betweenness and the normalized value of average 
betweenness divided by the order of the network. Considering the coefficient of variation, all the samples present 
an increase in such a metric. As we have described the average betweenness experiences a consistent reduction, 
while the standard deviation does not follow this trend. For the biggest samples in term of order, the standard 
deviation reduces, while for the smallest ones the standard deviation increases while more edges are added in 
the network. This tendency is due to the very nature of this strategy that aims at connecting different classes 
of nodes in terms of node degree and therefore enabling the creation of new paths between the nodes at the 
edge of the network (having node degree one or two) and the nodes at the core of the network with high node 
degree. Such nodes, previously unaffected by managing the shortest paths, are now involved. This more even 
distribution of shortest paths is shown by the mode of betweenness that already from the second evolution step 
provides non-zero values for the majority of the samples. Later stages show all samples (except the small sample 
^6) with a non-zero mode. 



Random evolution 



Table 25 contains for each sample of the Low Voltage Grid (column one) the values for the main topological 
quantities: order and size in columns three and four, respectively; average node degree in the fifth column; the 
characteristic path length is reported in column six; the clustering coefRcient follows in column seven; robustness 
is shown in eight column; the cost in term of redundant path length closes the data series (column nine). 

The random addition of edges in the network reduces the characteristic path length in the samples up to 50%. 
The first step in the evolution already provides a reduction of 30% compared to the original non-evolved samples. 
The two biggest samples are those which benefits the most of the increase of connectivity with a reduction of 
70% or more. Considering the clustering coefficient, we observe a tendency that we already noted in the Medium 
Voltage network: triangle creation in the small samples is easier and high values for this metric appear, while 
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the bigger samples have smaller values. The first evolution steps still show two samples with null clustering 
coefficient, but after this step all the networks have a significant value for this metric. Robustness almost doubles 
in the whole evolution process compared to the initial values. The addition of edges according to this strategy is 
beneficial both in facing random failures, and also against targeted attacks where actually there is most of the 
improvement. Eight of the eleven samples have values in the last step of the evolution that reach values around 
or higher than 70%. For the significant samples the best improvement is obtained in the second step where on 
average the redundant path length drops of about 50%. The last two steps add smaller benefits in order to reach 
a final reduction about 63%. 

Table 26 contains for each sample of the Medium Voltage Grid (column one) the values of metrics related to 
betweenness. In addition to order and size (columns three and four), average betweenness is provided in columns 
five, while a value of average betweenness normalized by the order of the graph is computed in column six in 
order to compare the different samples. A measure of the statistical variation of betweenness is the coefficient of 
variation which is shown in the seventh column. 

The random evolution strategy provides a consistent reduction in betweenness that reaches 60% of the initial 
value in the final stage of evolution and in the final step all the samples have a value of average betweenness 
to order that is below 3. In addition eight out of the eleven samples have a value that is below two. The 
improvement through the stages of the evolution is sub-linear with already 35% improvement in the first step. 
The random addition of edges provides also a reduction in the variability of betweenness with a coefficient of 
variation that in the last stage of the evolution process is on average 15% lower than in the original samples. In 
addition, eight out of the eleven samples show a coefficient of variation that is below one. Despite this benefits 
in the average betweenness and variation, the mode of betweenness stays zero for almost all samples through all 
the stages of evolution. 
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Sample ID 


Network 
type 


Order 


Size 


Avg. 
deg. 


CPL 


CC 


Removal 
robust- 
ness 

(.RobN) 


Redundancy 
cost 


original order +25% 


1 


MV 


444 


607 


2.734 


6.672 


0.00349 


0.402 


9.997 


2 


MV 


472 


632 


2.678 


6.945 


0.01308 


0.375 


10.496 


3 


MV 


238 


306 


2.571 


6.679 


0.00182 


0.423 


10.775 


4 


MV 


263 


360 


2.738 


6.065 


0.01742 


0.418 


9.945 


5 


MV 


217 


286 


2.636 


6.361 


0.00416 


0.409 


10.384 


6 


MV 


191 


258 


2.702 


5.726 


0.02524 


0.387 


9.448 


7 


MV 


884 


1323 


2.993 


6.557 


0.00711 


0.572 


9.231 


8 


MV 


366 


477 


2.607 


6.919 


0.00000 


0.381 


10.944 


9 


MV 


218 


290 


2.661 


6.014 


0.00282 


0.,387 


9.237 


10 


MV 


201 


255 


2.537 


6.595 


0.00936 


0.433 


11.033 


11 


MV 


202 


266 


2.634 


5.920 


0.00922 


0.408 


9.705 


12 


MV 


464 


623 


2.685 


6.996 


0.00603 


0.374 


10.053 


original order +50% 


1 


MV 


444 


729 


3.284 


5.367 


0.00354 


0.694 


7.872 


2 


MV 


472 


759 


3.216 


5.525 


0.01105 


0.684 


8.243 


3 


MV 


238 


367 


3.084 


5.135 


0.00154 


0.644 


8.086 


4 


MV 


263 


432 


3.285 


4.958 


0.01256 


0.704 


7.587 


5 


MV 


217 


343 


3.161 


4.944 


0.01076 


0.629 


7.701 


6 


MV 


191 


310 


3.246 


4.637 


0.02597 


0.645 


7.277 


7 


MV 


884 


1588 


3.593 


5.472 


0.00876 


0.728 


7.716 


8 


MV 


366 


573 


3.131 


5.519 


0.00209 


0.634 


8.431 


9 


MV 


218 


348 


3.193 


4.714 


0.00249 


0.509 


7.299 


10 


MV 


201 


306 


3.045 


5.120 


0.00843 


0.657 


8.385 


11 


MV 


202 


319 


3.158 


4.808 


0.01047 


0.607 


7.328 


12 


MV 


464 


748 


3.224 


5.711 


0.01216 


0.622 


8.642 


original order +75% 


1 


MV 


444 


850 


3.829 


4.729 


0.00597 


0.757 


6.932 


2 


MV 


472 


885 


3.750 


4.855 


0.01116 


0.736 


7.102 


3 


MV 


238 


428 


3.597 


4.498 


0.00294 


0.702 


6.920 


4 


MV 


263 


504 


3.833 


1.328 


0.01462 


0.719 


6.701 


5 


AiV 


217 


100 


:-!.()iSV 


l.:i()() 


0.0182() 


0.731 


(i.597 


6 


MV 


191 


362 


3.791 


4.089 


0.02379 


0.722 


6.595 


7 


MV 


884 


1853 


4.192 


4.909 


0.00921 


0.777 


6.788 


8 


MV 


366 


668 


3.650 


4.768 


0.00620 


0.743 


7.194 


9 


MV 


218 


406 


3.725 


4.270 


0.01157 


0.710 


6.463 


10 


MV 


201 


357 


3.552 


4.480 


0.01137 


0.744 


6.937 


11 


MV 


202 


372 


3.683 


4.221 


0.00635 


0.698 


6.672 


12 


MV 


464 


873 


3.763 


5.014 


0.01949 


0.748 


7.244 


original order +100% 


1 


MV 


444 


972 


4.378 


4.341 


0.00781 


0.769 


6.345 


2 


MV 


472 


1012 


4.288 


4.468 


0.01186 


0.769 


6.440 


3 


MV 


238 


490 


4.118 


4.025 


0.00504 


0.761 


6.230 


4 


MV 


263 


576 


4.380 


3.966 


0.01835 


0.765 


5.898 


5 


MV 


217 


458 


4.221 


3.949 


0.02150 


0.767 


5.890 


6 


MV 


191 


414 


4.335 


3.726 


0.02659 


0.774 


5.764 


7 


MV 


884 


2118 


4.792 


4.542 


0.01077 


0.786 


6.412 


8 


MV 


366 


764 


4.175 


4.349 


0.00773 


0.753 


6.424 


9 


MV 


218 


464 


4.257 


3.864 


0.01187 


0.740 


5.923 


10 


MV 


201 


408 


4.060 


3.965 


0.01376 


0.764 


6.254 


11 


MV 


202 


426 


4.218 


3.866 


0.00972 


0.760 


6.050 


12 


MV 


464 


998 


4.302 


4.559 


0.03318 


0.773 


6.495 



Table 15: Metrics for random strategy Medium Voltage samples evolution. 
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Sample ID 


Network 


Order 


Size 


Avg. be- 


Avg. 


Coeff. 




type 






tweenness 


betw/order 


vEiria- 














tion 


original order +25% 


1 


iVi V 


444 


607 


zzoy .0 / 


o.uoy 


1 99 


2 


MV 


472 


632 




(\ 07*^ 
O.U 1 o 


1 1 S7 


Q 
O 


MV 


238 


306 


1 909 097 


^ OP^I 


1 1 "iQ 
1. xoy 


4 


MV 


263 


360 


1 ^fil 074 


^ 1 7^1 


1 0fi7 


o 


MV 


217 


286 


1910 "^70 


o.o 1 y 


1 232 




MV 


191 


258 




^.yoo 


1 1 0*^ 
J., lyo 


7 


MV 


884 


1323 


AQA^ fi71 


o.oyu 


1 722 


8 


MV 


366 


477 


2153.055 


5.883 


1.306 


9 


MV 


218 


290 


1130.838 


5.187 


1.201 


10 


MV 


201 


255 


1190.609 


5.923 


1.016 


11 


MV 


202 


266 


981.574 


4.859 


1.145 


12 


.MV 


464 


623 


2897.982 


6.216 


1.204 






original order 


+50% 






1 


MV 


444 


729 


1825.762 


4.112 


1.029 


2 


MV 


472 


759 


2145.374 


4.545 


1.068 


3 


MV 


238 


367 


937.829 


3.94 


1.188 


4 


MV 


263 


432 


1026.759 


3.904 


0.956 


5 


MV 


217 


343 


877.673 


4.045 


1.185 


6 


MV 


191 


310 


696.794 


3.648 


1.094 


7 


MV 


884 


1588 


3985.532 


4.509 


1.455 


8 


MV 


366 


573 


1692.859 


4.625 


1.014 


9 


MV 


218 


348 


819.162 


3.758 


1.204 


10 


MV 


201 


306 


842.223 


4.19 


0.897 


11 


MV 


202 


319 


777.015 


3.847 


0.957 


12 


MV 


464 


748 


2249.061 


4.847 


1.058 


original order +75% 


1 


MV 


444 


850 


1582.139 


3.563 


1.021 


2 


MV 


472 


885 


1830.748 


3.879 


0.919 


3 


MV 


238 


428 


835.813 


3.512 


0.968 


4 


MV 


263 


501 


861.721 


3.288 


0.883 


5 


AIV 


217 


100 


720.953 


3.3.") 


l.()0() 


6 


MV 


191 


362 


599.589 


3.139 


0.972 


7 


MV 


884 


1853 


3465.824 


3.921 


1.275 


8 


MV 


366 


668 


1393.805 


3.808 


0.929 


9 


MV 


218 


406 


719.314 


3.3 


1.082 


10 


MV 


201 


357 


698.497 


3.475 


0.856 


11 


MV 


202 


372 


647.624 


3.206 


0.839 


12 


MV 


464 


873 


1914.675 


4.126 


0.939 


original order +100% 


1 


MV 


444 


972 


1433.794 


3.229 


0.921 


2 


MV 


472 


1012 


1626.893 


3.447 


0.84 


3 


MV 


238 


490 


725.687 


3.049 


0.875 


4 


MV 


263 


576 


752.496 


2.861 


0.842 


5 


MV 


217 


458 


642.514 


2.961 


0.953 


6 


MV 


191 


414 


533.958 


2.796 


0.852 


7 


MV 


884 


2118 


3132.693 


3.544 


1.092 


8 


MV 


366 


764 


1221.086 


3.336 


0.846 


9 


AlV 


218 


16 1 


626.181 


2.872 


0.91 


10 


Ai\" 


201 


108 


r,9 1.832 


2.9."')9 


0.817 


11 


MV 


202 


426 


573.313 


2.838 


0.781 


12 


MV 


464 


998 


1678.895 


3.618 


0.931 



Table 16: Betweenness for random strategy Medium Voltage samples evolution. 
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Sample ID 


Network 
type 


Order 


Size 


Avg. 
deg. 


CPL 


CC 


Removal 
robust- 
ness 
(RobN) 


Redundancy 
cost 


original order +25% 


1 


LV 


17 


22 


2.588 


2.563 


0.08627 


0.534 


7.903 


2 


LV 


15 


18 


2.400 


2.571 


0.08667 


0.379 


8.107 


3 


LV 


21 


27 


2.571 


3.200 


0.06803 


0.533 


8.136 


4 


LV 


24 


28 


2.. 333 


3.348 


0.05000 


0.426 


9.493 


5 


LV 


186 


236 


2.538 


6.286 


0.03886 


0.316 


9.226 


6 


LV 


10 


11 


2.200 


2.222 


0.12000 


0.370 


4.040 


7 


LV 


63 


77 


2.444 


3.839 


0.04743 


0.316 


5.958 


8 


LV 


28 


33 


2.357 


3.074 


0.04464 


0.413 


7.036 


9 


LV 


133 


175 


2.632 


4.750 


0.04700 


0.363 


8.185 


10 


LV 


121 


172 


2.774 


4.581 


0.05247 


0.349 


6.8,35 


11 


L\ 






2.:-!SV 


3,1()T 


0.0ir)88 


0. i<)5 


7.071 


original order -\-50% 


1 


LV 


17 


27 


3.176 


2.188 


0.29483 


0.555 


5.444 


2 


LV 


15 


22 


2.933 


2.429 


0.3,3926 


0.452 


5.589 


3 


LV 


21 


33 


3.143 


2.200 


0.22589 


0.526 


5.518 


4 


LV 


24 


34 


2.833 


2.522 


0.16982 


0.453 


5.931 


5 


LV 


186 


283 


3.043 


6.065 


0.06621 


0.312 


8.256 


6 


LV 


10 


13 


2.600 


2.111 


0.21333 


0.479 


5.360 


7 


LV 


63 


93 


2.952 


3.516 


0.05858 


0.368 


5.691 


8 


LV 


28 


40 


2.857 


2.778 


0.08413 


0.441 


4.709 


9 


LV 


133 


210 


3.158 


4.417 


0.08113 


0.378 


6.123 


10 


LV 


124 


207 


3.339 


4.419 


0.08386 


0.351 


6.188 


11 


LV 


31 


45 


2.903 


2.833 


0.11022 


0.437 


4.954 


original order +75% 


1 


LV 


17 


31 


3.647 


1.875 


0.51569 


0.540 


4.347 


2 


LV 


15 


26 


3.467 


1.929 


0.41657 


0.465 


4.554 


3 


LV 


21 


38 


3.619 


1.950 


0.43312 


0.552 


4.500 


4 


LV 


21 


40 


3.333 


2.370 


0.26865 


0.464 


5.028 


5 


L\ 


18() 


:-!:50 


3.5 W 


5.iS7() 


0.08132 


0.311 


8.210 


6 


LV 


10 


15 


3.000 


1.889 


0.30000 


0.585 


4.000 


7 


LV 


63 


108 


3.429 


2.726 


0.05858 


0.398 


4.828 


8 


LV 


28 


47 


3.357 


2.778 


0.13645 


0.422 


4.464 


9 


LV 


133 


245 


3.684 


4.333 


0.10548 


0.,399 


5.377 


10 


LV 


124 


241 


3.887 


4.346 


0.10261 


0.440 


5.988 


11 


LV 


31 


52 


3.355 


2.633 


0.15801 


0.454 


4.325 


original order +100% 


1 


LV 


17 


36 


4.235 


1.813 


0.58796 


0.493 


3.931 


2 


LV 


15 


30 


4.000 


1.857 


0.51378 


0.529 


4.089 


3 


LV 


21 


44 


4.190 


1.950 


0.47279 


0.532 


4.027 


4 


LV 


24 


46 


3.833 


2.304 


0.29664 


0.508 


4.368 


5 


LV 


186 


378 


4.065 


4.722 


0.12554 


0.343 


6.808 


6 


LV 


10 


18 


3.600 


1.556 


0.30000 


0.730 


3.880 


7 


LV 


63 


124 


3.937 


2.210 


0.05858 


0.408 


4.149 


8 


LV 


28 


54 


3.857 


2.685 


0.2,3988 


0.471 


4.214 


9 


LV 


133 


280 


4.211 


4.167 


0.13065 


0.397 


5.543 


10 


LV 


124 


276 


4.452 


4.191 


0.12194 


0.439 


5.465 


11 


LV 


31 


60 


3.871 


2.567 


0.19894 


0.438 


4.217 



Table 17: Metrics for assortative high node degree strategy Low Voltage samples evolution. 
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Sample ID 


Network 
type 


Order 


Size 


Avg. be- 
tweenness 


Avg. 

betw/order 


Coeff. 
varia- 
tion 


original order +25% 


1 


LV 


17 


22 


32.933 


1.937 


0.773 


2 


LV 


15 


18 


25.231 


1.682 


0.887 


3 


LV 


21 


27 


44.571 


2.122 


1.041 


4 


LV 


24 


28 


57 


2.375 


1.346 


5 


LV 


186 


236 


940.227 


5.055 


1.661 


6 


LV 


10 


11 


9 


0.9 


1.291 


7 


LV 


63 


77 


161.607 


2.565 


2.007 


8 


LV 


28 


33 


55.929 


1.997 


1.935 


9 


LV 


133 


175 


481.488 


3.62 


2.044 


10 


LV 


121 


172 


415.696 


3.352 


1.894 


11 


L\" 


31 


37 


():i.871 


2.()() 


1.822 


original order +50% 


1 


LV 


17 


27 


32.933 


1.937 


0.773 


2 


LV 


15 


22 


25.231 


1.682 


0.887 


3 


LV 


21 


33 


27.524 


1.311 


2.194 


4 


LV 


24 


34 


38.667 


1.611 


1.691 


5 


LV 


186 


283 


904.648 


4.864 


1.513 


6 


LV 


10 


13 


8 


0.8 


1.173 


7 


LV 


63 


93 


151.065 


2.398 


2.158 


8 


LV 


28 


40 


50 


1.786 


1.589 


9 


LV 


133 


210 


423.442 


3.184 


1.842 


10 


LV 


124 


207 


396.786 


3.2 


1.645 


11 


LV 


31 


45 


53.419 


1.723 


1.734 


original order +75% 


1 


LV 


17 


31 


32.933 


1.937 


0.773 


2 


LV 


15 


26 


22.308 


1.487 


1.094 


3 


LV 


21 


38 


20.381 


0.971 


3.138 


4 


L\' 


24 


40 


35.5 


1.479 


1.751 


5 


LV 


180 


330 


8H8. <)()() 


1.779 


1. 181 


6 


LV 


10 


15 


6 


0.6 


1 


7 


LV 


63 


108 


114.903 


1.824 


2.81 


8 


LV 


28 


47 


44.071 


1.574 


1.636 


9 


LV 


133 


245 


407.922 


3.067 


1.684 


10 


LV 


124 


241 


389.571 


3.142 


1.575 


11 


LV 


31 


52 


49.935 


1.611 


1.511 


original order +100% 


1 


LV 


17 


36 


20.533 


1.208 


0.83 


2 


LV 


15 


30 


18.923 


1.262 


1.434 


3 


LV 


21 


44 


19.619 


0.934 


2.66 


4 


LV 


24 


46 


33.5 


1.396 


1.409 


5 


LV 


186 


378 


882.58 


4.745 


1.482 


6 


LV 


10 


18 


3.75 


0.375 


0.702 


7 


LV 


63 


124 


76.548 


1.215 


4.604 


8 


LV 


28 


54 


40.643 


1.452 


1.692 


9 


LV 


133 


280 


395.039 


2.97 


1.66 


10 


L\' 


121 


276 


383.45 


3.092 


1.588 


11 


L\" 


31 


()0 


15.935 


1.182 


1. i3() 



Table 18: Betweenness for assortative high node degree strategy Low Vohage samples evolution. 
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Sample ID 


Network 
type 


Order 


Size 


Avg. 
deg. 


CPL 


CC 


Removal 
robust- 
ness 
(RobN) 


Redundancy 
cost 


original order +25% 


1 


LV 


17 


22 


2.588 


2.813 


0.15490 


0.473 


8.000 


2 


LV 


15 


18 


2.400 


3.000 


0.06667 


0.528 


10.179 


3 


LV 


21 


27 


2.571 


3.950 


0.19955 


0.529 


8.209 


4 


LV 


24 


28 


2.. 333 


3.239 


0.00000 


0.438 


10.868 


5 


LV 


186 


236 


2.538 


5.638 


0.00701 


0.318 


9.589 


6 


LV 


10 


11 


2.200 


2.111 


0.00000 


0.391 


5.760 


7 


LV 


63 


77 


2.444 


4.790 


0.13086 


0.275 


7.416 


8 


LV 


28 


33 


2.357 


3.704 


0.04932 


0.362 


9.301 


9 


LV 


133 


175 


2.632 


4.864 


0.01128 


0.356 


9.266 


10 


LV 


121 


172 


2.774 


4.622 


0.00899 


0.374 


8.052 


11 


L\ 






2.:-!SV 


:-!,S00 


0.00000 


0. ill) 


10,871 


original order -\-50% 


1 


LV 


17 


27 


3.176 


2.125 


0.36061 


0.517 


5.028 


2 


LV 


15 


22 


2.933 


2.500 


0.23111 


0.551 


5.982 


3 


LV 


21 


33 


3.143 


2.200 


0.25763 


0.604 


5.309 


4 


LV 


24 


34 


2.833 


2.804 


0.07262 


0.453 


6.417 


5 


LV 


186 


283 


3.043 


5.622 


0.08388 


0.371 


7.650 


6 


LV 


10 


13 


2.600 


2.000 


0.00000 


0.410 


6.640 


7 


LV 


63 


93 


2.952 


3.339 


0.13006 


0.329 


6.233 


8 


LV 


28 


40 


2.857 


2.630 


0.04801 


0.391 


6.760 


9 


LV 


133 


210 


3.158 


4.083 


0.08720 


0.415 


6.604 


10 


LV 


124 


207 


3.339 


4.073 


0.07965 


0.404 


6.253 


11 


LV 


31 


45 


2.903 


2.433 


0.04328 


0.460 


6.533 


original order +75% 


1 


LV 


17 


31 


3.647 


1.875 


0.58095 


0.519 


4.222 


2 


LV 


15 


26 


3.467 


2.000 


0.40296 


0.575 


4.732 


3 


LV 


21 


38 


3.619 


1.950 


0.42519 


0.610 


4.409 


4 


LV 


21 


40 


3.333 


2.739 


0.14484 


0.522 


5.326 


5 


L\ 


18() 


:-!:50 


:-!.5i<s 


5.(il() 


o.iiijyo 


0. i0(i 


(i.828 


6 


LV 


10 


15 


3.000 


1.778 


0.12000 


0.629 


5.360 


7 


LV 


63 


108 


3.429 


2.226 


0.12992 


0.374 


4.475 


8 


LV 


28 


47 


3.357 


2.426 


0.09464 


0.501 


4.959 


9 


LV 


133 


245 


3.684 


4.045 


0.14127 


0.505 


5.632 


10 


LV 


124 


241 


3.887 


4.041 


0.14584 


0.420 


5.958 


11 


LV 


31 


52 


3.355 


2.400 


0.14910 


0.534 


5.158 


original order +100% 


1 


LV 


17 


36 


4.235 


1.875 


0.59314 


0.544 


3.889 


2 


LV 


15 


30 


4.000 


1.857 


0.51091 


0.591 


4.179 


3 


LV 


21 


44 


4.190 


1.950 


0.44977 


0.628 


3.991 


4 


LV 


24 


46 


3.833 


2.674 


0.20079 


0.769 


5.090 


5 


LV 


186 


378 


4.065 


5.597 


0.14059 


0.404 


6.622 


6 


LV 


10 


18 


3.600 


1.556 


0.25000 


0.688 


4.320 


7 


LV 


63 


124 


3.937 


2.194 


0.16174 


0.394 


4.201 


8 


LV 


28 


54 


3.857 


2.389 


0.17619 


0.514 


4.577 


9 


LV 


133 


280 


4.211 


4.008 


0.17196 


0.508 


5.175 


10 


LV 


124 


276 


4.452 


3.992 


0.18500 


0.512 


5.460 


11 


LV 


31 


60 


3.871 


2.333 


0.20753 


0.508 


4.533 



Table 19: Metrics for assortative low node degree strategy Low Voltage samples evolution. 
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Sample ID 


Network 
type 


Order 


Size 


Avg. be- 
tweenness 


Avg. 

betw/order 


Coeff. 
varia- 
tion 


original order +25% 


1 


LV 


17 


22 


24.933 


1.467 


0.978 


2 


LV 


15 


18 


27.429 


1.829 


0.832 


3 


LV 


21 


27 


61.048 


2.907 


0.854 


4 


LV 


24 


28 


53.5 


2.229 


0.947 


5 


LV 


186 


236 


888.486 


4.777 


2.442 


6 


LV 


10 


11 


10.667 


1.067 


1.021 


7 


LV 


63 


77 


228.984 


3.635 


2.02 


8 


LV 


28 


33 


79.143 


2.827 


1.102 


9 


LV 


133 


175 


562.931 


4.233 


2.497 


10 


LV 


121 


172 


465.817 


3.757 


1.985 


11 


L\" 


31 




87. (il:-) 


2.820 


1.113 


original order +50% 


1 


LV 


17 


1 


14.267 


0.839 


2.177 


2 


LV 


15 


2 


21.429 


1.429 


0.898 


3 


LV 


21 


3 


26.762 


1.274 


2.128 


4 


LV 


24 


4 


42.583 


1.774 


0.912 


5 


LV 


186 


5 


886.077 


4.764 


1.781 


6 


LV 


10 


6 


6.889 


0.689 


0.867 


7 


LV 


63 


7 


163.377 


2.593 


2.184 


8 


LV 


28 


8 


47.214 


1.686 


1.696 


9 


LV 


133 


9 


414.901 


3.12 


2.379 


10 


LV 


124 


10 


404.085 


3.259 


1.946 


11 


LV 


31 


11 


49.484 


1.596 


2.238 


original order +75% 


1 


LV 


17 


1 


10.267 


0.604 


3.352 


2 


LV 


15 


2 


13 


0.867 


2.042 


3 


LV 


21 


3 


20.286 


0.966 


3.114 


4 


L\' 


24 


4 


.39.5 


1.646 


0.818 


5 


LV 


180 


o 


883.901 


i.752 


1. 197 


6 


LV 


10 


6 


5.556 


0.556 


0.652 


7 


LV 


63 


7 


83.645 


1.328 


4.621 


8 


LV 


28 


8 


41.857 


1.495 


1.636 


9 


LV 


133 


9 


409.145 


3.076 


1.912 


10 


LV 


124 


10 


402.305 


3.244 


1.627 


11 


LV 


31 


11 


47.742 


1.54 


1.977 


original order +100% 


1 


LV 


17 


1 


9.6 


0.565 


2.951 


2 


LV 


15 


2 


11.286 


0.752 


2.184 


3 


LV 


21 


3 


19.524 


0.93 


2.774 


4 


LV 


24 


4 


37.417 


1.559 


0.725 


5 


LV 


186 


5 


879.801 


4.73 


1.323 


6 


LV 


10 


6 


4.444 


0.444 


0.743 


7 


LV 


63 


7 


79.839 


1.267 


4.655 


8 


LV 


28 


8 


41.143 


1.469 


1.275 


9 


LV 


133 


9 


399.252 


3.002 


1.695 


10 


LV 


124 


10 


400.576 


3.23 


1.447 


11 


LV 


31 


11 


45.548 


1.469 


1.559 



Table 20: Betweenness for assortative low node degree strategy Low Voltage samples evolution. 
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Sample ID 


Network 
type 


Order 


Size 


Avg. 
deg. 


CPL 


cc 


Removal 
robust- 
ness 

(RobN) 


Redundancy 
cost 


original order +25% 


1 


LV 


17 


22 


2.588 


3.063 


0.11765 


0.489 


7.722 


2 


LV 


15 


18 


2.400 


2.857 


0.10667 


0.513 


8.000 


3 


LV 


21 


27 


2.571 


3.850 


0.13492 


0.484 


7.709 


4 


LV 


24 


28 


2.. 333 


4.087 


0.09583 


0.367 


6.611 


5 


LV 


186 


236 


2.538 


17.443 


0.08734 


0.142 


24.535 


6 


LV 


10 


11 


2.200 


2.056 


0.21333 


0.350 


4.280 


7 


LV 


63 


77 


2.444 


5.323 


0.13171 


0.249 


6.188 


8 


LV 


28 


33 


2.357 


4.222 


0.21548 


0.359 


5.607 


9 


LV 


133 


175 


2.632 


10.174 


0.09291 


0.232 


12.059 


10 


LV 


121 


172 


2.774 


6.943 


0.08714 


0.280 


9.905 


11 


L\ 






2.:-!SV 


l,<)(iT 


0.1.') 181 


().:-! 10 


().(i88 


original order -\-50% 


1 


LV 


17 


27 


3.176 


2.625 


0.22941 


0.500 


5.194 


2 


LV 


15 


22 


2.933 


2.643 


0.20889 


0.524 


5.. 304 


3 


LV 


21 


33 


3.143 


3.500 


0.33492 


0.500 


6.445 


4 


LV 


24 


34 


2.833 


3.630 


0.17817 


0.407 


5.382 


5 


LV 


186 


283 


3.043 


16.557 


0.15282 


0.161 


20.545 


6 


LV 


10 


13 


2.600 


2.000 


0.35667 


0.378 


5.400 


7 


LV 


63 


93 


2.952 


5.274 


0.14160 


0.282 


5.962 


8 


LV 


28 


40 


2.857 


3.815 


0.44643 


0.328 


6.209 


9 


LV 


133 


210 


3.158 


9.295 


0.21776 


0.345 


12.848 


10 


LV 


124 


207 


3.339 


6.683 


0.18182 


0.428 


9.559 


11 


LV 


31 


45 


2.903 


4.067 


0.29724 


0.355 


5.917 


original order +75% 


1 


LV 


17 


31 


3.647 


2.500 


0.30644 


0.581 


4.625 


2 


LV 


15 


26 


3.467 


2.357 


0.. 34222 


0.598 


4.554 


3 


LV 


21 


38 


3.619 


2.950 


0.41882 


0.549 


5.200 


4 


LV 


21 


40 


3.333 


3.348 


0.37718 


0.433 


5.069 


5 


L\ 


18() 


:-!:50 


3.5 W 


1.'). iOO 


0.2917:5 


0.23(i 


17.191 


6 


LV 


10 


15 


3.000 


1.889 


0.42333 


0.504 


4.440 


7 


LV 


63 


108 


3.429 


5.161 


0.20174 


0.487 


5.800 


8 


LV 


28 


47 


3.357 


3.648 


0.52262 


0.416 


5.643 


9 


LV 


133 


245 


3.684 


8.265 


0.. 34422 


0.375 


10.834 


10 


LV 


124 


241 


3.887 


5.972 


0.28127 


0.416 


7.653 


11 


LV 


31 


52 


3.355 


3.733 


0.38771 


0.424 


5.388 


original order +100% 


1 


LV 


17 


36 


4.235 


2.125 


0.41167 


0.596 


4.097 


2 


LV 


15 


30 


4.000 


2.143 


0.48937 


0.669 


4.018 


3 


LV 


21 


44 


4.190 


2.750 


0.56485 


0.565 


4.882 


4 


LV 


24 


46 


3.833 


3.043 


0.48585 


0.458 


4.660 


5 


LV 


186 


378 


4.065 


14.792 


0.44891 


0.250 


16.530 


6 


LV 


10 


18 


3.600 


1.667 


0.48333 


0.673 


3.720 


7 


LV 


63 


124 


3.937 


5.065 


0.29045 


0.550 


6.058 


8 


LV 


28 


54 


3.857 


3.426 


0.64345 


0.406 


5.230 


9 


LV 


133 


280 


4.211 


8.015 


0.39996 


0.383 


11.030 


10 


LV 


124 


276 


4.452 


5.224 


0.33963 


0.526 


6.988 


11 


LV 


31 


60 


3.871 


3.300 


0.51951 


0.391 


5.092 



Table 21: Metrics for triangle closure strategy Low Voltage samples evolution. 
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Sample ID 


Network 
type 


Order 


Size 


Avg. be- 
tweenness 


Avg. 

betw/order 


Coeff. 
varia- 
tion 


original order +25% 


1 


LV 


17 


22 


28.4 


1.671 


0.649 


2 


LV 


15 


18 


22.154 


1.477 


0.646 


3 


LV 


21 


27 


60.476 


2.88 


0.561 


4 


LV 


24 


28 


71.917 


2.997 


0.967 


5 


LV 


186 


236 


2821.193 


15.168 


1.214 


6 


LV 


10 


11 


7 


0.7 


1.429 


7 


LV 


63 


77 


252.885 


4.014 


2.074 


8 


LV 


28 


33 


88.786 


3.171 


1.475 


9 


LV 


133 


175 


1217.519 


9.154 


1.495 


10 


LV 


121 


172 


678.804 


5.471 


1.395 


11 


L\" 


31 


37 


120.258 


:-!.87<) 


1.313 


original order +50% 


1 


LV 


17 


1 


24.8 


1.459 


0.633 


2 


LV 


15 


2 


18.923 


1.262 


0.715 


3 


LV 


21 


3 


53.143 


2.531 


0.621 


4 


LV 


24 


4 


63.833 


2.66 


0.982 


5 


LV 


186 


5 


2747.807 


14.773 


1.25 


6 


LV 


10 


6 


5.5 


0.55 


1.226 


7 


LV 


63 


7 


247.443 


3.928 


2.111 


8 


LV 


28 


8 


75.357 


2.691 


1.704 


9 


LV 


133 


9 


1096.527 


8.245 


1.594 


10 


LV 


124 


10 


653.232 


5.268 


1.424 


11 


LV 


31 


11 


94.258 


3.041 


1.46 


original order +75% 


1 


LV 


17 


1 


21.333 


1.255 


0.596 


2 


LV 


15 


2 


17.846 


1.19 


0.731 


3 


LV 


21 


3 


39.905 


1.9 


0.898 


4 


L\' 


24 


4 


55.25 


2.302 


1.03 


5 


LV 


180 


o 


2513.07 


13.511 


1.278 


6 


LV 


10 


6 


5 


0.5 


1.153 


7 


LV 


63 


7 


238.754 


3.79 


2.201 


8 


LV 


28 


8 


70.143 


2.505 


1.801 


9 


LV 


133 


9 


992.558 


7.463 


1.703 


10 


LV 


124 


10 


571.214 


4.607 


1.541 


11 


LV 


31 


11 


84.581 


2.728 


1.483 


original order +100% 


1 


LV 


17 


1 


16.933 


0.996 


0.972 


2 


LV 


15 


2 


13.077 


0.872 


0.85 


3 


LV 


21 


3 


35.429 


1.687 


0.874 


4 


LV 


24 


4 


44.833 


1.868 


1.373 


5 


LV 


186 


5 


2375.205 


12.77 


1.34 


6 


LV 


10 


6 


4.25 


0.425 


1.157 


7 


LV 


63 


7 


234.885 


3.728 


2.202 


8 


LV 


28 


8 


62.071 


2.217 


1.967 


9 


LV 


133 


9 


951.163 


7.152 


1.716 


10 


L\' 


121 


10 


509.321 


1.107 


1.603 


11 


L\" 


31 


11 


7:-!. 11!) 


2.3()8 


l.()59 



Table 22: Betweenncss for triangle closure strategy Low Voltage samples evolution. 
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Sample ID 


Network 
type 


Order 


Size 


Avg. 
deg. 


CPL 


CC 


Removal 
robust- 
ness 

(RobN) 


Redundancy 
cost 


original order +25% 


1 


LV 


17 


22 


2.588 


2.563 


0.11373 


0.539 


8.111 


2 


LV 


15 


18 


2.400 


2.571 


0.27143 


0.510 


9.214 


3 


LV 


21 


27 


2.571 


2.900 


0.05873 


0.440 


8.427 


4 


LV 


24 


28 


2.. 333 


3.326 


0.08426 


0.396 


7.313 


5 


LV 


186 


236 


2.538 


5.354 


0.02151 


0.356 


9.082 


6 


LV 


10 


11 


2.200 


2.000 


0.40952 


0.352 


4.360 


7 


LV 


63 


77 


2.444 


3.484 


0.00000 


0.319 


7.449 


8 


LV 


28 


33 


2.357 


3.241 


0.04827 


0.375 


9.204 


9 


LV 


133 


175 


2.632 


4.152 


0.01128 


0.372 


8.367 


10 


LV 


121 


172 


2.774 


4.195 


0.00899 


0.389 


7.962 


11 


L\ 






2.:-!SV 


3,5()T 


0.00000 


0. :-!<)() 


7. (Hi:-! 


original order +50% 


1 


LV 


17 


27 


3.176 


2.250 


0.31634 


0.519 


5.417 


2 


LV 


15 


22 


2.933 


2.500 


0.25873 


0.509 


5.893 


3 


LV 


21 


33 


3.143 


2.450 


0.28903 


0.541 


5.936 


4 


LV 


24 


34 


2.833 


2.326 


0.14274 


0.4.38 


6.333 


5 


LV 


186 


283 


3.043 


4.897 


0.02869 


0.362 


7.066 


6 


LV 


10 


13 


2.600 


1.778 


0.81111 


0.368 


6.520 


7 


LV 


63 


93 


2.952 


2.435 


0.02118 


0.365 


5.843 


8 


LV 


28 


40 


2.857 


2.426 


0.05030 


0.413 


6.230 


9 


LV 


133 


210 


3.158 


3.803 


0.01128 


0.413 


5.593 


10 


LV 


124 


207 


3.339 


3.951 


0.01976 


0.405 


6.717 


11 


LV 


31 


45 


2.903 


2.267 


0.00000 


0.448 


6.304 


original order +75% 


1 


LV 


17 


31 


3.647 


1.875 


0.51036 


0.574 


4.333 


2 


LV 


15 


26 


3.467 


2.214 


0.19222 


0.647 


4.446 


3 


LV 


21 


38 


3.619 


2.050 


0.50913 


0.489 


4.627 


4 


LV 


21 


40 


3.333 


2.304 


0.11278 


0. 176 


4.757 


5 


L\ 


18() 


:-!:50 


3.5 W 


1.502 


0.02511 


0. 109 


(i.G.jl 


6 


LV 


10 


15 


3.000 


1.778 


0.83333 


0.385 


5.040 


7 


LV 


63 


108 


3.429 


2.161 


0.03188 


0.388 


4.682 


8 


LV 


28 


47 


3.357 


2.315 


0.04837 


0.435 


4.418 


9 


LV 


133 


245 


3.684 


3.712 


0.05067 


0.443 


5.044 


10 


LV 


124 


241 


3.887 


3.569 


0.01706 


0.413 


5.372 


11 


LV 


31 


52 


3.355 


2.233 


0.02222 


0.457 


4.796 


original order +100% 


1 


LV 


17 


36 


4.235 


1.813 


0.59804 


0.580 


3.903 


2 


LV 


15 


30 


4.000 


1.857 


0.25603 


0.646 


4.036 


3 


LV 


21 


44 


4.190 


1.850 


0.62824 


0.490 


3.964 


4 


LV 


24 


46 


3.833 


2.283 


0.15157 


0.474 


4.431 


5 


LV 


186 


378 


4.065 


4.362 


0.03534 


0.409 


7.533 


6 


LV 


10 


18 


3.600 


1.667 


0.69167 


0.487 


4.000 


7 


LV 


63 


124 


3.937 


2.161 


0.02121 


0.402 


4.147 


8 


LV 


28 


54 


3.857 


2.130 


0.03622 


0.485 


4.276 


9 


LV 


133 


280 


4.211 


3.652 


0.07177 


0.477 


5.164 


10 


LV 


124 


276 


4.452 


3.301 


0.03482 


0.510 


4.964 


11 


LV 


31 


60 


3.871 


2.167 


0.02174 


0.552 


4.383 



Table 23: Metrics for dissortative node degree strategy Low Voltage saraples evolution. 
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Sample ID 


Network 
type 


Order 


Size 


Avg. be- 
tweenness 


Avg. 

betw/order 


Coeff. 
varia- 
tion 


original order +25% 


1 


LV 


17 


22 


26 


1.529 


0.99 


2 


LV 


15 


18 


22.143 


1.476 


1.315 


3 


LV 


21 


27 


40.571 


1.932 


1.204 


4 


LV 


24 


28 


55.25 


2.302 


1.617 


5 


LV 


186 


236 


843.381 


4.534 


2.559 


6 


LV 


10 


11 


6 


0.6 


2 


7 


LV 


63 


77 


180.129 


2.859 


2.143 


8 


LV 


28 


33 


65.643 


2.344 


1.31 


9 


LV 


133 


175 


431.169 


3.242 


3.003 


10 


LV 


121 


172 


393.897 


3.177 


2.39 


11 


L\" 


31 


37 


T9.0:-!2 


2..') 11) 


1.808 


original order +50% 


1 


LV 


17 


1 


19 


1.118 


1.57 


2 


LV 


15 


2 


20.429 


1.362 


1.11 


3 


LV 


21 


3 


31.714 


1.51 


1.788 


4 


LV 


24 


4 


36.583 


1.524 


2.225 


5 


LV 


186 


5 


800.796 


4.305 


1.92 


6 


LV 


10 


6 


5.556 


0.556 


2.828 


7 


LV 


63 


7 


102.871 


1.633 


4.007 


8 


LV 


28 


8 


40.571 


1.449 


2.307 


9 


LV 


133 


9 


362.397 


2.725 


2.76 


10 


LV 


124 


10 


377.458 


3.044 


2.045 


11 


LV 


31 


11 


43.355 


1.399 


2.934 


original order +75% 


1 


LV 


17 


1 


13.375 


0.787 


2.611 


2 


LV 


15 


2 


15 


1 


1.021 


3 


LV 


21 


3 


22 


1.048 


2.876 


4 


L\' 


24 


4 


33.833 


1.41 


1.751 


5 


LV 


180 


o 


737.318 


3.9(il 


1.751 


6 


LV 


10 


6 


5.556 


0.556 


2.828 


7 


LV 


63 


7 


77.226 


1.226 


4.985 


8 


LV 


28 


8 


36.857 


1.316 


2.114 


9 


LV 


133 


9 


349.13 


2.625 


2.351 


10 


LV 


124 


10 


328.119 


2.646 


2.069 


11 


LV 


31 


11 


41.613 


1.342 


2.69 


original order +100% 


1 


LV 


17 


1 


11.125 


0.654 


3.411 


2 


LV 


15 


2 


12.429 


0.829 


1.127 


3 


LV 


21 


3 


15.81 


0.753 


3.915 


4 


LV 


24 


4 


32.167 


1.34 


1.523 


5 


LV 


186 


5 


705.503 


3.793 


1.659 


6 


LV 


10 


6 


5.333 


0.533 


2.763 


7 


LV 


63 


7 


74.871 


1.188 


4.318 


8 


LV 


28 


8 


33 


1.179 


2.256 


9 


LV 


133 


9 


338.931 


2.548 


2.14 


10 


L\' 


121 


10 


271.627 


2.215 


2.217 


11 


L\" 


31 


11 


39.220 


1.2()5 


2.102 



Table 24: Betweenness for dissortative node degree strategy Low Voltage samples evolution. 
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Sample ID 


Network 
type 


Order 


Size 


Avg. 
deg. 


CPL 


CC 


Removal 
robust- 
ness 
(RobN) 


Redundancy 
cost 


original order +25% 


1 


LV 


17 


22 


2.588 


2.563 


0.12157 


0.541 


7.556 


2 


LV 


15 


18 


2.400 


2.786 


0.00000 


0.530 


9.679 


3 


LV 


21 


27 


2.571 


3.150 


0.00000 


0.422 


9.036 


4 


LV 


24 


28 


2.. 333 


3.413 


0.06250 


0.487 


9.979 


5 


LV 


186 


236 


2.538 


6.038 


0.00269 


0.402 


9.931 


6 


LV 


10 


11 


2.200 


2.222 


0.06000 


0.365 


5.040 


7 


LV 


63 


77 


2.444 


3.823 


0.04321 


0.352 


8.029 


8 


LV 


28 


33 


2.357 


3.759 


0.08095 


0.445 


7.408 


9 


LV 


133 


175 


2.632 


5.205 


0.01078 


0.435 


8.984 


10 


LV 


121 


172 


2.774 


4.805 


0.0.3015 


0.407 


8.297 


11 


L\ 






2.:-!SV 


3, 1()T 


0.0118:5 


0. KiO 


S.5()<s 


original order -\-50% 


1 


LV 


17 


27 


3.176 


2.375 


0.29412 


0.677 


5.181 


2 


LV 


15 


22 


2.933 


2.357 


0.05556 


0.743 


5.964 


3 


LV 


21 


33 


3.143 


2.500 


0.01587 


0.598 


5.609 


4 


LV 


24 


34 


2.833 


2.652 


0.05278 


0.613 


6.347 


5 


LV 


186 


283 


3.043 


4.830 


0.00520 


0.546 


8.064 


6 


LV 


10 


13 


2.600 


2.056 


0.34667 


0.441 


6.080 


7 


LV 


63 


93 


2.952 


3.435 


0.04780 


0.403 


6.180 


8 


LV 


28 


40 


2.857 


3.481 


0.16378 


0.438 


5.903 


9 


LV 


133 


210 


3.158 


4.318 


0.02016 


0.683 


7.078 


10 


LV 


124 


207 


3.339 


4.045 


0.02746 


0.580 


6.771 


11 


LV 


31 


45 


2.903 


2.833 


0.01183 


0.690 


6.379 


original order +75% 


1 


LV 


17 


31 


3.647 


2.188 


0.34146 


0.745 


4.528 


2 


LV 


15 


26 


3.467 


2.071 


0.17778 


0.689 


4.786 


3 


LV 


21 


38 


3.619 


2.300 


0.12748 


0.633 


4.564 


4 


LV 


21 


40 


3.333 


2.457 


0.15615 


0.621 


5.347 


5 


L\ 


18() 


:-!:50 


3.5 W 


1.251 


O.OllK) 


0.700 


(i.807 


6 


LV 


10 


15 


3.000 


1.889 


0.28667 


0.508 


5.480 


7 


LV 


63 


108 


3.429 


3.081 


0.09340 


0.438 


5.438 


8 


LV 


28 


47 


3.357 


2.870 


0.22636 


0.528 


5.541 


9 


LV 


133 


245 


3.684 


3.886 


0.03401 


0.715 


6.126 


10 


LV 


124 


241 


3.887 


3.589 


0.03719 


0.706 


5.783 


11 


LV 


31 


52 


3.355 


2.667 


0.06237 


0.712 


5.554 


original order +100% 


1 


LV 


17 


36 


4.235 


2.063 


0.38277 


0.750 


4.083 


2 


LV 


15 


30 


4.000 


1.929 


0.14000 


0.722 


4.214 


3 


LV 


21 


44 


4.190 


2.100 


0.12611 


0.635 


4.255 


4 


LV 


24 


46 


3.833 


2.283 


0.15258 


0.744 


4.757 


5 


LV 


186 


378 


4.065 


3.886 


0.02847 


0.754 


6.114 


6 


LV 


10 


18 


3.600 


1.722 


0.31667 


0.737 


4.320 


7 


LV 


63 


124 


3.937 


2.935 


0.10115 


0.451 


5.012 


8 


LV 


28 


54 


3.857 


2.556 


0.25003 


0.620 


4.679 


9 


LV 


133 


280 


4.211 


3.538 


0.04299 


0.731 


5.646 


10 


LV 


124 


276 


4.452 


3.289 


0.04395 


0.744 


5.179 


11 


LV 


31 


60 


3.871 


2.533 


0.12734 


0.696 


4.888 



Table 25: Metrics for random strategy Low Voltage samples evolution. 
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Sample ID 


Network 
type 


Order 


Size 


Avg. be- 
tweenness 


Avg. 

betw/order 


Coeff. 
varia- 
tion 


original order +25% 


1 


LV 


17 


22 


20.667 


1.216 


0.78 


2 


LV 


15 


18 


23.429 


1.562 


0.854 


3 


LV 


21 


27 


42.19 


2.009 


1.256 


4 


LV 


24 


28 


59.833 


2.493 


0.841 


5 


LV 


186 


236 


898.716 


4.832 


1.213 


6 


LV 


10 


11 


8 


0.8 


1.256 


7 


LV 


63 


77 


191.129 


3.034 


1.724 


8 


LV 


28 


33 


82 


2.929 


1.263 


9 


LV 


133 


175 


558.279 


4.198 


1.195 


10 


LV 


121 


172 


452.318 


3.648 


1.116 


11 


L\" 


31 


37 


T5.1()l 


2.125 


1.22() 


original order +50% 


1 


LV 


17 


1 


16.933 


0.996 


0.858 


2 


LV 


15 


2 


17.857 


1.19 


0.716 


3 


LV 


21 


3 


32 


1.524 


0.968 


4 


LV 


24 


4 


41.833 


1.743 


0.949 


5 


LV 


186 


5 


693.125 


3.726 


0.923 


6 


LV 


10 


6 


8.889 


0.889 


1.497 


7 


LV 


63 


7 


150.548 


2.39 


1.469 


8 


LV 


28 


8 


72.286 


2.582 


1.384 


9 


LV 


133 


9 


431.308 


3.243 


0.916 


10 


LV 


124 


10 


375.254 


3.026 


1.065 


11 


LV 


31 


11 


60 


1.935 


0.85 


original order +75% 


1 


LV 


17 


1 


14.533 


0.855 


0.899 


2 


LV 


15 


2 


14.857 


0.99 


0.764 


3 


LV 


21 


3 


26.571 


1.265 


1.087 


4 


L\' 


24 


4 


37 


1.542 


1.061 


5 


LV 


180 


o 


()03. 138 


3.241 


0.859 


6 


LV 


10 


6 


7.778 


0.778 


1.301 


7 


LV 


63 


7 


129.71 


2.059 


1.528 


8 


LV 


28 


8 


56.786 


2.028 


1.061 


9 


LV 


133 


9 


381.847 


2.871 


0.915 


10 


LV 


124 


10 


324.644 


2.618 


1.048 


11 


LV 


31 


11 


51.032 


1.646 


0.899 


original order +100% 


1 


LV 


17 


1 


12.4 


0.729 


0.97 


2 


LV 


15 


2 


12.143 


0.81 


0.646 


3 


LV 


21 


3 


23.238 


1.107 


0.979 


4 


LV 


24 


4 


31.75 


1.323 


0.808 


5 


LV 


186 


5 


534.156 


2.872 


0.812 


6 


LV 


10 


6 


6.667 


0.667 


1.095 


7 


LV 


63 


7 


120.097 


1.906 


1.442 


8 


LV 


28 


8 


42.5 


1.518 


1.269 


9 


LV 


133 


9 


339.893 


2.556 


0.86 


10 


L\' 


121 


10 


288.251 


2.325 


0.989 


11 


L\" 


31 


11 


15.181 


1.1()7 


0.891 



Table 26: Betweenness for random strategy Low Voltage samples evolution. 
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6 Discussion 



To determine which strategy, from a topological point of view, is more rewarding and provides the higher 
benefits one needs to compare the results presented in Section [sj If determining a "best" strategy is unfeasible, 
one should aim at identifying trade-offs between one strategy compared to another one. We evaluate all the 
strategies considering average values over all the samples and considering the satisfaction of the topological 
metrics required for the Smart Grid explained in |59] and summarized in Appendix A We consider, as we have 
done in the presentation of results, the Medium Voltage Grid and the Low Voltage Grid separated and we compare 
each type of Grid in a separate subsection. 



6.1 Comparison of the Evolution Strategies for the Medium Voltage Distribution 
Grid 

It is difficult and perhaps even wrong to establish "the winner" between the various evolution strategies proposed. 
Also because the improvement in the various phases of evolution are not linear, a strategy that scores best in 
one evolution step could result worse compared to another when more edges are added. 

As a general remark we reinforce that the addition of links is beneficial to the topological metrics analyzed 
and one might speculate that the optimal solution in order to minimize (or maximize) the metrics is to realize a 
completely connected network (i.e., clique). With a real infrastructure such as the Power Grid this is impossible 
to achieve from the economical point of view, but also from a technical point of view (e.g., substations receiving 
thousand of connections). 

The comparison of the averages over all samples in each evolution step for the various evolution strategies is 
shown in Table[27] Take the triangle closure strategy, for instance, it is the one that in every step of the evolution 
process leads to the maximal values of the clustering coefficient, but it is also the one that has the worst results 
together with the least distance strategy for the metrics related to path length (characteristic path length and 10*'' 
redundant average path length). In addition, considering the robustness metric, the triangle closure offers the 
smallest improvements. Also considering the betweenness-related metrics, the triangle closure strategy results the 
worst compared to all the others with the second set of highest values for average betweenness and an increasing 
coefficient of variation. The same considerations apply to the least distance strategy that to a certain extent is 
very similar to the triangle closure. This strategy is the worst concerning the path length, having a characteristic 
path length more than double than the best evolution strategy (i.e., random); the same also applies for the 10*'' 
redundant average path length. Concerning redundancy aspects at the end of the evolution process, the least 
distance strategy is the one (except for random) that have the best robustness results. In terms of metrics related 
to betweenness, the least distance strategy shows results that are even worse compared to the triangle closure one. 
Considering the high degree assortative evolution strategy the results show that, overall, it is the best evolution 
strategy in the first two steps of the evolution; in the first step it is even better than the random one (with the 
exception of robustness). Considering betweenness metrics, this strategy is second only to the random strategy, 
with the exception of coefficient of variation which is worse than all other metrics. In general, this is a good 
strategy of evolving the network especially concerning the path-related aspects, with just robustness that lacks 
compared to other evolution strategies. The assortative low degree strategy is not particularly appealing and it 
is outrun by the assortative high degree strategy in every metric except robustness and coefhcient of variation of 
betweenness, therefore we do not find it particularly interesting for evolution purposes of the Medium Voltage 
networks. The dissortative strategy does not excel particularly in any of the metrics considered, however its 
values are quite fair, especially in the initial stage of the evolution. In particular, since the values of the metric 
do not improve particularly in the following steps of the evolution process, such evolution strategy could be used 
as a slightly more robust alternative compared to the assortative high degree in those scenario where the number 
of edges to be added is minimal and there is no need of special excellence in one topological parameter. When 
enough connectivity is added, the evolution strategy that scores best compared to all the others for all the metrics, 
with the exception of the clustering coefficient metric, is the random one. It has already been noted by Casals 
et al. [65j that some randomness in the network is beneficial especially for aspects related to robustness. In our 
analysis we have the same general results, with characteristic path length about 4 and a 10*'' redundant average 
path length of just 6. In addition, the networks evolved according to this strategy are the most robust with a 
value higher than 0.7 when the maximal connectivity is reached. The same considerations apply for betweenness 
metrics that obtain the best results when the connectivity is enhanced reaching a betweenness to order ratio 
smaller than three and a coefficient of variation below one at the last evolution step. However, it is difficult to 
propose for a Distribution provider to improve his Grid in a random fashion, even if part of the weaknesses and 
inefficiencies come from the lack of such randomness. Considering rational and evolution strategies that come 
with a motivation we consider the strategy assortative high degree as the one that, by evaluating the various 
topological metrics, scores best in the evolution tests that we have performed on the Dutch Medium Voltage 
Grids. 
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Table 27: Comparison of evolution strategies for Medium Voltage network. 



6.2 Comparison of the Evolution Strategies for the Low Voltage Distribution Grid 

As for the Medium Voltage evolution strategies, one cannot declare a winner. Also in this case some improvements 
between the various phases of evolution are not linear, therefore a strategy that scores best in one evolution step 
could result worse compared to another one when more edges are added. We point out once again that the 
addition of links is beneficial to the topological metrics analyzed and one might speculate that the optimal 
solution to minimize (or maximize) the metrics is to realize a completely connected network (i.e., clique). With a 
real infrastructure such as the Power Grid this is impossible to achieve due to economical and techncial obvious 
considerations. 

The comparison of the averages over all samples in each evolution step for the various evolution strategies is 



shown in Table 28 Once again, the strategy that scores best in comparison with the others is the random one for 
robustness and for metrics related to betweenness. For robustness in the last stage of evolution the value reaches 
on average 0.7. For betweenness-related metrics the two final stages of evolution are the best ones for average 
betweenness and average betweenness to order ratio. This strategy is the only one for which the coefficient of 
variation results below one (on average for all the samples) in the final stage of evolution. However, the random 
strategy has its weak point in the clustering coefficient metric that score below all the others evolution strategies. 
Considering the assortative strategies, one sees that the two strategies have quite similar scores for the metrics. 
In the very first stage, the high degree assortative strategy is better, but later in the evolution the assortative 
low degree slightly outperforms the assortative high degree one for clustering coefficient and characteristic path 
length, while for the robustness the difference is limited. The results between these two metrics for betweenness 
are quite similar and here the only interesting difference is in the coefficient of variation which scores best for 
the high node degree strategy. These two strategies are comparable to the random one for characteristic path 
length matters, while scoring worse for robustness, but better for clustering. Considering the clustering coefficient 
metric alone, the strategy that outperforms the others is the triangle closure strategy with values that are double 
compared to the others at the end of the evolution process. However, this strategy is worse than all the others 
for the characteristic path length. Although it has the highest values for average betweenness, the coefficient of 
variation ranks second compared to the other strategies considered. The strategy that slightly outperforms the 
others (except the values for robustness of the random one) is the dissortativc strategy. In fact, the path length 
is just smaller than 2.5 even better than the random edge addition. The clustering coefficient is in line with 
the assortative strategies and scores around 0.5. Concerning betweenness, this strategy scores best considering 
the average betweenness, but for the coefficient of variation this strategy scores worst. In the Low Voltage 
Distribution Grid it is even more difficult to propose a strategy to follow for the evolution given the similarity for 
the values of the evolution strategies considered. Considering non-random evolution strategies we consider the 
dissortative strategy that scores best in the evolution tests that we have performed on the Dutch Low Voltage 
Grids. 



49 





Evolution step 


CPL 


CC 


rob 


10'" red. path 


avg. bet. 


bet./ order 


coeff. var. 


Assortative high degree 


-\-^D /O 


o.DUU 


n (1(^9/1 Q 


n zinc 


( .000 


907 Qf^O 

zu ( .yoy 


9 f^l 
Z.OlO 


1 c;i fl 
l.Olo 


+ OU70 


3.225 


U.io7U2 


0.432 


5.797 


191.974 


2.227 


1.564 


+ / o/o 


z.y / 


n OQ/1 00 
U.zo4zZ 


U.40 / 


o.uoy 


ioz.yo4 




1.0 (O 


-|-1UU /c 


TOO 




n /I Q T 

U.4ol 


A 077 
4.y i ( 


T 7/1 c^OQ 


1 QQO 

l.ooU 


1 779 
1. M Z 


Assortative low degree 


-\-2b/o 


3.866 


0.05714 


0.405 


8.865 


zzbALb 


2.869 


1.436 


+ OU70 


3.074 


U.iZD (O 


0.446 


6.310 


-100 00 /I 


2.093 


1.746 


+75 7o 


2.827 


n 000 1 ^? 


0.509 


5.187 


177.925 


1.825 


2.113 


+ iUU/o 


Z. (DO 




U.OOl 


4. / oU 


1 i o.oiz 


1. ( OO 


i.yvSy 


Triangle closure 


-(-^0 /o 


o.yuo 


u. ioU 1 1 


U.04D 


Q 098 


/I QQ 1 97 
4oc5. IZ / 


4.oyy 


1 909 

1 .zuz 


+0U/0 


5.463 


U.z4yDi 


0.383 


8.070 


4Di.oyo 


2.044 


1.247 


-\- 1 0/0 


O.UzU 


U.o04o / 


U.400 


07Q 

D.y i 


/I 1 11/1 

4iy . 114 


1 Qf^O 

l.ooU 


1 Q1 


H-iUU /o 




U. 40104 


u.4y i 


D.O i 


•309 7Q1 
Oi)Z. i ol 


1 7Q9 
1. 1 OZ 


1 /1 9fi 
1 .4Zo 


Dissortative 


1 or:(V 


3.396 


n noQ /I Q 
U.Uyv54o 


0.403 


7.923 


194.500 


2.413 


1.849 


-\-o\J /o 






U.4oO 


0. IDo 


1 fl7 '^'^Q 

10 1 .ooy 


1.0/0 


9 1 s 
Z.olo 


( /o 


Z.U-LU 


n 91 "^09 




4 924 


itju. y i^L 


fifi7 




+100% 


2.476 


0.23151 


0.501 


4.618 


140.275 


0.619 


2.502 


Random 


+25% 


3.748 


0.03852 


0.440 


8.404 


219.250 


2.650 


1.159 


+50% 


3.171 


0.09465 


0.583 


6.323 


172.730 


2.113 


1.054 


+75% 


2.841 


0.14127 


0.636 


5.450 


149.836 


0.642 


1.038 


+100% 


2.621 


0.15564 


0.689 


4.831 


132.417 


0.567 


0.978 



Table 28: Comparison of evolution strategies for Low Voltage network. 



6.3 Desiderata parameters satisfaction for Medium Voltage networks 

We now consider the satisfaction of the quantitative metrics that we have proposed in our previous work |59j to 
analyze the appropriateness of synthetic topologies to improve the Grid in reducing losses, facilitate local energy 
distribution, and increasing network robustness. The metrics that we require to be satisfied in the future Smart 
Grid can be can be categorized into three macro categories with respect to how they affect a Power Grid: efficiency 
in the transfer of energy, resilience in providing alternative path if part of the network is compromised/congested, 
and robustness to failures for network connectivity. The satisfaction of these metrics, provides topologies where 
the cost of electricity should be reduced since the considered metrics directly influence the cost of distributing 



electricity j58j. For a through insight of the metrics, we refer to the Appendix A or to |58j where a detailed 



explanation is provided. The results are shown graphically in Tables [29] and 30 where for each evolution strategy 
and step of evolution the satisfaction of the parameter is marked with a "tick" sign (/), while the dissatisfaction 
is marked with a cross sign (X), and a value closer to satisfaction by an approximation sign (w). 



Assortative high degree The assortative high degree evolution strategy satisfies the desired values for the 
metrics concerning the characteristic path length already from the second step of evolution for ten out of twelve 
samples. All samples satisfy the requirements over the redundant path. Considering the clustering coefficient, 
the majority of samples satisfies the requirement at the third step of evolution where seven samples have a 
CC > 5 X CCrg', in the final stage of evolution this metric is satisfied by all the samples except two. The real 
drawback of this metric, as shown in the table, is represented by the robustness that never reaches the goal of 
0.45, but stops around 0.37. The unsatisfaction of the metric is also present for betweenness ones: all samples 
are one unit larger even in the last step of evolution. In addition, the coefficient of variation never reaches the 
target for the samples at any stage of evolution. 



Assortative low degree The assortative low degree evolution strategy satisfies the desired values for the 
metrics concerning the clustering coefficient already from the second step of evolution with eleven out of twelve 
samples. All the samples satisfy the requirements over the redundant path already from the first step. Considering 
the clustering coefficient, the majority of samples satisfies the requirement at the second step of evolution where 
seven samples have a CC > 5 x CCrq', in the final stage of evolution this metric is satisfied by all the samples 
except one. This metric almost provides the satisfaction of the robustness requirement from the third step on 
by having six samples fully satisfying it and three very close to the 0.45 threshold. The metrics concerning 
betweenness are not satisfied both for betweenness to order ratio and for the coefficient of variation which never 
reach the target for the samples at any stage of evolution. 



Triangle closure The triangle closure evolution strategy focuses on the clustering coefficient, therefore reaching 
the target already in the first step of the evolution. The evolved graphs do not only satisfy the requirement posed 
of having CC > 5 x CCrg, but also a more restrictive requirement of CC > 10 x CCrg which can be considered 
the condition for satisfying the small- world requirement for this property cf. |83| . Concerning the path length 
properties, the characteristic path length is never lower than the logarithm of the order of the graph. The 
requirement over the 10*'' average path length is satisfied already from the first addition of edges. On all other 
requirements this strategy is weak. For robustness just few samples reach values around 0.4, while the majority 
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is about 0.3 when the most of the edges are added. For metrics that involve betweenness, this evolution strategy 
scores poorly and never even close for any single sample to the target for both average betweenness to order ratio 
and for the coefhcient of variation. 



Dissortative node degree The dissortative node degree strategy scores quite poorly. As mentioned before, 
this strategy does not excel in one specific metric, but all the metrics are slightly improved. Such aspects result 
in a limited crossing of the threshold for the desiderata parameters imposed. Only after the third evolution step 
the characteristic path length is almost satisfied: nine of the twelve samples reach the target, while one more 
reaches the target when even more edges are added. A similar condition is true for the robustness metric. In fact, 
six samples satisfy the target after the third step and other three are above the 0.4 value. Clustering coefficient is 
always well below the threshold, having values that often are just slightly higher than the clustering coefficient of 
a random graph with same order and size. The only metric fully satisfied is the one related to redundant paths 
in the network. While the metrics related to betweenness are not satisfied even at later stages of evolution. 



Least distance The results of the least distance strategy are similar to those of triangle closure. The clustering 
coefficient desiderata is reached already in the first step of the evolution. From the second step on, the evolved 
graphs do not only satisfy the requirement posed of having CC > 5 x CCrG: but also a more restrictive 
requirement of CC > 10 x CChg which can be considered the condition for satisfying the small- world requirement 
for this property cf. [53]. Concerning the path length properties, the characteristic path length is never lower 
than the logarithm of the order of the graph. However, the easy requirement related to the lO"* average path 
length is satisfied already from the first addition of edges. Robustness is another positive note of this evolution 
strategy. From the third step on the seven samples reach the 0.45 target and other three are close to 0.4; in 
the last step basically all the Medium Voltage evolved samples satisfy the desiderata. The situation is not so 
positive for betweenness related metrics since they score worst compared to the others strategies for both average 
betweenness to order ratio and for the coefficient of variation, therefore having two crosses for these cells in the 
table. 



Random The random evolution strategy is the one that satisfies most of the desiderata parameters not only 
after the first step of evolution, but already after the second one. The metrics concerning the characteristic path 
length are satisfied basically already from the second step of evolution. The same applies to the redundant path 
whose goal is met already in the first evolution step. Robustness is really impressively achieved after the addition 
of 50% of more edges, with all the samples above the 0.45 threshold and the majority of them scoring even higher, 
above 0.6. The metrics related to betweenness that usually fail for the other samples are here met completely on 
the coefficient of variation side after the third evolution stage (eight out the twelve samples meet the target). In 
addition, this evolution strategy is the one that also goes closer to the satisfaction of the betweenness to order 
ratio. The only problem of such a strategy is on the clustering coefficient side. With such a strategy the formation 
or closure of topological triangle structures is really difficult, therefore the cross sign for this parameter. 



As a general remark about the satisfaction of the desiderata parameters is what we have in part already 
noted in the comparison of the result for topological quantities. The random evolution strategy is the one that 
satisfies most of the desiderata requirements (three parameters fully satisfied and one almost satisfied) already 
at the second stage of evolution and from the third step on, four requirements are fully satisfied. Concerning the 
other strategies that have an evolution with a specific goal, those that satisfy most of the requirements are the 
assortative high degree and the least distance one. Both strategies satisfy three out of the six parameters and 
the difference being that the assortative one has good performance for the characteristic path length parameters, 
while the least distance strategy reaches the target for the robustness aspects. Therefore, these sub optimal 
strategies could be used where such different requirements are most needed. 

From the comparison of the different evolution strategies and their steps in adding new edges as shown in 
Tables 29 and |30j there are evolution strategies that tend to satisfy the majority of the metrics for the Smart 
Grid that we have defined (cf. Appendix A or [59j). The evolution strategy that already in the second step 
fully satisfies three requirements is the random evolution. The requirements satisfied become four when 75% of 
additional edges are added. This strategy is in line with the finding of our previous work [SH] where a small-world 
model was the best solution for modeling a Smart Grid topology. The addition of random edges goes into that 
direction: the network with the rational structure planned by the power engineers is modified by the addition 
of random links, thus in something between the rational topology and the fully random network. The other 
strategies that are quite successful are the assortative high degree and the least distance but only when the 
edges are doubled. This is the situation considering a pure topological decision. When also the costs of the 
evolution are taken into the picture the optimal evolution strategy might change due to economical constraints. 
The extended analysis that consider economic aspects is performed in Section [Tj 
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6.4 Desiderata Parameters Satisfaction for Low Voltage Networks 



We now consider the satisfaction of the quantitative metrics that we have apphed for the Medium Voltage 
networks for the Low Vohage evolved samples. The aspects that these metrics assess deal with losses, local 
energy distribution, and increasing network robustness. These are the aspects that we require to be satisfied in 
the future Smart Grid. The results are shown graphically in Tables [31] and [32] where for each evolution strategy 
and step the satisfaction of the parameter is marked with a "tick" sign (/), while the dissatisfaction is marked 
with a cross sign (X), and a value closer to satisfaction by an approximation sign («). In Tables 31 and 32 we 



consider the satisfaction of the clustering coefficient metric, therefore assigning a "tick" sign, if the clustering 
coefficient desiderata requirement is satisfied by the three biggest samples in the network (i.e., sample #5, #9, 
and #10). This is done in consideration of the extreme high values of clustering coefficient that the desiderata 
parameter impose when samples have a small order. 



Assortative high degree The assortative high degree evolution strategy almost satisfies the desired values for 
the metrics concerning the characteristic path length from the second step of evolution with eight out of eleven 
samples; while fully satisfaction is in the last evolution step. However, the requirement for the redundant path 
are never fully satisfied. Considering the biggest samples, they never satisfy the condition CC > 5 x CCro- Also 
robustness requirement are not satisfied, only six samples in the last two evolution stages reaches the goal of 
0.45. The betweenness requirement is only partially satisfied: eight out of the eleven samples have a betwenness 
to order ratio around the target of 2.5. Considering the coefficient of variation, only the two smallest samples 
satisfy this property. 



Assortative low degree The assortative low degree evolution strategy satisfies or is close to satisfying the 
metrics in the last stage of the evolution process. For characteristic path length, already from the second step, 
ten of the eleven samples reach the target; on the contrary, the target is never reached for the redundant path 
requirement. The three biggest samples score sufficiently concerning the clustering coefficient requirement, while 
the others have high values, but never satisfy the high- demanding requirements. Robustness is close to satisfaction 
in the last stage of evolution where all samples reach the target except two that stop to a value around 0.4. The 
metrics concerning betweenness show an almost satisfaction of the average betweenness to order ratio with eight 
of the samples well below the target; however the coefficient of variation never reaches the target for the samples 
at any stage of evolution. 



Triangle closure The triangle closure evolution strategy focuses on the clustering coefficient, therefore reaching 
the target already after the second step of the evolution. Two of the three biggest samples satisfy also a more 
strict condition than the desiderata requirement such as CC > 10 x CCrqi that can be considered the condition 
for satisfying the small- world requirment for this property cf. |83j. Concerning the path length properties, the 
characteristic path length is never lower than the logarithm of the order of the graph. However, the requirement 
over the 10*'' average path length is satisfied already after the second step of evolution. In fact, the redundant 
paths are smaller than twice the characteristic path length given the quite high characteristic path length that the 
networks have. For robustness, seven out of the eleven samples reach the target in the last step of evolution. For 
metrics that involve betweenness, only seven samples are below the target for the betweenness to order metric, 
actually the smallest in order, while the three biggest samples are far from the target. The coefficient of variation 
increases in the evolution, therefore not allowing for this metric to gain a "tick" sign. 



Dissortative node degree The dissortative node degree strategy scores quite well for evolving Low Voltage 
networks. Especially in the last stage of evolution half of the metrics are satisfied or very close to satisfaction. 
Already in the second step of evolution, the characteristic path length is satisfied; actually already in the first 
step of evolution eight samples satisfy the condition. Robustness is satisfied with nine out of the eleven samples 
that fully comply with the desiderata and the other two that are above 0.4. The real weak point is the clustering 
coefficient which even for the biggest samples is never close to the target. Betweenness-related metrics have a 
dissimilar behavior: the betweenness to order ratio satisfies the desiderata for nine samples, while the coefficient 
of variation has an increasing trend and in the final step does not comply with the requirements for any sample. 



Random The random evolution strategy proves to be one of the strategies that satisfies most of the parameters. 
The metrics concerning the characteristic path length are satisfied basically already from the second step of 
evolution with ten of the eleven samples compliant. A different behavior applies to the redundant paths whose 
goal is not met, but in the final step eight of the eleven samples reach the target. Robustness is almost achieved 
already after the second step: all the samples score above 0.4 and just three do not reach the target in that step, 
but later in the evolution. In the last step, the satisfaction of the redundant path is almost satisfied with eight of 
the samples (containing also the three biggest ones). The metrics related to betweenness that usually fail for the 
other types of evolution are here almost entirely met concerning the coefficient of variation and the betweenness 
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to order ratio in the last stage of evolution. The main drawback is on the clustering coefficient side. With such 
a strategy the formation or closure of topological triangle structures is really difficult, therefore the cross sign for 
this parameter. 

The random evolution strategy is the one that satisfies most of the desiderata (two parameters fully satisfied 
and three almost satisfied) requirements at the final stage of evolution. Concerning the other strategies that 
have an evolution with a specific goal those that satisfy most of the requirements are the assortative low degree 
followed by the dissortative one which almost satisfy four and three requirements respectively. Therefore these 
sub optimal strategies could be used too. In particular, the assortative low degree satisfies the requirements 
regarding the clustering coefficient, that are not satisfied neither by the random, nor by the dissortative strategy. 



A commonality between the evolution for Medium Voltage and Low Voltage is the best result achieved by 
the random evolution of network that scores best among the strategies considered in this study. An interesting 
difference lies in the sub-optimal strategies that score high in for the Medium Voltage network and Low Voltage 
network evolution. For the Medium Voltage network, it is best to provide more connectivity between those nodes 
that already hold a high node degree, thus reinforcing their role as key components of the network. On the 
other hand, for the Low Voltage network, it seems that adding more connectivity between the nodes that have a 
small connectivity provides better performances for evolution. The assortative low degree aims at creating more 
connections and (therefore hubs) where is less connectivity in the current samples. Also the other strategy that 
is sub-optimal for Low Voltage networks, the dissortative, aims at giving a more important role to the nodes in 
the periphery of the network by connecting them to the more connected nodes. A synthesis of the topological 
performance of the different strategies is provided in Table 33 For the layer of the Power Grid considered (column 
one), each strategy (column two) is assessed with an optimality level (column three) based on the satisfaction 
of the topological metrics. The step during the evolution process in which the most of the metrics are satisfied 
is also provided in column four. The achievement of the good results depends on the network layer and on the 
strategy used. The optimality level is assigned based on the full satisfaction of the topological metrics described 
earlier in this section and more thoroughly in [Appendix A[ One can see a distinction between the Medium Voltage 
and the Low Voltage networks: the Low Voltage achieve less topological optimality (less metrics are satisfied) 
although requiring on average more evolution steps. Another aspect to note is that some evolution strategies (i.e., 
dissortative and triangle closure) do not have an improvement between the various steps despite the addition of 
more connectivity. On the other hand, other strategies (i.e., assorative and random) benefit more from additional 
connectivity by improving the satisfaction of the metrics. 
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Network 
level 


Strategy 


Optimality level (0 to 6 *) 


Step for optimality 


MV 


Assortative HD 


★ ★ ★ 


1 


MV 


Assortative LD 


*★ 


2 


MV 


Triangle closure 


★it 


1 


MV 


Dissortative 


it 


1 


MV 


Random 


it it *★ 


3 


LV 


Assortative HD 


it 


4 


LV 


Assortative LD 


it 


4 


LV 


Triangle closure 


itit 


2 


LV 


Dissortative 


it 


2 


LV 


Random 


itit 


4 



Table 33: Parameter satisfaction and evolution step of achievement. 



7 Economic Considerations 

Traditionally the problem of evaluating the expansion of an electrical system is a complex task that involves 
both the use of modeling, usually based on operation research optimization techniques and linear program- 
ming SB) IS] , and the experience and vision of experts in the field aided by computer technologies [37] . In this 
latter case, computers acquire knowledge based on previous expert decision and, based on the electrical physical 
constraints of the domain, are then able to support Power Grid evolution decision finding the most suitable 
technical and economical solution [73]. With more distributed generating facilities at local scale, traditional 
methods have limits and need to be modified or updated to take into account the new scenarios of the Smart 
Grid. The strategies that we have so far analyzed as being candidates for evolving the current Grid in the future 
Smart Grid need also to be evaluated from the economic point of view. How much will it cost to deploy electrical 
infrastructures according to these models? What is the actual cost of adding an edge to the topology of the Grid, 
physically? 

To answer these questions, we need to characterize the cables used for the upgrading process in terms of 
their physical properties and costs. The information of length for each cable in the existing sample networks 
is available. For the majority of the nodes in the samples belonging to the Medium Voltage we know the 
geographical coordinate information, so with some approximation it is possible to compute the length of the new 
cable connecting nodes that do not have a line connecting them yet. Therefore, each new added line is assigned 
a length. In order to define the properties characterizing the new cables, we adopt a k-nearest neighbor (KNN) 
classification based on the length of the cable [221 US] ■ The classification is performed using the Java Machine 
Learning Library (Java-ML) vO.1.5 (http : // j ava-ml . sourcef org e . net/[ ) . For each new cable, we consider the 
k cables (we choose /c = 5) in the original network used as the training set for the classifier that have the length 
closer to the added one and then the type of cable is available. We choose this method to identify the cable 
since it is simple and it has proven successful in several practical contexts |68| . Once the properties of the new 
cables are identified, all the information we need for the economic analysis are then available: resistance per unit 
of length, cost per unit of length, maximal supported current. We remark that our proposed analysis does not 
aim at being a comprehensive investment analysis, for which other techniques are well established [72], simply 
an economic evaluation of the proposed evolution strategies to confirm their feasibility or unfeasibility. 

Once again, as in our previous work [59], to assess these costs in the Medium Voltage infrastructure, we 
consider a simple relation where the cost of cabling and cost of substations are added: 

N M 

Cirnpl ^^SsCj + ^CCi (1) 

3=1 i=l 

where C'impi stands for implementation cost, Sscj is the adaptation cost for the substation j and Cci is the cost 
for the cable i. The cost of the cable can be expressed as a linear function of the distance the cable i covers: 
Cci = Cuci ■ h, where Cua is the cable cost per unit of length and k is the length of the cable. There are several 
types of cables used for power transmission and distribution with varying physical characteristics and costs. In 
addition, the cost for installation can vary significantly |55j . In the present work, to provide an initial estimate, 
we simply consider cabling costs and ignore substation ones. While the former are directly tied to the topology 
and length of the links, the latter pricing is too dependent on other factors to be useful in the present analysis 
(e.g., different equipment in the substation). 

For the cost analysis, we limit our investigation to the three strategies that have scored best in the pure 
topological analysis (i.e., random, assortative high node degree and least distance) described in Section |4j In 
addition, the economic analysis can be applied only to the Medium Voltage network samples since there is no 
geographical information about the location of the nodes. The results of the cost for evolution of the Dutch 
samples are shown in Tables [34] [35] and [36] The first column of each table contains the sample ID, while 
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the second provides the information about the evolution step considered, the third column has the information 
concerning the cost of the evolution of the network according to the specified strategy. The fourth column contains 
the information on the fraction of the cost that the evolution of the infrastructure impacts on the whole cost of 
the infrastructure. From the tables the most interesting result is the difference in the cost that the three methods 
of evolution require. On average the costs of network improvement is similar when considering the assortative 
and the random strategy: the cost of adding more edges is about 75% of the cost for cabling of the whole Grid 
infrastructure. The situation is radically different for the least distance strategy whose development impacts only 
marginally in the total cost of the infrastructure, by adding just less than 13% of the cost of the infrastructure. 
Purely from the point of view of cabling costs, the most promising evolution strategy in economic and topological 
terms for the Medium Voltage Grid is the strategy that connects the nodes that are geographically closer (i.e., 
least distance strategy). 





Evolution step 


Evolution cost (euro) 


Fraction of evolution 

cost on the whole infrastructure 


Sample #1 


+zo/o 




U. / o 


+50% 


59334032 


0.73 


+75% 


59334032 


0.73 


+ 100% 


119505699 


0.85 


Sample #2 


+25% 


107698129 


0.73 


+50/0 




0.85 


+75% 


255301286 


0.87 


+ 100% 


298792505 


0.88 


Sample #3 


+2570 


o o r\r\ T cr a a 

23007544 


0.72 


+50% 


44739934 


0.83 


+75% 


72515143 


0.89 


+100% 


105567690 


0.92 


Sample #4 


+25% 


45762304 


0.65 


+50/0 


r r; 1 KCiQCi'v 

ooloyoy / 


0.69 


+75% 


55159897 


0.69 


+100% 


55159897 


0.69 


Sample #5 


+2o/o 


d5v59oDo1 


0.74 


+50 /o 


11 / 984090 


0.84 


+ /5/0 


15 /8/D214 


0.87 


+ 100% 


194838593 


0.90 


Sample #6 


+25% 


31063470 


0.67 


+50 /o 


629/ 1do8 


0.81 


+ /5/0 


104 ( 11d8 / 


0.87 


+100% 


165938342 


0.92 


Sample #7 


+25 /o 


'70 1 /I o 1 o 
r91481o2 


0.61 


+50% 


199600692 


0.80 


+75% 


229365781 


0.82 


-t-iUU/o 


oou t oZyoo 


U.o / 


Sample #8 


+25% 


32826377 


0.67 


+50% 


36976365 


0.70 


+75% 


83305090 


0.84 


+ 100% 


123241632 


0.88 


Sample #9 


+25% 


11739091 


0.46 


+50% 


42894090 


0.76 


+75% 


72834999 


0.84 


+ 100% 


103147797 


0.88 


Sample #10 


+25% 


14899834 


0.60 


+50% 


32136672 


0.77 


+75% 


45646927 


0.82 


+100% 


60546613 


0.86 


Sample #11 


+25% 


13094364 


0.50 


+50% 


26377193 


0.67 


+75% 


38709007 


0.75 


+ 100% 


40774034 


0.76 


Sample #12 


+25% 


93369553 


0.70 


+50% 


93369553 


0.70 


+75% 


138091611 


0.77 


+ 100% 


253056714 


0.86 



Table 34: Cost of Medium Voltage network evolution for assortative high degree strategy. 

One may then wonder if such investments are beneficial for the end users and the distribution companies in 
reducing the cost of electricity flows. We resort to a set of metrics that we have developed and already applied 
to the Northern Netherlands Distribution Grid [58] and to synthetically generated networks [59]. The goal is to 
consider, from a topological perspective, those measures that are critical in contributing to the cost of electricity 
as elements in the Transmission and Distribution Networks as described in economic studies such as the one of 
Harris and Munasinghe [211 [H] : 
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PiVollli'lOTl ^1"PT> 
-1— J V d vlWlX O 
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J_> V \J± lAvLKJLl y^KJ^K/ \^UIXV^J 


Praption of pvnliTfimi 

cost on the whole infrastructure 


Sample #1 


+25% 


49585634 


0.70 


+ 5U/C 


99550100 


0.82 


+7570 




0.87 


+ iUU/o 


iyoUoUolo 


n on 
U.yU 


Sample #2 


+25% 


81443232 


0.67 


+5U7C 


173835139 


0.81 


+75% 


270738228 


0.87 


+ 100% 


352501220 


0.90 


Sample #3 


+25% 


20401026 


0.70 


+50% 


39401799 


0.82 


+75%. 


58698139 


0.87 


+ 101) A, 




n on 


Sample #4 


+25% 


38550574 


0.61 


+5070 


74405909 


0.75 


+75% 


112507292 


0.82 


+ 100% 


144092688 


0.85 


Sample #5 


+257C 


34665986 


0.60 


+50% 


63522559 


0.74 


+7570 


n n "1 T "1 n ,i 

9U1719o4 


0.80 


+ iOO/o 


i iO i 4Do10 


U.o4 


Sample #6 


+25% 


32756917 


0.68 


+507C 


D /90DOOy 


0.82 


+75% 


102189138 


0.87 


+ 100% 


1 /in/^onooo 

14UDo2z3o 


0.90 


Sample #7 


+25% 


85124154 


0.63 


+5070 


1 oi nncTTcn 

181095750 


0.78 


+75% 


262020773 


0.84 


+100% 


352245673 


0.87 


Sample #8 


+25% 


40434688 


0.72 


+50% 


79435946 


0.83 


+ /O/o 


iU( f54oDOU 


\J.O t 


+ 100% 


139828749 


0.90 


Sample #9 


+25% 


31096020 


0.70 


+50% 


57036433 


0.81 


+75% 


85190089 


0.86 


+ 100% 


115442294 


0.89 


Sample #10 


+25% 


12677249 


0.56 


+50% 


26126828 


0.73 


+75% 


42325971 


0.81 


+ 100% 


53357492 


0.84 


Sample #11 


+25% 


1214608 


0.08 


+50% 


2553383 


0.16 


+75% 


3117824 


0.19 


+100% 


4494600 


0.25 


Sample #12 


+25% 


105240341 


0.72 


+50% 


223059810 


0.85 


+75% 


317337397 


0.89 


+ 100% 


417627345 


0.91 



Table 35: Cost of Medium Voltage network evolution for random strategy. 
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cost on the whole infrastructure 


Sample #1 


+25% 


635014 


0.03 


+ 5U/C 


1944cJ5o 


0.08 


+75% 


3536060 


0.14 


+ iUU/o 


oooyU / i 


n on 
U.zU 


Sample #2 


+25% 


1832132 


0.04 


+5U7C 


5362070 


0.12 


+75% 


9454081 


0.19 


+ 100% 


15488444 


0.28 


Sample #3 


+25% 


357093 


0.04 


+50% 


1143678 


0.11 


+75%. 


2075073 


0.19 


+ 100'X, 


320(>l()t} 


(J. 26 


Sample #4 


+25% 


831466 


0.03 


+5070 


2936596 


0.11 


+75% 


5026907 


0.17 


+ 100% 


7688212 


0.24 


Sample #5 


+25% 


801764 


0.03 


+50% 


2128363 


0.09 


+75% 


3412290 


0.13 


+ iOO/o 


oio / zoi 


u.iy 


Sample #6 


+25% 


543105 


0.03 


+507C 


1433992 


0.09 


+75% 


2588723 


0.15 


+ 100% 


662262':) 


0.20 


Sample #7 


+25% 


1104605 


0.02 


+5070 


3521507 


0.07 


+75% 


6534753 


0.11 


+100% 


10042275 


0.17 


Sample #8 


+25% 


584039 


0.03 


+50% 


2127943 


0.12 


+ /o/o 


4Ui0oy4 


n on 
U.zU 


+ 100% 


5840352 


0.27 


Sample #9 


+25% 


661698 


0.05 


+50% 


1861427 


0.12 


+75% 


3476464 


0.20 


+ 100% 


5524643 


0.29 


Sample #10 


+25% 


314451 


0.03 


+50% 


987993 


0.09 


+75% 


1971402 


0.17 


+ 100% 


2912718 


0.23 


Sample #11 


+25% 


478534 


0.03 


+50% 


1230921 


0.09 


+75% 


2224697 


0.14 


+100% 


3262358 


0.20 


Sample #12 


+25% 


1069466 


0.03 


+50% 


3143898 


0.07 


+75% 


5723602 


0.12 


+ 100% 


8686081 


0.18 



Table 36: Cost of Medium Voltage network evolution for least distance strategy. 
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• losses both in line and at transformer stations, 

• security and capacity factors, 

• line redundancy, and 

• power transfer limits. 

For each of these elements we associate topological quantities that are representative, so that we can assess their 
values and therefore compare them in the original networks and in the evolutions that we propose in this paper. 
Here we assess the 

• Losses on the transmission/distribution line can be expressed by the quotient of the weighted characteristic 
path length : 

WCPLn (2) 



The average resistance of a line (a weighted edge in the graph): 



1 

i=l 



• Losses at substation level are expressed as the number of nodes (on average) that are traversed when 
computing the weighted shortest path between all the nodes in the network: 

-^substations 

= NodeswcPLM (4) 



• Robustness is evaluated with random removal strategy and the weighted-node-degree-based removal by 
computing the average of the order of maximal connected component between the two situations when the 
20% of the nodes of the original graph are removed. It can be written as: 

\MCCjiandom20% \ + \MCCNodeDegree20%\ /r\ 

RoOn ^ (5) 

• Redundancy is evaluated by covering a random sample of the nodes in the network (40% of the nodes 
whose half represents source nodes and the other half represents destination nodes) and computing for 
each source and destination pair the first ten shortest paths of increasing length. If there are less than ten 
paths available, the worst case path between the two nodes is considered. To have a measure of how these 
resilient paths have an increment in transportation cost, a normalization with the weighted characteristic 
path length is performed. We formalized it as: 

^^^^ = WCPL 



• Network capacity is considered as the value of the weighted characteristic path length, whose weights are 
the maximal operating current supported: 

CapN = WCPLcurrentN (7) 



The average supported current of a line (a weighted edge in the graph) : 



I M 

'^current ^ ^ '^{current (8) 

1=1 



The aspects considered are just some of the factors (the ones closely coupled to topology) that influence the overall 
price of electricity. Naturally, there are other factors that influence the final price, e.g., fuel prices, government 
strategies and taxation, etc., as illustrated for instance in the economic studies of Harris and Munasinghe [5^155] . 

Figures [16] [17] [18] and [19] show the improvement as a fraction of the original values in the metrics related 
to the cost of electricity distribution respectively in the four step of evolution we consider (-1-25%, -1-50%, -1-75%, 
and -1-100% of edges) compared to the original samples. To give a general idea of the improvement, we use the 
average results over the 12 Medium Voltage networks samples. For each strategy, we look at the improvement, in 
percentage, for the same metrics compared to the initial samples. Already from a small increase in connectivity 
(addition of 25% of edges) the resistance of a path decreases with every strategy with the best result obtained by 
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the assortative high degree strategy. Adding more lines in general promotes a reduction of the average resistance 
of edges compared to the initial situation. This is even more true in the assortative and least distance strategy 
where short distance connections (and therefore less resistive) tend to dominate or more efficient cables are 
chosen in the KNN procedure selection. In addition, a considerable reduction in the losses that are experienced 
in traversing substation can be avoided with more dense networks: almost a 30% less station traversed in the 
graph enhanced with assortative connections; of course less benefits take place with the least distance strategy 
that reduces the average number of traversed stations in a path by 7%. This higher connectivity provides 
benefits to the robustness aspects: resilient paths that are less lossy (about 30% reduction for the least distance 
strategy). An improvement is also experienced in the resilience of the network to node disruptions: about 80% 
more resilience with the random strategy. Also the amount of current that is supported in a path (the weighted 
Characteristic Path Length with maximal supported current as weight of the edges) is definitely higher (more 
than 2.5 times for the least distance strategy) compared to the initial samples. Similar considerations can be 
done for the second evolution step (i.e., +50% cables). Considering the loss aspect the benefits are for the three 
strategies on average about 30% in the reduction of the weighted path. For the average weight (i.e., resistance) 
of a cable the random addition creates networks with much more cables with an high resistance. This is an 
indication that the random strategy creates long distance connections that act as shortcuts in the network from 
a topological point of view. However, in a physical system there are no significant benefits since the resistance 
grows with the distance. These limits of the benefits are shown by the weighted path analysis. Considering 
robustness oriented metrics, the most of the benefits take place in resilience where the networks on average stay 
twice as much as connected than the initial evolution step, except for the assortative strategy that lags behind. 
We see a small decrease for all the strategies in the characteristic path when the weight is the maximal current 
that can flow in the cables; this is not difficult to understand since the characteristic path (cf. Definition [l0| ) 
uses the most efficient (i.e., with smallest weight paths). When even more connectivity (+75% and +100%) is 
added one sees even better results and improvements but smaller in their magnitude. In these last two stages, 
the benefits for the losses reach values of reduction above 30%, the same for the nodes to be traversed on average 
on a path. Only the least distance strategy has a small reduction, but this is due to the very nature of such 
evolution strategy that avoids long distance (topological shortcuts) cables. Considering network robustness, the 
main improvements are for resilience to node failures that for the least distance strategy allow the network to be 
twice as much as robust then the initial samples. Concerning the capacity of the network we assist to a slightly 
decrease for the same reasons that we have explained above. However, the networks with the least distance 
evolution are able on average to transport twice the amount of current than the initial samples. 

Given the cost analysis for cable addition and the benefits that such an higher connectivity brings to the 
networks considering the different strategies in the topological aspects related to price of electricity, we provide 
a suggestion on how to evolve the analyzed samples. Using our topological-based method, we consider that the 
least distance strategy with the addition of 75% or 100% of edges is a very good way of making the network more 
connected. Such an evolution strategy provides benefits from a topological point of view, and it keeps the costs 
low. For the samples examined it requires an investment in cabling costs about 25% of the value of the cables 
already on the ground. Such investment can provide benefits in the loss reduction of the network as well as in its 
robustness. The energy economic studies show [311 ISHj that losses and robustness are directly related to the cost 
of electricity. These factors are tight to topological parameters as shown above, therefore such evolution would 
provide less cost in electricity distribution. It is difficult to translate in exact monetary terms for the end user 
the savings in the energy bill that an improvement in the topology parameters that influence the cost of energy 
distribution might bring. For the realization of a Smart Grid less costs in the distribution of electricity and a 
local reliable and robust network are the essential ingredients to enable a paradigm where energy is produced 
and distributed locally such as neighborhood or city level. 
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(a) Losses-related metrics. (b) Reliability-related metrics. 

Figure 16: Comparison to the original networks of properties influencing electricity cost at 1** stage of evolution 
(i.e., 4-25% edges). 
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(a) Losses-related metrics. 
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(b) Reliability-related metrics. 



Figure 17: Comparison to the original networks of properties influencing electricity cost at 2""^ stage of evolution 
(i.e., +50% edges). 
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(b) Reliability-related metrics. 



Figure 18: Comparison to the original networks of properties influencing electricity cost at S'"'* stage of evolution 
(i.e., -1-75% edges). 
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Figure 19: Comparison to the original networks of properties influencing electricity cost at 4*'' stage of evolution 
(i.e., +100% edges). 



8 Related Work 

To the best of our knowledge there are no works that use the principles of Complex Network Analysis to design 
and consider evolutions of Distribution Grid networks. Thus, in considering related work, we look into three 
main related areas: Evolution of Complex Networks, Complex Network Analysis studies applied to the Power 
Grid, and electrical engineering power lines design and adaptation. 



Evolution of Complex Networks 

The problem of the evolution of networks by using Complex Network Analysis techniques has been widely 
investigated. The focus is mainly to assess and analyze existing networks that have evolved over time and 
understand which principle is underlying such evolution. More in general, the evolution of networks is considered 
in defining abstract algorithms that can be used to build networks that mimic the same behavior of natural or 
synthetic networks. For a general introduction to Complex Network Analysis generation techniques, we refer to 
the extensive works such as [121 EZ]- The work of Liben- Novell and Kleinberg [ITj discuss link prediction. The 
authors consider the evolution of the ArXiv library co-authorship network (i.e., a social network) and define several 
quantitative metrics to infer the future links to be formed in the network (i.e., new co-authored publications). 
The metrics that the authors use for predictions are based on evaluating the paths, the neighborhood of nodes, 
rank of the matrix representing the graph, and clustering. Compared to a random prediction of new links the 
methods used in the paper are far better: the random predictions is correct less than 1% of the times, while on 
some topics of the ArXiv library the predictions reach a correct outcome 50% of the time. The ideas used in 
the paper are interesting and valuable, and the metrics used in the prediction of new connections between nodes 
prove valuable for the co-authorship social network. Even if at abstract level network are the same (i.e., graphs) 
it is difficult to apply these metrics and prediction strategies in a real infrastructure such as the Power Grid. In 
fact, many of the metrics considered in the paper are inspired by the underlying social aspects of the co-author 
relationship such as neighbors in common. This aspect is reinforced by the different results that the same metrics 
used on different topics of the ArXiv library provide sometimes far different results in the predictions following 
the different metrics. Given the small average node degree of the Power Grid networks, the absence of samples 
at different time stamps, it is difficult to apply the metrics and methods of the paper to Distribution Networks 
and expect meaningful and interesting results. 

The survey of Lii et al. |49j provides a good overview on link prediction. The authors propose several methods 
and algorithms based on similarity properties of the nodes at local level (e.g., the neighborhood of a node), global 
similarity (i.e., over the whole network), and quasi-local that combines the best of the two approaches. They 
also propose other strategies such as probabilistic models to evaluate the likelihood of the formation of a link. 
One of the networks that is considered for link prediction in [35] is the Western U.S. Power Grid (data already 
used in |82j) and in this case the methods that give more stress on the local interactions between nodes perform 
better in the prediction than other approaches. The authors also state that for the infrastructural networks the 
geography plays an important role where the long distant connections are discouraged and very rare. The paper 
describes three main application of link prediction in the field of Complex Network Analysis: reconstruction of 
networks, evolution mechanism of networks and classification of partially labeled networks. We consider these 
prediction mechanisms valuable for networks with rich neighborhood structures where differences and intersection 
between neighborhoods are meaningful (e.g., social networks) or where the discovery of links is particularly costly 
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or long in time and efforts (e.g., protein interactions). The current Power Grid networks have low average node 
degree and the structure is usually close to radial, especially at Medium and Low Voltage level, thus we consider 
the predictions less interesting. However, these techniques and link predictions might be more applicable in the 
further evolution and prediction of new links when the networks are more rich in connectivity. In addition, some 
techniques for link prediction would benefit from additional information available to characterize the nodes, i.e., 
the number of prosumers that are connected to a node (i.e., substation) and the available production power. 

The concept of network evolution of a technological infrastructure is considered in where a model for 
the evolution of the Internet Autonomous System (AS) is provided. As most works (e.g., [28l [6]), this one 
focuses on creating a model of the existing evolution of the Internet rather than proposing new ways of evolving 
networks. The interesting aspect on the modeling of the Internet is the better results of the proposed model 
Parallel Addition and Rewiring Growth (PARG) compared to other models in the literature to capture specific 
features of the AS topology. In particular, the model is able to capture the dissortativeness of hubs (i.e., the 
most connected nodes in the network) and linear local assortativeness of nodes. The model proposed proves 
successful, however one must recognize that the AS connection are mainly logical and are due more to business 
and contractual agreements than physical interconnections, therefore missing part of the problem that is essential 
in the case of the Power Grid topology. 

Complex Network Analysis Applied to the Power Grid 

Complex Network Analysis works take into account the Power Grid at the High Voltage level usually to analyze 
the structure of the network without considering in detail the physical properties of the power lines. In our 
previous work |60j . we have analyzed several scientific papers that investigate Power Grid properties using 
Complex Network Analysis approach. There are two main categories: 1) understand the intrinsic property of 
Grid topologies and compare them to other types of networks assessing the existence of properties such as small- 
world or scale-free [SJ [HH [Ml E] ; 2) better understand the behavior of the network when failures occur (i.e., edge 
or node removal) and analyze the topological causes that bring to black-out spread and cascading failures of 
power lines [III El |24] . Few studies in the Complex Network Analysis landscape consider the possibility of using 
the insight gained through the analysis to help the design in terms of the reliability of the Power Grid. Examples 
are the study of the Italian High Voltage Grid ^5] and the study of improvement by line addition in Italian, 
French and Spanish Grids |57j- Also Holmgren [?D] uses the Complex Network Analysis to understand which 
Grid improvement strategies are most beneficial showing the different improvement of typical Complex Network 
Analysis metrics (e.g., path length, average degree, clustering coefficient, network connectivity) although in a 
very simple small graph (less than 10 nodes) when different edges and nodes are added to the network. Broader 
is the work of Mei et al. [5T] where a self-evolution process of the High Voltage Grid is studied with Complex 
Network Analysis methodologies. The model for Power Grid expansion considers an evolution of the network 
where power plants and substations are connected in a "local-world" topology through new transmission lines; 
overall the Power Grid reaches in its evolution the small- world topology after few-steps of the expansion process. 
Concerning the topic of Smart Grid with focus on the High Voltage network, Wang et al. [80l HI] study the Power 
Grid to understand the kind of communication system needed to support the decentralized control required by 
the the new Power Grid applying Complex Network Analysis techniques. The analyses aim at generating samples 
using random topologies based on uniform and Poisson probability distributions and a random topology with 
small-world network features. The simulation results are compared to the real samples of U.S. Power Grid and 
synthetic reference models belonging to the IEEE literature. These works also investigate the property of the 
physical impedance to assign to the generated Grid samples. Complex Network Analysis is not generally used as 
a design tool to propose adaptation strategies of current Medium and Low Voltage samples for the future Smart 
Grid as we use in this paper where we also assess the benefits in terms of economical improvement. 

Power Engineering Approach for Power Lines Design and Adaptation 

A power engineering problem that resembles our present proposal is that of Distribution Grid topological re- 
configuration. The early works on the topic of Merlin and Back |52j considered the problem with heuristic 
techniques and algorithms (e.g., branch and bound) to reconfigure the switches of the distribution infrastruc- 
ture. Later works of Cinvilar et al. |17j provide computationally tractable methods to modify the topology of 
the Distribution Grid by opening/closing switches to reduce losses. The authors provide a simplified formula 
to estimate the reduction in losses between the situation before and after the reconfiguration in the topology. 
Baran et al. [S] have also the idea of minimizing the losses and guaranteeing the load balancing. Basically, the 
problem deals with finding the minimal spanning tree, since the radial configuration needs to be preserved, that 
minimizes the objective function (i.e., system losses) while satisfying the constraints on voltage, capacity of lines 
and transformers, and reliability. For this problem, Baran et al. use a simplification of the power flow equations 
to compute the power flow in the network to be optimized. 

More recently, with the new interest in the Distribution Grid due to the Smart Grid and distributed generation 
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opportunities to be placed in the Distribution layer of the Power Grid, new interest has grown concerning the 
reconfiguration problem together with islanded Grids and micro-generation plants. The work of Ramesh et al. |63) 
focuses on the minimization of losses in the Distribution Grid. The authors provide several options that have 
been simulated and tested on the field to reduce losses on real and reactive power. The three solutions proposed 
concern three areas: distributed generation, capacitor placement and restructuring of feeders. The first proposal 
consist in placing local energy generation closer to the end users. In certain buses of the simulated IEEE 37 
Bus provides loss reduction of up to 9 MW. The second solution proposes the installation at optimal locations 
of capacitors. The third proposed technique concerns the restructuring of the topology of the network. This last 
aspect is closer to the topic of the current paper. The solution proposed in |63j is basically a reconfiguration of the 
opening and closing of the switches of the distribution network to maintain the network radial and satisfy its load 
requirements. This approach although valuable from the practical point of view, lacks the vision in imagining 
the future energy scenario and there is no investigation of which topology is best or according which principles 
the new pieces of infrastructure (e.g., cables) should be interconnected. The aspect of network reconfiguration 
of power networks to achieve a minimization of losses is the topic of [54] . The authors state that the only way 
to improve efficiency is by altering the topology. They focus on two aspects: (i) finding an optimal switching 
scheme of the switches connecting the lines and minimize losses, (ii) adding lines. The case study used is a 
IEEE 14 Bus where three more lines are added. The authors do not give any details why only three lines are 
considered and the motivation to choose specific nodes to attach the new lines. In addition, the only argument 
provided concerns the economic aspect of adding only the three lines, but no quantitative evidence is provided. 
The idea of intervening on topology to improve the efficiency of the Distribution Grid is as we consider in this 
paper, but we provide quantitative values, metrics and an economic analysis to decide which and how many lines 
need to be added in the Distribution Grid. We also consider another key ingredient that in the various works 
on reconfiguration is missing that is the added robustness that a more connected Grid may provide. In [TO] the 
problem of network reconfiguration in a Smart Grid environment is addressed. The idea of the authors is to 
reconfigure the topology of the network by operating switches in order to minimize the overload in the branches 
of the network. The scenario they consider takes also into account the higher penetration rate that distributed 
generation will have in the future. The numerical evaluation of the proposed genetic algorithm to reconfigure the 
network is realized on a simple 33 Bus network where the system configuration problem is solved to minimize the 
objective function (i.e., losses) and without violating the voltage and current constraints on the lines. Even in 
that work the importance of topology is claimed and the benefits in terms of reduced losses are shown. However, 
no additional lines or investment of the Distribution Grid are proposed and the benefit come just from the 
reconfiguration; no motivation or order of magnitude of the avoided investment are provided and the effects of 
topological changes available in the switching are anyway limited. The work of Xiaodan et al. [SS] is in the spirit 
of studying the reconfiguration in a micro-grid environment with micro-generation plant based on renewables. 
Basically, the problem is considered as two optimization problems: one related to the determining the capacity 
of each island of the Grid that has micro-generation capabilities, and the other is a problem of reconfiguration of 
the Distribution Grid with the objective of minimizing the power losses. The optimization algorithm is validated 
against two test Grids: the IEEE 33 Bus and the PG&E 69 node system. The reconfiguration for the system 
in an island and micro-generation uses techniques and approaches similar to traditional Distribution Systems. 
The problem that is considered is an optimization problem and the authors constrain the Distribution Grid to 
be always operated in a radial configuration. Usually, the reconfiguration problem boils down to the definition 
of an objective function to minimize the losses of the system and to establish constraints to satisfy the load. The 
function is then solved resorting to some heuristic. Our approach in the present work is novel: we do not deal 
with the congestion problem and our configuration of the network does not take into account how the switches 
are operated. In fact, in our dataset, we do not have the information of the state of the switches and we consider 
the full topology of the network for the Medium and Low Voltage samples, that is, the paths in the physical 
samples are already the best (i.e., shortest) achievable in the optimal operation of the switches (i.e., enabling all 
paths) . 

Power system engineers do resort to graph theory principles |7SJ US] , though GNA is not typically considered 
as a design tool. The traditional techniques involve the individuation of an objective function representing 
the cost of the power flow along a certain line which is then subject to physical and energy balance. This 
problem translates in an operation research problem. These models are applied both for the High Voltage 
planning [Ml and the Medium and Low Voltage [751 HSj since long time. Not only operation research, but 
also expert systems [SHI have been developed to help in the process of designing grounding stations based on 
physical requirements as well as heuristic approaches from engineering experience. The substation grounding 
issue is approached as an optimization problem of construction and conductor costs subject to the constraints 
of technical and safety parameters, its solution is investigated through a random walk search algorithm [36j . 
In [32| , a pragmatic approach using sensitivity analysis is applied to a linear model of load flow related to various 
overloading situations and a contingency analysis (N-1 and N-2 redundancy conditions) is performed with different 
grades of uncertainty in medium and long term scenarios. In practice, the planning and expansion problem is 
even more complex since it implies power plants, transmission lines, substations and Distribution Grid. In |37j 
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all these aspects are assessed separately and several challenges appear. For instance, in the planning of a High 
Voltage overhead transmission line, specific clearance code must be followed and load is not the only driver, but 
also topography and weather/climatic conditions (above all wind and ice) play an import role. For substation 
planning, the authors of [37] emphasize, in addition to the need for upgrading the Grid (e.g., load growth, system 
stability) and budgeting aspects, the multidisciplinary aspects which involve environmental, civil, electrical, and 
communications engineering. A more general approach proposed in |37| to deal with power system planning 
might be regarded as a multi-objective (e.g., economics, environment, feasibility, safety) decision problem thus 
requiring the tools typical of decision analysis [42] . 

Most of the works cited take into account mainly the High Voltage end of the Grid while not least important 
is the Distribution Grid especially in the vision of the future electrical system as proposed in this work where 
the end-user plays a vital role. The integrated planning of Medium and Low Voltage networks is tackled by 
Paiva et al. [61] who emphasize the need of considering the two networks together to obtain a sensible optimal 
planning. The problem is modeled as a mixed integer-linear programming one, considering an objective function 
for investment, maintenance, operation and losses costs that need to be minimized satisfying the constraints 
of energy balance and equipment physical limits. Even more challenges to the Electrical system planning are 
posed by the change in the energy landscape with several companies running different aspects of the business 
(generation, transmission, distribution). In addition, by accommodating more players in the wholesale market, 
transmission expansion should follow (as it is already for generation) a market based approach i.e., the demand 
forces of the market and its forecast should trigger the expansion of the Grid [14]. The same consideration 
regarding the need of a different approach in planning in a deregulated market are expressed in [69] where an 
optimization of an objective function in the market environment is applied. 

9 Conclusions 

The Smart Grid promises to revolutionize the energy sector providing the Grid with more efficiency, more real- 
time information and more renewable power sources. In such vision, the end-users will be able to produce their 
own energy and exchange the surplus in a completely free market. Based on the experience that we have accrued 
in our previous study of distribution Grid evolution considering Complex Network models, we have gone forward 
considering several types of possible evolution of networks by increasing connectivity. We have used physical 
samples from the Distribution Grid of the Northern Netherlands and we have applied six distinct strategies. 
We have performed two types of analysis of the results. First, a pure topological analysis considering a set of 
specifically defined metrics where random evolution, assortative high degree and least distance evolution strategies 
resulted as best satisfying the desiderata of the metrics. In our second set of analysis, we have also taken into 
account the physical properties of the cables, therefore realizing a weighted graph analysis, and considering the 
cost of evolving the networks assessing the cost of cables. In this weighted analysis, we have seen that the random 
and assortative high degree are the best strategies in a pure topological analysis, but they have costs that are 
extremely high and unrealistic to realize in practice. In addition, the weighted analysis gave us the possibility 
to investigate the elements of the Power Grid that influence the cost of electricity distribution (i.e.. Grid losses 
and Grid reliability) and the results of this analysis suggest that evolving the network by adding connections 
between the nodes with smallest distance is beneficial and provides in many cases better results compared to the 
other evolution strategies. Therefore, with an investment about 25% of the actual costs in cables already on the 
ground the Distribution Grid can improve consistently in reducing transportation costs. 

In this work we have for the first time realized a study on how to evolve the Distribution Grid by using the 
principles of Complex Network Analysis. We have studied the feasibility of this approach taking into account the 
costs and benefits for realizing different evolution approaches. This study sets the basis for a multi-flow study 
of the Distribution Grid with enhanced prosumer contribution to local energy generation. Moreover, it is a step 
forward in realizing a decision support system tailored for the Smart Grid. Our future efforts will be devoted to 
the further development of this decision support system. 
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Appendix A Network topological requirements 



In |59| we proposed a number of metrics useful for analyzing synthetic network topologies for the Power Grid 
having in mind decentralized energy trading. We recall them here since we use them to measure the adequacy 
of the proposed evolution strategies of the physical samples of the Dutch Distribution Grid. 

Quantitative requirements 

As a global statistical tool, quantitative requirements are even more useful as they give a precise indication of 
network properties. Here are the relevant ones when considering efhciency, resilience and robustness of a power 
system. 

• Characteristic Path Length ( CPL) lower or equal to the natural logarithm of order of graph: CPL < ln{N). 
This requirement represents having a general limited path when moving from one node to another. In the 
Grid this provides for a network with limited losses in the paths used to transfer energy from one node to 
another. 

• Clustering Coefficient ( CC) which is 5 times higher than a corresponding random graph with same order 
and size: CC > 5 x CCrc- Watts and Strogatz [53] show that small- world networks have clustering 
coefficient such that CC CCrq- Here we require a similar condition, although less strong by putting a 
constant value of 5. This requirement is proposed in order to guarantee a local clustering among nodes since 
it is more likely that energy exchanges occur at a very local scale (e.g., neighborhood) when small-scale 
distributed energy resources are highly implemented. 

• Betweenness-related requirements: 

— A low value for average betweenness in terms of order of the graph v = jj, where a is the average 
betweenness of the graph and N is the order of the graph. For the Internet Vazquez et al. [75] have 
found for this metric v sa 2.5. Internet has proved successful to tolerate failures and attacks [T51[3], 
therefore we require a similar value for this metric for the future Grid. 

— A coefficient of variation for betweenness Cy = = < 1 where s is the sample standard deviation and x 
is the sample mean of betweenness. Usually distribution with Cy < 1 are known as low-variance ones. 

The above two requirements are generally considered to provide network resilience by limiting the number 
of critical nodes that have a high number of minimal paths traversing them. These properties provide 
distributions of shortest paths which are more uniform among all nodes. 

• An index for robustness such that Rob^ > 0.45. Robustness is evaluated with a random removal strategy 
and a node degree-based removal strategy by computing the average of the order of the maximal connected 
component (MCC) of the graph between the two situations when the 20% of the nodes of the original graph 
are removed [55]. It can be written as Robiv = . Such a requirement is 
about double the value observed for current Medium Voltage and 33% more for Low Voltage samples [55] . 

• A measure of the cost related to the redundancy of paths available in the network: APLiQth < 2 x CPL. 
With this metric we consider the cost of having redundant paths available between nodes. In particular, 
we evaluate the 10*'' shortest path (i.e., the shortest path when the nine best ones are not considered) by 
covering a random sample of the nodes in the network (40% of the nodes whose half represents source nodes 
and the other half represents destination nodes). The values for the paths considered are then averaged. 
In the case where there are less than ten paths available, the worst case path between the two nodes is 
considered. This last condition gives not completely significant values when applied to networks with small 
connectivity (i.e., absence of redundant paths). 



Metric 


Efficiency 


Resilience 


Robustness 


CPL 


/ 






CC 


/ 






Avg. Betweenness 




/ 




Betw. Coeff. of Variation 




/ 










/ 


APL^Qth 


/ 


/ 





Table 37: Metrics classification related to properties delivered to the network. 
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The above quantitative metrics can be categorized into three macro categories with respect to how they affect a 
Power Grid: efficiency in the transfer of energy, resilience in providing alternative path if part of the network is 
compromised/congested and robustness to failures for network connectivity. Table 37 summarized the property 
each metric assesses. 
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